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Abstract 

We investigate finite-size corrections to anomalous dimensions of large-spin twist-two operators in 
the planar maximally supersymmetric Yang-Mills theory. We develop a framework for analysis of these 
corrections, that is complementary to the conventional spin-chain approach, by making use of the hole 
rather than the magnon picture. From the dual string theory perspective where the large-spin operator is 
identified with the Gubser-Klebanov-Polyakov (GKP) string, our approach is equivalent to constructing 
the first Liischer correction to the energy of the GKP string by incorporating the contribution of virtual 
excitations propagating on it. It allows us to propose a formula that controls a particular class of 
large-spin corrections to the twist-two anomalous dimension and holds at any value of the coupling 
constant. Compared to wrapping corrections computed with magnons propagating on the spin chain, 
the finite-size corrections that are encoded in our formalism start at a lower loop level. Our formalism 
thus calls for modification of the asymptotic contributions which are conventionally incorporated within 
the Asymptotic Bethe Ansatz. An educated guess allows us to remedy this pitfall and successfully 
confront our predictions with known results up to five loop accuracy at weak coupling. Finally, our 
formula sheds light on the weak-to-strong coupling transition for the subleading large-spin corrections 
under study and confirms stringy expectations at strong coupling where they are found to be identical 
to the first Liischer correction to the vacuum energy of the 0(6) sigma model. 
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1 Introduction and discussion 



AdS/CFT correspondence [T] provides a powerful framework for geometrical reformulation of gauge 
theory dynamics in the planar limit in terms of free strings propagating on a nontrivial curved back- 
ground. For maximally supersymmetric Yang-Mills (SYM) theory it allowed one to establish a precise 
dictionary between field and string theory observables, in particular, with the spectrum of scaling di- 
mensions of single trace gauge invariant operators in J\f = 4 Yang- Mills theory identified with energies of 
string configurations on AdSsxS 5 target space. The weak/strong nature of the duality was hampering 
direct tests of the conjecture since the weak coupling regime on the gauge side corresponded to a highly 
quantum regime for the string sigma model, and inversely the domain when quantum fluctuation on 
the strings are suppressed required techniques for strong coupling expansion of anomalous dimensions 
of composite operators. The difficulty in pursuing the endeavor of testing the duality as a function of 
parameters defining both sides of the story was overcome by identifying hidden integrable structures 
of the system. On the string theory side, the classical sigma model was found to possess an infinite 
number of conserved charges [2]. On the gauge theory side, the dilatation operator whose eigenvalues 
define anomalous dimensions was found to correspond to a certain long-range magnet [31 SI E] , — a 
generalization of an earlier discovery of the XXX Heisenberg spin chain for maximal helicity operators 
at one loop in QCD [6, 7J. 

In spite of this revelation the extraction of exact anomalous dimensions at any value of 't Hooft 
coupling g 2 = gy M N c /(4:7r) 2 from the underlying integrable models is still highly nontrivial. However 
since the energies depend on quantum numbers of corresponding states, like spin, one possesses an extra 
lever arm which makes the calculations manageable in certain limits. Namely, significant simplifications 
of underlying calculations occur when one considers operators with large values of spin on either compact 
or non-compact parts of the AdSsxS 5 background, or both 0[9]. Thus one is on the quest of finding 
systematic techniques which provide interpolating formulas between the regimes of weak and strong 
coupling. Being interested in generic properties intrinsic to supersymmetric gauge theories and more 
importantly not, since isometries of the noncompact manifold corresponds to space-time symmetries 
shared by both types of field theories, we will focus on the latter in the present work. This configuration 
corresponds to the so-called Gubser-Klebanov-Polyakov (GKP) string [9], which is a classical folded 
closed string rotating with a large angular momentum S on an AdS3 part of the AdSsxS 5 space. In 
this limit, the string of proper length R = 2 InS + O(l) becomes homogeneous and stretches all the 
way to the boundary when S — > oo. Then its proper energy, as measured in the rotating frame of the 
string, scales with the length R, 



where the energy density w(g) is determined classically by the string tension g, w(g) = 2g + . . . . This 
scaling behavior is akin the large-spin asymptotics of scaling dimensions of Wilson operators in all gauge 
theories. While the GKP string determines the leading Regge trajectory of the string with a given spin 
S, its gauge theory dual corresponds to the twist-two operators, schematically, 



where Z is one of the complex scalars of N = 4 SYM and D + is the light-cone projection of the covariant 
derivative. The large-spin asymptotics of their scaling dimension A is driven by the cusp anomalous 
dimension HOI 1 ! 



This behavior is known to persist to all orders in 't Hooft coupling constant and reflects the copious 
emission of soft gluons at the edge of phase space [I2j[l3]. Thus, the identification |GKP string) = 02,s, 

x The cusp anomalous dimension A(g) is a function of the coupling constant g only that controls the UV divergences of 
a cusped Wilson loop [TT| with a large hyperbolic angle [10], i.e., in the light-like limit. 



E'gkp = E GKP (g, S) = w(g)R + ... . 




2 ,s = trZD^Z + ... , 




E = A - S = 2A(g) In S + ... . 
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and its corollary -Egkp = E, provided a stringent test for the conjectured duality [9]. The scaling law 
(|l,ip was shown to hold at higher loop orders within the 1/g expansion of the string sigma model, with 
explicit computations performed at one |14| and two loops |15j . 

The interpolating formula arising from the weak coupling expansion on the gauge theory side was 
deduced by studying the large-magnon limit of the long-range spin chain. In the language of the 
spin chain encoding the integrable structure of the dilatation operator in N = 4 SYM, this scaling 
dimension corresponds to the energy of the length L = 2 chain with S magnon excitations of momenta 
Pj propagating on it. The integrability implies that the total spin-chain energy is a sum of individual 
excitations with dispersion relation [16] E(p), 

S I 

A = 2 + Y, EijPj) , E(p) = V 1 + !% 2 sin 2 (;p/2) , (1.4) 

while the momenta obey Asymptotic Bethe Ansatz (ABA) equations [5] 

e iLpk = Y[S( Pj , Pk ), (1.5) 

where S(j>j,pk) is the two-particle scattering matrix. These considerations yielded an integral equation, 
the so-called Beisert-Eden-Staudacher (BES) equation [IT], for the density of Bethe roots 

Uj = \ cot {p 3 /2)yfl + 16g 2 sm 2 ( Pj /2) (1.6) 

which generates the cusp anomalous dimension. The solution at weak coupling [T7j matched the result 
of tedious calculation to fourth order in g 2 [13l HEJ HSJ [20j [21], while much trickier strong-coupling 
analyses [22[ [2"3"[ [24"[ [25] were found in agreement with the energy density w(g) to two loop order in 
1/g expansion. Together with the numerical interpolation of Ref. [26], this confirmed the conjectured 
duality for non-protected observables, yielding w(g) = A(g). 

This success in interpolating between the weak and strong coupling for the cusp anomalous dimension 
raises the question of how far we can go in unravelling the transition between the gauge and string theory 
at the level of subleading corrections in large spin. It has been known already for some time that the 
approach developed for the cusp anomalous dimension is amenable to generalizations [27] for the first 
subleading contribution ~ log S controlled by the so-called virtual scaling function, B = 2 + 0(g 4 ). 
The twist-two energy is then found to assume the form 

E = 2AlogS + B + ... , (1.7) 

where log S = log <S + 7e> with 7e the Euler-Mascheroni constant, and where both A and B are functions 
of the coupling only, interpolating between the gauge and string theory predictions. The ellipsis in (II. 7j) 
stands for the corrections we are after in this paper that are explicitly suppressed at large spin. Both at 
weak and strong coupling, our understanding of their structure is mostly based on explicit expressions for 
the first few orders of corresponding expansions. At weak coupling, for instance, the energy receives an 
infinite series of corrections suppressed by inverse powers of the spin, each possibly receiving logarithmic 
enhancements from powers of log S. The details of this expansion are not entirely arbitrary and display 
an interesting pattern [28]. In particular, leaving aside the corrections ~ 1/5, which can be eliminated 
by an appropriate redefinition of the expansion parameter S —> rj(S) = S + . . . , to which we shall 
come back later, the first new, 'dynamically independent' contributions appear at order 0{1/S 2 ). One 
can say to some extent that the typical scale of corrections at weak coupling is ~ 1/S 2 . A similar 
pattern was observed at strong coupling [29J, at least classically, suggesting that a direct interpolation 
along the lines as explained before is possible. However, the fact that this program cannot be realized 
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becomes transparent at the one-loop level on the string side. Here, the situation changes drastically 
with corrections of the type ~ 1/ log S prevailing over the aforementioned power suppression. This was 
uncovered both from direct string theory computations [30 |, I3 H [32] and integrability methods [33]. To 
one-loop order, the first stringy corrections to the scaling (|1.7|) is found to be driven by the Casimir 
energy [3Q1 ED M M 

5E~--^-, (1.8) 
12 log b 

that merely reflects the contributions of the five massless modes of the GKP string. The reconcilia- 
tion of the string and gauge theory predictions thus entails considering an interpolating function that 
dependents in an entangled manner on both the 't Hooft coupling and the spin. 

The weak-to-strong transition a la BES seems furthermore doomed from the outset as soon as one 
takes a first glance on the above spin suppressed contributions. The reason is that this strategy relies 
solely on the ABA equations, which, though defined for any value of the coupling, are known not to 
account entirely for all gauge and/or string contributions |34[ 135], This implies, in particular, that the 
identities quoted in (jl.4p and (|1.5p are not correct beyond certain accuracy which is determined by the 
size of the so-called wrapping corrections!! The exact energy should then be written as 

E = E ABA + £ wrap , (1.9) 

with .E ABA , only, being determined by (jl.4l 1 1 . 5 [) . To account for the new contribution, i.e., E WTap , 
one has to start from the integrable model on the string side and invoke the Thermodynamic Bethe 
Ansatz [36]. On this side, indeed, the whole problem is formulated within the framework of a two- 
dimensional quantum field theory, compactified on a cylinder of circumference ~ L, and whose ground 
state, the Berenstein-Maldacena-Nastase (BMN) reference state \8\, is identified with the spin-chain 
vacuum tv Z L . In this set up, the decomposition (jl.9p can be understood qualitatively as the splitting 
of the energy of an excited state into two pieces: one which, under the assumption of integrability, 
simply sums up the mechanical energies of each particle constituents and the other one incorporating 
the energy stored in the virtual fluctuations that they source. The latter contributions are conveniently 
understood as originating from the exchange of virtual excitations circulating around the cylinder |34] 
and scattering with the real particles that build up the state under study. An ingenious trick [36J, that 
goes under the name of double Wick rotation and which interchanges space and time, allows one to 
analyze the dynamics of these virtual excitations which mirrors in many aspects the original one. At 
this end, these considerations can be used to advocate the leading wrapping contributions, also known 
as the (first) Liischer correction, to the energy of the twist-two operators |37[ 138] . The energy E WTap is 
then found to admit the form 

E™*P~ [ dp e- LE *W (■•■), (1.10) 



exhibiting a typical exponential decay with the system size L and its rate controlled by the energy |34] 
E* = 2 arcsinh w 1 + pi/ (Ag) of the particle exchanged in the mirror kinematics. For our discussion, it 
matters that the wrapping corrections appear delayed at large spin: they are expected indeed not to 
affect the validity of the leading asymptotics (|1.7p such that the coefficients A and B are free from wrap- 
ping and are governed entirely by the ABA equations. This is confirmed by the large spin asymptotics 
of the leading wrapping correction for the twist-two operators [38] that is found to scale like [39J 

wran S ( , 64vT 2 b , 



2 The terminology comes from the fact that in the spin-chain language these contributions are associated to wrapping 
interactions whose range exceeds the length of the chain. 
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It immediately confirms however that the contributions ~ 1/S 2 that we want to elucidate in our current 
analysis are contaminated by the wrapping effects. The same wrapping corrections appear also crucial 
at strong coupling for the recovery of the stringy expectation (|1.8p as was shown in |33j . 

The computation of this wrapping energy is amenable to a systematic treatment proposed in [4U\ I3T] . 
But this framework comes with a hefty pricetag: one has to deal with an excited state with a very large 
number of quanta, namely, S'> 1 magnons accounting for the light-cone derivatives in (11. 2\ . This makes 
the analysis of the wrapping energy rather difficult, and, so far, only the first two leading contributions, 
i.e., the four- [38] and five- loop [42] expressions, have been extracted]! Though they possess a very 
remarkable feature of being applicable at any spin, not necessarily large. 

The AdS/CFT correspondence challenges us however to unravel the details of the transition between 
the gauge and string theory and already hints at a possible new path. Namely, getting back to the 
string theory result (jl.8p and its interpretation as a Casimir energy, one immediately realizes that the 
contributions we are after are naturally associated with the finite-size corrections induced by virtual 
particles now circulating around the GKP string itself, and thus over a length R ~ 2 log S. This shift 
in paradigm, which places the GKP string at the center of the analysis and treats it as a vacuum on its 
own, yields at a qualitative level an immediate realization of the change of the scale 1/S 2 to l/logS" as 
one passes from the gauge to string expectations. The main observation is the remarkable feature of the 
dynamics on top of the GKP string, which was elucidated by Alday and Maldacena |43| : the theory is 
gapped at weak coupling with a mass gap m ~ 1 but becomes asymptotically gapless at strong coupling 
with m ~ exp (—irg). At strong enough coupling, this phenomenon is moreover governed by the 0(6) 
sigma model which stands for the low-energy effective theory on the GKP vacuum. A direct analogy 
with a relativistic QFT would then lead us to conclude that the two previous scales merely reflect these 
two extreme regimes, with the finite-size corrections estimated as ~ exp (—mR) ~ 1/S 2 in the massive 
'phase' and as ~ 1/R ~ 1/logS* in the massless one. The validity of this analogy is far from being 
obvious and actually relies on the double Wick rotation invariance of the GKP background, which was 
first articulated in [33] . 

The purpose of this paper is to establish precisely, at a quantitative level, the above picture by 
advocating and constructing a representation for the leading finite-size correction which transparently 
incorporates the physics of excitations on top of the GKP string. It will imply that the loop-by-loop 
order large-spin predictions on the gauge and string theory side are the two ends of the stick emerging 
under the appropriate conditions, i.e., g 2 <C l/logS" and g » ^ log log S, respectively. Our first step, 
however, is to set up the following ansatz for the complete energy of the string. Viewed as a vacuum 
energy, we aim at decomposing it as 

E = E hulk + E FS , (1.12) 

where the finite-size energy E FS is always suppressed with the system's length, in the large volume 
limit log S 3> 1. The details of this suppression depend, however, on certain features of the theory, 
in particular, on its spectrum. The bulk energy, on the other hand, scales with the length and incor- 
porates therefore the leading asymptotic behavior (jl.7p . From the world-sheet point of view, the bulk 
energy should resum corrections associated with diagram of trivial topology contrary to the finite-size 
corrections which accommodate those with nontrivial windings Q They are shown in Fig. Q](d) and (c), 
respectively. The decomposition (11. 12ft immediately raises several questions: How do we compute the 
bulk/finite-size energy and what are their relation to the AB A/wrapping expression? Probably the most 

naive and conservative attempt is to set that E hulk = _e; aba an d E FS = E wrap . We shall see that this is 

3 This approach has indeed the drawback of involving the summation over an infinite number of virtual excitations, bound 
states of elementary ones, which scatter with an 'equally large' number of real particles. The five-loops expression [42] 
further requires to into account the incorporation of the back-reaction of these virtual contributions on the quantization 
conditions for the set of momenta {T3}. 

4 The former contribution can also be understood as renormalizing the GKP background. 
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not possible however. To realize why this is not the case, let us focus on the finite-size corrections and 
derive a Liischer-like formula that encodes them. 

As we already said, our approach relies on a complimentary formalism which operates with a different 
set of degrees of freedom than the magnons. Information about these dual excitations can be deduced 
directly from the ABA equations, up to corrections that are suppressed by at least the system size 
~ logS. Integrable models, indeed, are often endowed with a particle-hole transformation which allows 
one to exchange particle and vacuum modes. While the hole degrees of freedom are invisible within 
the ABA approach they become manifest within the formalism of the Baxter equation [45 } [46] . The 
long-range Baxter equation is a second-order finite-difference equation for the Q-function 

A + {u + i)Q(u + i) + A_(u - i)Q(« - i) = t(u)Q(u) . (1.13) 

where A±(u) are known dressing factors and t{u) as the eigenvalue of an auxiliary transfer matrix. The 
solution to (|1.13p identical to ABA is sought in the polynomial form 

s 

Q(u) = Y[(u- Uj ), (1.14) 

i=i 

where the set of rapidities Uj coincide with the Bethe roots (|1.6|) . The equivalence between the two 
formalisms is achieved by setting u = Uj in both sides of Eq. (j!.13j) and observing that the latter admits 
the form 

, , A_(u-±)Q(u-i) 

1 + 3^=0, Yh (u) = . , 1.15 

A + (u + ^)Q[u + i) 

where Y^u), introduced here for convenience, is a quantity that we shall encounter on many different 
occasions in the following discussion. An easy inspection confirms its agreement with Eq. ([2]), as 
anticipated. 

However, now in addition to the Bethe roots, the Baxter equation encodes another set of data that 
enters through the transfer matrix. Namely, at one-loop order t(u) = t^°\u) + 0(g 2 ) is parametrized 
by its roots 

L 

*<°>(u)= 2 !](«-**)• (1-16) 

k=l 

The idea behind the particle-hole transformation is to trade the dynamics of the Bethe roots Uj for the 
one of the holes 5k- In the large spin limit, two of the holes are large and somehow decouple, leaving L — 2 
(small) holes identified with the dynamical degrees of freedom [47] . These holes encode the dynamics 
of excitations, here, the scalar Z-fields propagating on the background of derivatives. 

The state with no small holes corresponds to the twist-two operator and is identified with the GKP 
vacuum, see Fig. [T] (a)H Adding a hole amounts to inserting a Z field in this background. This is not 
the only flavor of excitations however. To find all of them one should extend consideration to the full 
system of ABA equations [5]. The whole set of elementary excitations was elucidated in Refs. [48 \ 149]: 
they are associated to the twist-one partons of the gauge theory, i.e., 6 scalars, 4 + 4 fermions and 1 + 1 
gauge fields! For illustration, a one-particle excitation X on top of the ground state is found as 

s 

|one-particle excitation) = ^ c s tr ZD S + XD^T S Z + . . . , (1.17) 

s=0 



5 The SO(6) symmetry of the GKP vacuum, which is not obvious for the operator (II. 2[) . is achieved by means of isotopic 
degrees of freedom in the form of composite Bethe roots which do carry zero energy and momentum [48] . 



They carry so (6) ~ su(4) isotopic degrees of freedom which are easily inferred from their corresponding degeneracy. 
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(a) (b) (c) (d) 



Figure 1: (a) Vacuum corresponds to twist-two operators, (b) Real one-particle excitation represents 
twist-three operators, (c) Corrections from finite size of the system, (d) Virtual bulk correction renor- 
malizing the background. 

where X is one of the previously mentioned parton inserted in the twist- two background (II. 2p . This 
is a twist-three operator, which on the string theory side corresponds to a particle propagating on the 
closed-string world sheet, see Fig. [1] (b). The generalization to higher twists follows the same line of 
reasoning with possible additional states in the form of twist > 1 partons which are bound states of 
gauge fields. This pattern of excitations and their associated dispersion relations are in good agreement 
with string theory predictions |14 4 150} [5T]. though a few mismatches still await their reconciliation. 

Equipped with the spectrum of excitations, one can make an educated guess for the form of the 
leading finite-size corrections. This one can be argued to take a universal form, known as the Liischer 
formula. The logic here is identical to the one for magnons |34^ 137] . Our proposal for the Liischer 
formula reads 

^ FS = " E £ / * y * mirr ° r (p) + • • • » (1-18) 

where the function l^ mirror (p) controls the amplitude for propagation of a ^-particle around the GKP 
string with the momentum p in the mirror kinematics. The ellipsis correspond to higher Liischer 
corrections associated with multiparticle exchanges and n* is the number of particles of a given ^-flavor, 
e.g., rih = 6 for scalar, etc. We shall find that the double Wick rotation symmetry is valid for all 
excitations and thus makes the analysis akin to the one typically performed in relativistic theories, 

yr rror (p) ~ e-^ p ) — ► E FS ~Y, [dpe- RE * ip) , (1.19) 

where E±(p) is given by the dispersion relation of a real excitation and R ~ 2 log S. This result 
immediately implies that the contribution of a ^-particle to the energy f 1 1 . 1 8 1) scales as 

E™ ~ 1/S 2m * , (1.20) 

with m* = E*(p = 0) being the mass of the excitation. As a consequence, within our accuracy at weak 
coupling, i.e., ~ l/S 2 , we can restrict the summation in Eq. (11 . 18f) to the set of twist-one excitations 
only, i.e., * = h (hole = scalar), f (fermion), gf (gauge field). This also explains why at strong coupling 
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only the scalar sector, which is controlled by the 0(6) sigma model, remains relevant: while all other 
excitations maintain their masses of order one |14[ I43j. the holes' mass becomes exponentially small. 
These conclusions survive a more careful analysis given below, but it is interesting to point out that 
this consideration requires a modification of (|1.19|) by an extra ~ R° ~ log S contribution, 

RE+ip) -> 2E*(p) log S + 25E* . (1.21) 

This seems to reflects a curious feature of the propagation of excitations in the GKP background, 
probably related to the loss of its homogeneity, observed explicitly at strong coupling [14, 32J. At 
weak coupling, this new factor determines the dependence of the finite-size corrections on the coupling 
constant, which according to our analysis conducted below is 5E+ ~ —logg 2 resulting in E FS ~ g A . 
Supersymmetry somewhat delays the onset of the finite size corrections by one order of the coupling, 
however, and as a consequence one finds E FS ~ g 6 Hence, compared to the wrapping effects introduced 
in the magnon picture the finite size corrections arise starting already from three loop order for twist-two 
operators. This immediately exhibits the fundamental difference of our finite size correction compared 
to the conventional Liischer formula (ll.lOp that one add to the ABA energy. Therefore, blindly adding 
our (|1.19p to the ABA energy would not produce the correct energy. This is not totally surprising 
since after all, the Liischer formula fjl . 191) that we are proposing does not quite sum up effects due to 
propagation of virtual magnons, as it is the case for the wrapping energy (jl.lOp . but instead of holes 
which, though somewhat dual to them, are nevertheless distinct. 

The Liischer formula was given schematically in (|l,18p . An important part of our analysis will be 
therefore to make it work at a quantitative level. To this end we shall analyze the transformation to the 
mirror kinematics separately for all our excitations and construct the relevant Y functions. We shall see 
that they can all be expressed in terms of the solution to the BES equation and its siblings [27] at the 
given accuracy in the spin. These can be though as fixing the asymptotic data, directly constructible 
from the ABA equations. This has a major benefit that the formula (|1.18p can be univocally determined 
at any value of the coupling. It will then be possible to construct a weak coupling expansion of (|1.18p 
to high orders of perturbation theory and find, in particular, that it admits the expansion 

^ FS = ^E^ 2n +°( 1 ^ 2 ), (1-22) 

^ n>2 

where Vm is a coupling independent polynomial in the variable log S = log S + 7e or degree n — 2. 
It thus predicts that the higher-loop corrections are dressed by higher powers of log S with no bound 
as the loop order n becomes infinite. This can be seen as the first evidence that the formula (|1.18p is 
capable to reconcile the gauge and string theory predictions through its intricate dependence both on 
the coupling g and the length ~ log S. 

Our proposal for the twist-two energy would remain incomplete without full control on the second 
halve of the equation (|1.12p given by the bulk energy. It is much more difficult to advocate its general 
expression beyond the fact that it should incorporate the leading asymptotics (]1.7p . Our approach will 
remain unfortunately empirical and is advocated by the comparison to the ABA predictions as follows. 
The ABA equations, or equivalently the asymptotic Baxter equation, can be solved at any coupling and 
to high accuracy at large spin. A thorough analysis demonstrates that the energy stemming from ABA 
can be decomposed in terms of sign preserving (regular) _E reg and alternating E 3 ^ contributions 

E ABA = E Tes + E alt , (1.23) 

7 The cancellation between bosonic and fermionic contributions may not be as effective if one considers ~ 1/S 4 and 
higher suppressed corrections. Bound states that contribute at this order do not have super-partners, for instance. It 
makes us believe that the onset of the finite-size corrections starts at two loops, though not apparent at the level of the 
~ 1/S 2 contribution. 
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with S alt oc e , which depends on the momentum t? of the state. The latter expression evaluates to 
(— l) s on twist-two states with integer spin. This requires including into our consideration states with 
nonzero momentum ■& = ir, realized for odd (unphysical) spins S. This continuation was surreptitiously 
performed in the finite-size energy above, see Eq. (pQ2l) . where the factor = (-l) s can be thought as 
twisting the boundary conditions for excitations around the GKP stringH The physical configurations 
corresponding to even values of the spin are then associated with the periodic boundary conditions. 
A direct analysis shows that the quantity £' alt is suppressed by two powers of the spin at large S, to 
any order at weak coupling. The general structure of its expansion is furthermore akin to the finite- 
size energy: both are given by a series of terms with increasing logarithmic enhancement in S as one 
passes to higher and higher orders in the coupling. Though their general forms appear to be similar, 
the accompanying expansion coefficients are different. Nevertheless, the former smilarity, together with 
the expectation that the Liischer formula (|3.92[) should complement in some way the ABA prediction, 
inspires us to promote the substitution rule 

E aXt _^ e fs ^ ( L24 ) 

as a mean of exacting the ABA prediction. It is equivalent to the statement 

E hulk = E res , (1.25) 

which together with (|1.18p determines the (full) energy (| 1 . 1 2 1) . This, of course, should be understood as 
an empirical rule at the currently considered accuracy in large spin. Some physical evidence for it will be 
presented shortly. However, before turning to it, we would like to stress that by means of Eq. (|1.24p our 
proposal for exacting of the ABA prediction with the help of the Liischer formula becomes complete: at 
the given accuracy in the spin, the bulk energy can be obtained from the ABA equations and these same 
equations can be used to construct all ingredients involved in the Liischer formula fll . 18j) . It is found, for 
instance, that in distinction with the alternating or finite-size energy the large spin expansion of the bulk 
energy is much more rigid. With the help of an appropriate expansion parameter r] = i](S) = S + . . . 
alluded to before, it takes the form 

£ bulk = 2A (log fj - +B + C/r? + o(l/r, 2 ) , (1.26) 

where C is a new coefficient function of the coupling only. After re-expansion in the spin, the bulk 
energy displays a logarithmic enhancement at order ~ 1/S 2 but this one is bounded by log 2 S for any 
coupling, contrary to the alternating or finite-size energy. 

Perhaps more importantly, the rule fjl.24f) makes it possible to compare the content of (II. 18ft with 
existent results for the large spin expansion of twist-two anomalous dimension up to five loops, including 
wrapping corrections. These results based on the Liischer formula 'on the BMN vacuum' are built on 
a much firmer basis than ours and it is thus more than a requirement to successfully compare the two. 
To this end, one merely observes that 

E FS _ E M + ^wrap ( L27 ) 

holds as a corollary of our previous discussion. This should be the case at order 0(1/ S 2 ) at large 
spin and to any order at weak coupling. After extracting the relevant large spin expression for E wra,p 
from [381 [391 02], we will be able to verify Eq. (jl.27p up to five loopsj^l We notice finally that, through 
the identity (jl.27j) and the common structure of the finite-size and alternating ABA energy at large spin, 



3 The linearity in e 11 ' of the first Liischer correction originates directly from the linearity in the Y functions. 
5 This agreement is actually observed only for the physical configuration, that is, for even spin. 
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our proposal allows us to make immediate qualitative predictions at higher loops, notably, suggesting 
that 

K-L P ~ -^T^ ■ (1-28) 

This discussion sets up the stage of our analysis. Before plunging into detailed computations of the 
bulk and finite-size energy, let us come back to our problem and let us see how our proposal helps us 
uncover the weak-to-strong transition in coupling. 

To some extent, the explanation of the weak-to-strong transition is already embodied in the ABA 
prediction. To our accuracy in the spin, the main role is played by the alternating component of the 
energy that displays the enhancement of the logarithmic growth as the loop order increases. This 
behavior is shared by the Liischer formula (jl . 18H . which also motivated the substitution rule (11.241) . 
What happens next is that the weak coupling expansion breaks down as soon as g 2 ~ 1/logS*, calling 
for a resummation. All of these corrections can be properly encoded into the formula 

£ alt = I j> d Ph (u - i/2)YC" m (u - §) + . . . , (1.29) 

with ellipsis standing for corrections suppressed stronger than ~ 1/S 2 at weak coupling. This identity 
looks pretty similar to the Liischer formula (11.18P if not for a different choice of integration contour 
and the absence of contributions from other GKP excitations than the holes. The similarity is enough 
however to understand the behavior of E 3 ^ at strong coupling. Namely, when the coupling gets large, the 
energy (|1.29h gets higher due to the reduction of the mass m = of the scalar, since y^ rLirror ~ \/S h . 
It then takes a very suggestive form 

E* lt = -—{-l) s K 1 {mR) + ... , (1.30) 

7T 



where K\(z) ~ exp(— z) for z S> 1 is the modified Bessel function, m ~ exp(— ng) the mass gap |43j 
and R ~ 2 log S the length. For even spin, this equation gives the contribution to the vacuum energy 
of two massive bosons subject to periodic boundary conditions on a cylinder of size R 3> 1/m. For 
comparison, the finite-size energy (|1.18p . which is also dominated by the contribution from the lightest 
modes, i.e., 6 scalars, at strong coupling, reads 

E FS = - — (-l) s K 1 (mR) + ... . (1.31) 

It suggests that the ABA equation accounts only for 2 out of 6 scalars excitations. The former may be 
viewed as the two scalar fields Z, Z which are not treated as independent dynamical excitations in the 
BMN picture, since Z is part of the background and Z is a BMN mass 2 excitation that decays into two 
BMN scalara 10 !. In light of it, the proposed substitution I? alt — > E FS can be understood as correcting for 
the proper number of scalar excitations not captured by ABA. Adding further contributions of fermions 
and other bosons, the formula (|l,18p together with (jl.24p appears as a natural attempt for exacting the 
quntitative predictions for the energy of the GKP string. 

The results (|1.3ip is transparent enough to allow one to guess the final step towards the string per- 
turbative result. It reveals simultaneously the limitation of our expression at strong coupling. Namely, 
as was already said earlier, the result (|1 . 18[) accounts only for the first Liischer correction to the vacuum 
energy computed within the 0(6) sigma model. The latter however predicts an infinite number of higher 
Liischer corrections, with all of them encrypted into a system of TBA equations, see, e.g., [52J for a 
recent discussion. From this more established approach, it becomes obvious that the asymptotics (jl.3ip 
can only prevail provided mR S> 1. This requires an extraordinarily large length that the perturbative 



10 This is the traditional interpretation of the mixing ZZ ~ YY ~ XX between the various scalars in the BMN 
background. 
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string analysis cannot afford. The latter always operates in a regime where the Compton wavelength 
~ 1/m of the 0(6) particles is much bigger than the length R, i.e., mR <C 1. In this regime, the natural 
excitations appear to be 5 massless Goldstone bosons, instead. The stringy perturbative prediction 
necessarily resums an infinite tower of Liischer corrections, which from the weak-coupling standpoint all 
originate from highly subleading contributions suppressed as ~ 1/S 2n with n > 1. Tracing back these 
corrections would constitute a considerable improvement of our approach and would entail working out 
a system of TBA equations for the dynamics of excitations on top of the GKP string. This is however 
beyond the scope of our present analysis. 

A reader would notice that our previous discussion is, of course, a well understood phenomenon in 
relativistic quantum field theory, which can be studied in great details corresponding system becomes 
integrable. The result for the finite size energy can be encapsulated into the formula [531 



£ FS = -^ + ..., (1.32) 

where the omitted terms stand for corrections ~ 1/R d with d > 2, associated with irrelevant defor- 
mations of the 0(6) sigma model, and 'UV contributions ~ 1/S 2 from fermions and gauge fields. 
Here c(mR) is the effective central charge of the 0(6) sigma model. It interpolates between c(mR) ~ 
exp (— mR) when mR S> 1 (massive regime) and c(mR) ~ 5 when mR <C 1 (massless regime). In 
between, it is controlled by the full set of TBA equations for the vacuum state of the 0(6) sigma model. 
Our expression (11.180 is only applicable in the massive regime while the perturbative string analysis 
works in the UV regime. In this latter domain, which is far beyond our reach, the 0(6) sigma model 
can be used to shed some light on the structure of the large spin expansion. It is known, for instance, 
that in the 0(n) sigma model the effective central charge can be found to take the form |52| [55j^l 

. . . 3(n — 1) „ /log log (1 /mR) \ 

c{mR < 1) = (n- 1 - - V — , + O S 9 S V / . 1.33 

v ' y 1 2 log (I/mi?) V log 2 (1/mR) J v ; 

This renormalization-group improved expansion resums the standard (renormalized) perturbative ex- 
pansion taken at fixed renormalization scale ~ 1/R in the 0(n) model. Specializing to n = 6 and 
re-expanding in terms of the string tension y/\ = Aug 3> 1 leads to 

c(mfl)=5 -™+o(^), (1.34) 

indicating that the stringy loop expansion will receive enhancement by powers of log R ~ log log S at 
higher loops. This can only make sense at large spin if log R <C vA, which eventually is just a rephrasing 
of the condition mR <C 1, since m ~ exp (— vX/4) . 

Our above discussion illustrates how rich and complicated is the structure of the subleading large 
spin corrections to the twist-two anomalous dimension and how subtle the interpolation between the 
gauge and string predictions becomes. 

Our subsequent consideration is organized as follows. In the next section we will formulate the 
long-range Baxter equation and construct its solution at large spin for the twist-two operator, up 
to the relevant accuracy. As a follow up, we shall derive the all-loop expression for the bulk (i.e., 
regular) and alternating contribution to the energy, both encapsulated in the solution to a system of 
integral equations. Then, in Section [3l we advocate and postulate the Liischer formula for excitations 
propagating on the GKP background. We discuss there the rotation to the mirror kinematics for all 
flavors of excitations and provide all-order representation for their associated Y functions. Later on, 
we perform the analysis up to five loops and demonstrate the equivalence of our formalism to the one 
available in the literature. Several appendices summarize technical details. 



11 This result was derived from the leading [52] and subleading |55] large n expression by assuming that the two-loop 
coefficient, obtained after re-expansion in the renormalized coupling e 2 of the O(n) model, is a polynomial in n of degree 2 
vanishing for n = 2 (i.e., for free theory). The latter re-expansion is found by means of 1/ log (f /mR) ~ /3 e 2 ~ (n- 2)e 2 . 
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2 ABA solution 



In this section, we shall derive the ABA expression for the twist-two anomalous dimension in the 
large-spin expansion up to the order 0(1/ S 2 ). We will make extensive use of the long-range Baxter 
equation [351 SS] j whose content is equivalent to the Asymptotic Bethe Ansatz equations of Ref . [5j , and 
apply the strategy developed in Ref. [47J for construction of corresponding large-spin solutions. After 
exposing the main ingredients of the Baxter equation in the next section, we engineer two non-polynomial 
solutions, — the main building blocks of our framework, — out of which the sought polynomial solution 
(|1.14p can be obtained. Then we address the next step of our formalism and derive a system of equations, 
similar to the ones worked out for the cusp anomalous dimension [T7] , that allow us to extract appropriate 
scaling dimensions to all orders in 't Hooft coupling. 

Before we move on to solving the Baxter equation, let us comment on the parameter we shall use 
to carry the large spin analysis. As we will establish later in this section, the large-spin expansion is 
naturally carried out in terms of the parameter rj that is determined by the conserved charge c\2 = f] 2 
entering the transfer matrix t(u), 

r, 2 = C 2 + S. (2.1) 
The latter is essentially defined by the renormalized quadratic Casimir 

C 2 = (S+i 7 + l)(5 + i 7 ) (2-2) 

of the underlying sl(2) algebra. This is the case up to the additive contribution 5 which scales log- 
arithmically with the spin 5 = x log fj + y + 0(77 ), with log 77 = log rj — ip(l) = log 7/ + 7e, with 7e 
being the Euler-Mascheroni constant, and where x and y are functions of the coupling only. Since our 
analysis relies on the ABA equation, the anomalous dimension 7 = A — Aq in (|2.2p stands for its ABA 
prediction, 7 = 7 ABA . The ABA equation, or equivalently the Baxter equation, also provides us with 
the coefficient 5, related to some higher spin-chain charges. 

At the level of the ABA equation, the relation (|2.ip between r\ and S can be unambiguously computed 
for any S and/or expanded at large spin to any accuracy. For many reasons, however, we might prefer 
to restrict ourself to this part of the relation (|2.ip that is not contaminated by wrapping effects. At 
weak coupling, this happens if the relation r\ = rj(S) = S + . . . is worked out at large spin up to (and 
excluding) l/S 2 corrections. This stands for the fact that the main ingredients entering (|2.ip . namely 
7 and 5, are expected to be uniquely determined at large spin up to corrections ~ I/77 2 which is the 
typical scale for wrapping effects [381 EH] ■ I n our discussion of contributions ~ l/S 2 to the anomalous 
dimension, the suggested accuracy for determination of rj = rj(S) is good enough^ We stress that 
these remarks apply to any order at weak coupling (only). At strong coupling, wrapping effects enter 
earlier [33] with the scale ~ 1/logS, truncating its expansion at large S to the first asymptotic term, 
tj(S) = S + . . ., i.e., right before this scale enters its right-hand side. This limitation will not affect 
however our discussion of the strong coupling regime. 

2.1 Long-range Baxter equation 

The Baxter equation that we will be solving here was given above in Eq. fjl . 13[) . Below, we will discuss 
in turn two important ingredients determining it, the dressing factors and the transfer matrix. The 
notations that we will be using presently will be however different from the ones adopted in Refs. 

nans] 

12 This statement follows from the logarithmic scaling of the anomalous dimension, 7 ~ log 77 = log (S + i^S) ~ logS + 
SS/S + . . ., such that contributions ~ l/S to the parameter 8S cannot affect 7 at the working accuracy. 
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2.1.1 Dressing factors 



To start with, let us notice that the ratio of the dressing factors entering the left-hand side of the Baxter 
equation (11,13|) could serve as a generating function for the spin-chain conserved charges, 



A-(u) 



1/2 



1 — g 2 jxx- 



g 2 /xx~J 



exp 



which are defined as [16J, 



i=i 



1 

27 



n>l V 7 



log ■ 



(2.3) 



(2.4) 



These are defined in terms of the renormalized rapidity x = (u + -^/m 2 — (2g) 2 )/2 which is the large 
solution, i.e., x(u) ~ u as u — > oo, of the Zhukowski map u = x + <7 2 /x. Here and elsewhere in the 
paper we will be using the notations x^ = x (u ± |) . The contour of integration in (|2.4p encloses the 
cut along u 2 < (2g) 2 present in the x and goes in the anticlockwise direction. The lowest charge, with 
1, determines the ABA anomalous dimension or equivalently the energy, 



E 



ABA 



L + 2g 2 q 2 . 



(2.5) 



Recall that the energy is determined by the scaling dimension A minus the spin S. The latter stands 
for the number of magnons propagating on the spin chain of length L, where according to the spin- 
chain/operator dictionary the integer L reflects the twist of the operator, or the number of scalar fields 
Z building it up. 

On the other hand, the individual dressing factors contain more information than their ratio. They 
can be written as 



logA+(«) =Llogx + 2J2 n fn ( — 



n=l 



OO 

logA_(n) = Llogx + 2^7„ ( |-J 

n=l ^ 7 



(2.6) 



where the expansion coefficients 7 n , j n , are complex conjugate of each other. The latter quantities are 
dynamical: they depend on the state under study, that is on the solution Q(u) we aim at constructing, 
and thus on the coupling constant. Their precise form shall be given shortly. 

In view of the representation (|2.6p . the dressing factors A±(u) do not define analytic functions of u. 
At leading order in weak coupling, the expression (|2.6p reduces to 



A + (u) = A-(u) 



+ 0(g 2 ) 



(2.7) 



and, after plugging it back to (|l,13p . one recovers the Baxter equation for the short-range sl(2) Heisen- 
berg spin chain of conformal spin s = 1/2. However beyond certain order in the weak coupling expansion, 
the dressing factors develop poles at u ~ which resum into cuts with branch points at u = ±2g at 
finite value of the coupling. For this reason one needs to be careful about the meaning of the expressions 
A±(u ± |) entering ([1.130 : namely, the long-range Baxter equation is understood to be evaluated on 
the sheet where (|2.7|) hold true at large rapidity, u ~ oo. 

An equivalent representation for the dressing factors (12.60 can be found after introducing the gener- 
ating function ^(t) as an infinite series 



7 (t) = ^(2n) 7n J n (t) , 



(2.8) 



n>l 
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in Bessel functions J n (t), with well-known asymptotic behavior J n {t) ~ t n for t ~ 0. Similar equation 
holds for 7(i) via the obvious replacement j n — > j n above. We note that since for t real the Bessel 
functions are real, then j(t) = ("f(t))*. In terms of these functions we find 



f°° dt 

log A + (u) = L\ogx + -7(±2 fl t) e ±iut , 

r It (2 - 9) 

logA_(u) = Llogx+ / — j{^2gt)e ±lut , 



where the upper (lower) sign refers to u in the upper (lower) half plane, respectively. This immediately 
follows from the identity 

■in \ n r°° lit 

- ) =n -J n (±2gt)e ±iut , n=l,2,..., (2.10) 
x J Jo t 

where the db corresponds to 3?m(ti) ^ 0. 
Now, we can write the coefficients j n as 

j r (j^n^ / n \ n . if dx ( x \ n 



(2-11) 



Here the contour of integration goes anti-clockwise and surrounds the interval v? < (2g) 2 in all cases. 
The coefficients d n ,d n , are complex conjugate of one another and account for the modification stemming 
from the BES/BHL dressing phase [17\ I56j . It is important to realize that dm-\ is purely imaginary 
while d2n is real as it follows from their definitions. This implies that dressing phase only affects the 
imaginary part of the coefficients 72«-i and the real part of the coefficients 72 n - Therefore, in the 
representation for d's, that will appear to be convenient in our subsequent discussion, 

= r ^m hMfjhm , rf2 „ = r ? AM t(*i(a , (2 . 12) 

Jo t e —l Jo t e —l 

only the first term in Eq. (|2.11|) involving the Baxter polynomial will contribute to the right-hand sides of 
Eqs. (|2.12p . This property becomes obvious since the integrand of Eqs. (|2,12p involves the combinations 
7_l-(i) = [y(t) ± 7 (— 1)]/2 (and similarly for 7±(t)) possessing specific parity under the reflection t — > —t, 
which project out the (^-dependence on the right-hand side of Eqs. (|2.1ip . The representation of the 
dressing phase given in (|2.12p is just one among many. The reason to select it here is that it becomes 
convenient in large-spin analysis when one aims to cast the Baxter equation into the form of an infinite 
system of equations for the defining coefficients 7 n ,7n- 

We note that for a polynomial solution (|1.14p . the contour integrals in (|2.1ip are just a fancy way 
of rewriting of the finite series representation for -f n 
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n r-f \ ix 



d n , (2.13) 



and similarly for j n by complex conjugation. Making use of the reality properties of d n 's for even/odd 
values of its index, one finds that we have the following relations between the conserved charges (|2.4p 
and 7 's 

q2n = (-ir+y- 2n ( 7 2n-l + 7 2n-l) , </ 2 n+l = (" lf +1 9'^ (j2n ~ l2n) ■ (2.14) 

The relations (I2.14j) give us a clue on the typical scaling of the coefficients 7n,7n at weak coupling. 
Namely, assuming the charges q n are non-vanishing at one loop, i.e. q n = 0(g°), we get that 7 n ~ 
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In ~ g n ■ This could have been inferred directly from the representation (|2.1ip as well. As a particular 
application of (|2.14|) . we can find the ABA energy as 

E ABA = L + 2g (71+71) • (2-15) 

We see that the same quantities 7nj7n> parameterizing the dressing factors of the Baxter equation, 
immediately determine the energy we are interested at constructing. A fruitful strategy, which is easily 
implementable at large spin, is then the following: first, obtain a representation of the (large-spin) 
solution Q(u) to the Baxter equation with unspecified 7 n ,7n, and then derive a system of equations 
for these coefficients by imposing the relations (|2.11|) . The last step shall be to compute the energy 
via ipiS) . 

2.1.2 Transfer matrix t(u) 

The second ingredient in the Baxter equation is the function t(u). It can be decomposed into two 
contributions 

t(u) = p(u) + s(u) , (2.16) 

where p is a polynomial of degree L in u and s(u) is a non-polynomial function of u, whose presence is 
needed to match the analytical structure of the left-hand side of f|l. 13j) . 
The polynomial component, p(u), can be parameterized as 

p{u) = 2u L - qiu^ 1 - q 2 u L - 2 - . . . - q L , (2.17) 

where the coefficients q n generalize the auxiliary conserved charges of the short-range Heisenberg spin 
chain alluded to above, to which they reduce when g = 0. The singular part, s(u), admits another type 
of expansion, which can be chosen as 

,(«) = 2 £ sJ^L\ n + 2 £ s n (JL\ n , (2.18) 

n>l n>l 

with s n , s n , being complex-conjugate coefficients, both independent on u. An alternative representation, 
that holds for —1/2 < 3 : m(n) < 1/2, is given by 

s(u)= ™ s ( 2 gt)e iut - t / 2 + -s(2gt) e"™'-'/ 2 , (2.19) 
Jo t Jo t 

where s(t) = ^2 n> i(2n)s n J n (t) , and similarly for s(t) with s n — > s n . 

For a polynomial solution, Q(u) ~ u s , the two first coefficients in p(u) can be written as 

qi = -2^( 7l - , q 2 = j(j - 1) + \L + 2g 2 L + 2g 2 ( l2 + j 2 ) + g 2 ( 7l - 7l ) 2 . (2.20) 

where j is the renormalized (total) sl(2) spin of the state, 

j = S + g( 11 +^ 1 ) + lL, (2.21) 

since 5(71 +71) is exactly half of the ABA anomalous dimension, see Eq. (|2.15p . The relations (|2.20p are 
a direct consequence of the asymptotics Q(u) ~ u s for relevant solution to the Baxter equation (11 . 13|) . 

At weak coupling, the polynomial part p(u) gives the leading contribution to t(u). The non- 
polynomial component s(u) appears suppressed, with the typical perturbative scaling for the expansion 
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coefficients s n being s n = 0(g n ). Their precise scaling actually depends on the state and on the spin- 
chain length under consideration. A general expression for these coefficients, in terms of the Baxter 
polynomial and dressing factors, reads 



1 



2(2ivr) 
1 

2(2ivr) 



dx 
x 

dx 

x 



A-(u) 



Q(u-i) 
Q{u-i) 



(2.22) 



Q(u + 



where the contour of integration goes counterclockwise and encompasses the cut along the interval 
u 2 < (2g) 2 in the u-plane as we discuss in Appendix [Aj The latter cut is hidden in the relation 
x = x(u). The contour of integration in (|2.22p should be chosen in such a manner that it encircles all 
singularities emanating from u = i| as the coupling increases. This includes the previously addressed 
cut, for finite values of the coupling. It is moreover not excluded that some of the singularities in A±(u), 
that originally resided on the second sheet of the Riemann surface \x\ 2 < g 2 , will migrate to the first 
sheet. It appears that this is not the case in our problem, as could be inferred from the analysis of the 
analytical properties of the solution to the BES equation offered in [24]. It is therefore safe to assume 
that the contour of integration in (|2.22[) surrounds closely the cut in u plane, that it equivalent to the 
circle of radius g for the variable x. 

Note finally that we can introduce a set of secondary, i.e., complementary to Eqs. (|2.14p . charges 
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Q2n = (-1) 



ir+y- 



2n 



2n-l 



72n-l, 



q 2n +l = (-l) n 5" (72n+72n) , 



which allow us to derive an alternative representation for s n 's 



(2.23) 



1 / dx 
2(2ivr) J ~~x 



x L exp 



-i$-iy j q n+ ix n + q n+ i ( — J 

n>l n>l \ X / 



(2.24) 



where $ is the spin-chain total momentum, 



s T + 
x j 



n 



3=1 J 



(2.25) 



The expression (|2,24p is convenient for the weak coupling expansion. Indeed, the contour integral 
in (I2.24p is easily performed then, after using the known behavior of q n ~ q n ~ g° and expanding the 
integrand at weak coupling (and fixed x). 



2.1.3 Comments on twist-two operators 

Since we are considering a length-two chain only in the current work, let us indicate the simplification 
that apply to this particular case. 

The polynomial twist-two solution at a given (integer) spin S is unique and is described by a (real) 
symmetric distribution of roots, such that Q{—u) = (— l) s Q(u), with the change of sign for odd spins 
reflecting the presence of the Bethe root at u = 0. This implies that the coefficients 7 n are real, 7 n = 7 n , 
and, in particular, the twist-two anomalous dimension is merely determined by one of them, 

E ABA = 2 + 457i _ (2 2 ^ 

13 Note that the complementary charges are real, by their definition. They do not seem to have a simple spin-chain 
interpretation, however. 
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It also follows that the polynomial part of t(u) is completely fixed, with in particular qi = 0. Therefore 
we can write 

p(u) = 2u 2 - rj 2 , (2.27) 

where 

V 2 = q 2 = j(j - 1) + \ + V + 4 5 2 72 , j = S + 2 91l + 1 . (2.28) 

This is the definition of the quantity r\ = S + . . . that we shall employ as an expansion parameter at 
large spin, in place of the spin S itself. It is of type (|2.1j) suggested before, since C2 = j(j — 1), with 
5 = t; + 4g 2 + 4g 2 72. The latter coefficient scales logarithmically at large spin, 72 ~ log??, as do all the 
quantities j n , including in particular 71 governing the energy (|2.26p . 

In the same vein, the coefficients s n that control the non-polynomial part of the transfer matrix are 
all real, again as a consequence of the parity of the twist-two Baxter polynomial. We stress furthermore 
that s(u) vanishes identically in the first two orders in perturbation series and starts from three- loop 
order for twist two operators. With regards to the large-spin limit, it is important to realize that the 
coefficients s n are of order ~ rp ~ 5° at weak coupling, up to logarithmically enhanced factors whose 
maximal power increases with the loop order. Finally, another distinguished feature of these coefficients 
that shows up at the leading ~ r/° order at large spin is that they acquire the (— l) s spin-alternating 
character. This is due to the presence of the overall phase exp (i?9) in their expression, see Eq. (|2.24p . 
which, as is well known, evaluates to exp = (— l) s for twist two operators. 

Keeping these remarks in mind, we are now equipped with all required machinery to attack the 
construction of the large-spin solution to the Baxter equation that will eventually determine all the 
coefficients 7„. 



2.2 Non-polynomial solutions and large-spin expansion 

The Baxter equation (|1 . 13[) is a second order finite-difference equation which admits two linearly inde- 
pendent solutions. When solving the equation order by order in the coupling constant for finite values 
of the spin S, one of the solutions is chosen as polynomial (j!.14j) while the other one is not. Only the 
first one corresponds to the eigenstate of the eigenvalue problem for the anomalous dimension mixing 
matrix. However, when considering the large-spin expansion, it appears fruitful to look for a pair of 
solutions Q±(u) both of which are non-polynomial. They enter in our analysis through the ansatz 

Q{u) = Q + (u) + Q_(u), (2.29) 

where the left-hand side is of course the sought polynomial solution. We shall also demand in the 
following that the relative normalization in (|2.29p is fixed by the extra condition Q-{u) = (Q + (u*))* . 

The rationale behind the decomposition (|2.29p is the observation [37j that the sought solution to 
the Baxter equation, Q(u), admits two different asymptotic expansions at large spin depending on the 
sign of the imaginary part of u. Namely, to arbitrary order in 't Hooft coupling the two functions Q± 
possess the behavior 

Q±(u)^^ 2iu , (2.30) 

for rj ~ S 3> 1 and fixed u ~ 0{rf). As a consequence, the function Q + (u) dominates in (|2.29p for u 
in the upper-half plane, and analogously for Q-(u) when u is in the lower-half u-plane0 Moreover, 
due to their specific scaling (I2.30P the functions Q± satisfy simplified equations at large spin. These 
are obtained by neglecting one of the two terms in the left-hand side of the Baxter equation (|1.13|) and 
read [47] 

A ± (u±£)Q±(it±i) =t{u)Q±{u). (2.31) 

14 Notice however that for real values of u both Q-(u) and Q+(m) are of the same order. 
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In fact, since Q±(u =F i)/Q±(u) ~ V 7 ? 2 while Q±{u ± i)/Q±{u) ~ rj 2 and t{u) ~ — rj 2 , not only do the 
simplified equations (|2.31|) apply at large spin, but they are also valid up to corrections of order 0(l/r] 4 ). 
For our purpose of constructing the large-spin solution up to order 0(l/r] 2 ) it is therefore sufficient to 
consider (|2.3ip . 

The construction of the large-spin solution along these lines was proposed in [37j at one loop, to 
leading order at large spin. It was further extended to higher loops in [57] and to all loops in |58| 
where it was shown to be equivalent to the BES equation [T7] , which is based on the Bethe-root density 
approach [59]. To two-loop order, but more accurately in the spin, the analysis of the large spin Baxter 
equation was refined in [60] and proved to reproduce correctly the large spin expansion of the two-loop, 
twist-two anomalous dimension of [18]. Our analysis below essentially consists in a generalization to 
higher order, either in the loop or spin parameter, of all previous considerations. It reproduces, in 
particular, the BES equation and its siblings [27] for the leading order solution and consistently reduces 
at two loops to the solution found in ]60] for the subleading ~ 1/rj 2 contributions. 

We put forward this agreement with established analysis because the equation (|2.29p on which it 
is based is to some extent ambiguous. Though, in practice, the functions Q±(u) are easily constructed 
in a iterative fashion, following [13 [60], the decomposition (|2.29|) and asymptotics (j2.30j) at large spin 
does not seem to fix them unequivocally. It is always possible, for instance, to find new solutions, given 
a pair of solutions Q±(u), as Q±(u) + aQ(u). We did not identify a convenient criterium that would 
pick out uniquely two such solutions among all possible pairs. Nevertheless, since we shall explicitly 
construct such a pair of solutions at large spin, it is always possible at the end to test if they sum up 
correctly to a polynomial solution, which is in principle unique. In this regard, we shall verify that our 
functions define a solution via (|2.29p that is holomorphic at the working accuracy in the spin. To ensure 
that it is the correct one would require testing that it possesses the expected asymptotics ~ u s at large 
u. This would turn out to be difficult because it would entail considering the regime u ~ rj ~ S, and 
ultimately u S> rj, which are not covered by our analysis. We shall therefore content ourselves with 
partial agreement alluded to above, which together with our ability to reproduce all known results up to 
four loops, make us confident about the validity of the above solutions. It would be interesting, however, 
and perhaps crucial for analysis at higher order at large spin, to address more seriously the problem of 
uniqueness that we emphasized before. 

Prior to presenting our results for the functions Q± , let us comment a bit further on the magnitude 
of corrections stemming from dropped terms of in the Baxter equation. The above estimate ~ 1/rj 2 
was based on the weak coupling expansion. However, due to the logarithmic enhancement displayed 
by the power suppressed contributions, the hierarchy observed between various corrections within the 
context of the perturbative expansion, which follows from the behavior (|2.30p . does not survive at 
strong coupling. For instance, corrections of the type 1/rj 2 and l/ry 4 both receive some log-enhanced 
contributions g 2n log™ r/ that change the analytical properties of the result when resummed to all orders 
in 't Hooft coupling. In fact, it is expected that all corrections of the form ~ l/r] 2k shall become of 
the same order when (7—7-00. Then, as one undergoes a transition from weak to strong coupling, the 
typical scale could change from I/77 to l/logry, the latter standing for the size of corrections found at 
strong coupling within the ABA framework [33]. For this reason, when referring to ~ 1/rj 2 corrections 
in the following discussion, we shall always mean contributions which are of this order to all loops, i.e., 
we shall use the 'weak coupling scale' for the ordering of the large-spin corrections. 

Now we are in a position to give the solutions to the truncated Baxter equations (|2.3ip . Namely, up 
to ~ 1/rj 2 corrections, we find 

, „ , , , , /, 1 \ 2iu 3 f°° dt e iut J (2gt) - iut - 1 

logQ + (u + ) = logc+-2m log 77-— T + " ' 



V J t e*-l 



00 dte iut j(2gt) - 2 gil t 2 f°° dt e iut s + (2gt) 1 r°° dt 
o T e*^l ^ J o J -rzi ^ J o J 
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where u + = u + i/2 and c+ is an arbitrary constant. Due to the last term in the right-hand side of this 
equation, the representation (|2.32|) holds for 9m(ti) > 0. Notice that, in distinction, the other integrals 
in (|2.32p only assume that 3 : m(u) > —1. Analogously, 

1 \ 2iu 3 [°° dt e~ iut J (2gt) +iut-l 



( 1 \ 2iu f 

logQ_(u~) = logc_ + 2iu ( log 77 - J ~ g-j- ~ 2 J 



t e'-l , 

(2 33) 

f°° dt e- iut 1 (2gt)-2g ll t 2 f°° dt e~ iut s + (2gt) 1 [°° dt _ iut 

where u~ = u — i/2 and with now 3m(tt) < and c_ = c* + . In both equations, we used the definition 
log 7/ = \ogr] + 7#, with V( z ) being the digamma function and 7e the Euler-Mascheroni constant, 
7e = — V'(l)- Also we introduced a parity-even combination s+(i) = [s(i) + s(— i)]/2 of s(t) defined in Eq. 
(|2.19p . We notice finally that the functions (|2.32|, [2733]) have the property that c_Q + (ii) = c + Q_(— u), 
wherefrom we find 

r- D (—,A 

(-1) 5 , (2.34) 



c + Q+(+u) 

after making use of the condition Q(—u) = (— l) s Q(u) obeyed by the twist-two polynomial Q(u) = 
Q+{u) + Q-(u). This identity will be useful later on in our discussion. 

It would not be difficult to verify that (I2.32h and (I2.33P are solutions to the equations (|2.31h at the 
given accuracy in inverse spin expansion. We shall see later that they also obey the scaling (|2.30p at 
weak coupling. Before that however, we would like to emphasize that even though the sum Q(u) = 
Q+(u) + Q-(u) should be regular over the entire complex u-plane the individual functions Q±(u) are 
not. Indeed, it turns out that that these functions develop branch point singularities lying inside the 
strip —2g < 9ie(ii) < 2g, see Fig. [2j This is not evident in the realm of perturbation theory where, at 
any given order of the weak coupling expansion, both Q±(u) appear as meromorphic functions, though 
with increasing order of the poles. The presence of the cuts, however, implies that the equation (|2,3ip 
for the functions Q±(u) apply only upon the additional requirement that the expressions Q±(u ± i) are 
understood as argument shifts performed outside of the strip alluded to above. This is only under this 
restriction that the functions given below can be observed to solve (|2.3ip . 

Let us look at the weak coupling limit. In this case, we can neglect the second lines in Eqs. (|2.32p 
and (|2.33p which enter starting from two loops only. This is because the functions 7 and s are small 
at weak coupling. Then, after simplifying the first lines by keeping only the first term in the Taylor 
expansion for J$(2gt) = 1 + 0(g 2 ) and applying the integral 



if „iut _„•»./ _ 1 

=■ - t = 1 °S r ( 1 " iu ) + ^(!) » ( 2 - 35 ) 

t e r — 1 



one immediately finds that 



Q±(u ± ) = i C ^. , 2 exp 

r(i =f my 



T2iu \ lo SV-Y^)±^ + o(l/v 



2iu 3 

(2.36) 



We explicitely verify the behavior f|2.30j) which holds as long as?i< r/J^j After expansion at large values 
of the spin S, by using the one-loop relation rj 2 = S(S + 1) + |, this expression is find to be in complete 
agreement with the earlier results of Refs. [471 160] for the leading and 0(1/ rj 2 ) corrections, respectively. 
As one can see explicitly from this result both Q±(u) appear regular over the entire complex u-plane, and 
so is their sum Q(u). This property persists at two loop order, but breaks down starting from three loop 
and is reflected in the appearance of poles at u = ±(i/2— in), with n = 0, 1, 2, . . ., in Q±(u), respectively. 
The regularity of the solution is restored however in the linear combination Q(u) = Q+(u) + Q-(u) 
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Figure 2: Singularity structure of solutions Q±(m) of the truncated Baxter equation. The Q±(u) are 
holomorphic functions in the upper /lower half of the complex plane at leading order of the large spin 
expansion and possess an infinite number of (blue) cuts for 3m(u) 5s 0, respectively. However, they 
develop (red) cuts in physical region once subleading, i.e., 1/n 2 , effects are taken into account. 

within the given ~ 1/f] 2 accuracy. This important property can be verified, for any value of the 't Hooft 
coupling, as follows. 

Already at leading order of large-spin expansion, the functions Q + /_ (u) develop an infinite number 
cuts in the lower/upper half-plane, respectively, with branch points at — i/2 — in ± 2g and i/2 + in ± 2g, 
with n = 0, 1, 2, . . . , as demonstrated by the blue cuts in Fig. [2j These singularities do not compensate 
each another but instead add up such that Q(u) appears to have an infinite number of cuts in the strip 
— 2g < 3f?e(u) < 2g rather than being regular. However due to the fact that Q±(u) ~ n T2lu , it is easy to 
convince oneself that the singularities at ^sm(u) = ±i/2 are actually of order 1/n 2 in accord with the 
conclusion arrived from the weak coupling analysis sketched in the previous paragraph. Namely, if we 
consider the vicinity of 9m(u) ~ i/2 we find immediately that Q(u) is dominated by Q+(u) ~ n which 
is obviously regular to leading order of the large spin expansion while Q-(u) ~ l/?7, which is singular 
there, is actually suppressed by 1/n 2 compared to the former contribution. So in order to confirm the 
regularity of Q(u) for 3 : m(ti) ~ i/2 we should observe that Q+(u) has a cut of order 0(l/n 2 ) (see the 
red cut in Fig. [2]) which compensates the singularities introduced by Q-(u). This is indeed the case, 
thanks to the presence of the last contribution in Eq. (12. 32ft stemming from the non-polynomial effects 
in the transfer matrix. Let us see how the compensation mechanism works in more details. 

To order 0(l/n 2 ), which is the accuracy of our considerations, it suffices to consider the 3 : m(n) ~ i/2 
domain more carefully to define Q(u) in the upper-half plane. Analogous analysis will be applicable for 
the lower semi-plane with obvious modifications. We start by rewritting Q{u) making use of the Baxter 

15 Higher-loop contributions would not affect this behavior at any given order, though they will come with a logarithmic 
enhancement. 
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equation obeyed by Q-(u) as follows 



Q(u) = Q+(u) + Q_(«) = Q+(u) 1 + 



-(«) 



Q+(«) 



Q+(u) 1 



A_(u- i)Q. 



7] 2 Q+(u) 



(2.37) 



where the last equality holds only up to o(l/r/ 2 ) corrections. To the same accuracy, we can replace 
Q±(u) by Q(u) in the second term 



\ogQ(u + |) ~ logQ + (u 



A_(u)Q_(^-§) 



A_(n)Q(»-§) 
2j ^(n + l) ' 



(2.38) 



which is valid for —1/2 < Qin(it) < 1/2. Then from the condition of logQ(u + |) to be regular around 
u = 0we should require that (see Appendix [A] for clarification of the origin of this equation) 



1 

2in 



dx 
x 



logQ(u + i) = 0. 



(2.39) 



for any n = 1, 2, . . . . We can check the fulfillment of this condition by plugging (|2,38|) into it. Then we 
immediately find that 



rj f dx 
2iir T x 



log Q+(n+i) 



i r dx 

2iir J x 

2s-n i 



A_(u)Q(u-%) 



(2.40) 



where in the last equality we used Eq. (|2.22p defining the expansion coefficient s n together with the 
reality conditions^ = s n valid for the twist two sector. Now substituting the explicit solution (I2.32|) in 
the left-hand side of Eq. (|2,40p . we immediately find that the only non- vanishing contribution comes 
from the last term in Eq. (|2.32p 



- \ r t etut <- 2 ^ = --2 e s « (- v c q+( u + 

T Jo t ^ t>\ ^ %x ' 



(2.41) 



It leads to the identity s n = s n in Eq. (|2.4U|) , after evaluation of the contour integral in the left-hand side 
of (|2.40p by taking the residues. This completes the verification of regularity of the complete solution 
Q(u) up to 0{l/rf) accuracy of our solution for 3 : m(ii) ~ 1/2. A similar reasoning applies in the 
lower-half plane, for 9m(ti) ~ —1/2. 

It is clear however that we cannot compensate singularities that are further away from the real axis 
in the complex plane (see Fig. [2]) if we do not go beyond the l/r] 2 approximation of our framework. 
For Qm(ii) ~ ±(i + n), the singularities produce effects of order 0(l/rj 2+2n ) such that for n > 1 they 
fall outside of the accuracy of our consideration. In conclusion our solution Q(u) is indeed regular 
everywhere in the complex plane up to o(l/rj 2 ) corrections. The proof of polynomiality however is 
beyond reach due to the order-of-limits problem, alluded to before, as this would require sending u to 
infinity prior to performing the large r/-expansion. Our construction was done in the opposite order and 
holds for u fixed and for 77 S> 1. 



2.3 Bulk and ABA energy 

The final step towards the construction of the ABA energy S ABA is to derive a closed system of 
equations for the coefficients 7 n introduced in Eq. (|2.6p . This is easily done with the help of the 
representation (|2.32p for the Baxter solution. In the analysis that follows, it will be instructive to 



21 



decompose all quantities into the so-called regular, or bulk, and alternating components as we advocated 
in the introduction. This separation is echoed by the expansion for 7n 's 

7n = 7 r k + 7f , (2-42) 

with each individual contribution obeying an independent set of equations. In the above decomposition, 
the first and second terms admit the following scaling with large r\ 

= < log + K + 0{l/r, 2 ) , 7 f = 0{(-l) s /r, 2 ) , (2.43) 

respectively. In both cases, we will find it possible to cast the system of equations for 7n in the form 

. , r dt r to ^ 7+(2gt)-(-l)" 7 i(2gt) 

In + J jJn( 2 9t) ~TZl = n ' ( 2 - 44 ) 

where n n is a given set of coupling- and spin-dependent sources which possess asymptotics behavior 
with rj analogous to one for 7n 's as will be shown in the subsequent sections. As should be obvious, Eq. 
(|2,44p generalizes the system of equations for the solution 7^ determining the cusp anomalous dimension. 
The latter controls the leading large-spin asymptotics of the bulk coefficients, see Eq. (|2.43|) . and solves 

(gam for 

n n = 2g8 n>1 . (2.45) 

Here (and below) b n ,m is the Kronecker symbol with well-known properties, 5 n ,m = {1, n = m;0,n / m}. 
On the other hand, the subleading function <5 7 ^ solves Eq. (|2.44p with the following inhomogeneous term 

A0 [°° dt J n (2gt)J (2gt) - gtSn,! 

K = 2j o J ^ • (2-46) 

It turns out that the solution 7 ^ for the leading logarithmic asymptotics plays a distinguished role in 
the construction of energy eigenvalues for more general fi^'s. Indeed, it allows one to project out the 
n = 1 component 7 j from the general-n solution 7n to Eq. ([2.45P with a generic set of source coefficients 
K n via the identity [61J 

4<7 7i = -^(-l)"2n 7 ^ n . (2.47) 

n>l 

The demonstration of this identity can be obtained from |61j . So from the knowledge of the coefficients 
K n and 7 ^, where the latter is a function of the coupling only, one can determine 7 j which controls 
the energy E ABA . Before we move on to the construction of explicit results for ^ ulk / alt obtained by 
means of the representation (12.47P , we will explain in the following few paragraphs how one achieves the 
separation between the bulk and alternating components and derives the system of equations (|2.44p . 
Here are the main steps. 

To begin with, let us decompose the solution Q+(u) in the upper half-plane as 

Q + (u) = Q b + ulk (u)Qf(u), (2.48) 

where the two factors have different scaling with 77, namely (u) ~ r]~ ~ 2iu and Qf{u) = 1 + 

0((— l) s /r] 2 ), respectively. Substituting the decomposition (I2.42P to (I2.32P and attributing non-po- 
lynomial parts of the transfer matrix to Q^^u), one deduces 

1 \ 2iu 3 f 00 dt e iut J (2gt) - iut - 1 



log Q b + uU > + |) = - 2iu (log fj - JLy) + ^ - 2 jf 



+ 



dt e iut 7 bulk (2 g t) - 2g^ ulk t 
T e*-l 



' (2.49) 
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and 

dte iut j ! " lt (2gt)-2gjf t t 2 [°° dt e iut s + (2gt) 2^ /g 



Jo t e — 1 V Jo 



V 2 Jo * e* — 1 rj 2 4^ V ix 



(2.50) 



where in the last term we made use of the equality 



., ^^^s(-2gt) = -^s n (^-)\ (2.51) 



to get an alternative representation for log Q^(u + |) which turns out more convenient for our purposes. 
At this point the following remark is in order on the analytical properties of the above two components. 
While the alternating contribution (^^(u) inherits the cut in the upper half plane with the branch 
points at u = i/2 ± 2g from the one of the full Q+(u), the Q+ ulk (ti) does not and is thus holomorphic in 
the upper half-plane. However, both Q^ ulk (n) and Q^{u) have infinitely many cuts in lower-half plane, 
with branch points at u = —in/2 ± 2g, n = 1, 2, . . .. 

Note that in order to evaluate 7 alt from Eq. (j2.11|) . we need to account for the presence of Q-{u) in 
Q(u) to stay within the accuracy of our approximation. To this end we first rewrite 

log Q(u) = log Q+(u) + log (1 + Y h (u)) , Y h (u) = ^4 ~ V" > (2-52) 

and observe that in the upper-half plane, 9m(ti) > 0, and particularly in the vicinity of u = i/2, 

Y h (u + §)~(-l) s A? 2 - (2.53) 

So Yh(u + |) is a small correction to log Q(u + §) in the upper-half plane that contributes at this order 
in the large-77 expansion to the alternating part only. This counting allows one to write 

logQ(u + §) = log Q + {u+r) + Y h {u + §) + ... , (2.54) 

where the ellipsis stand for 0(l/n 4 ) corrections. The second term Y^(u + |) contributes to the final 
expression for 7^* after plugging (j2.54p into (|2.1ip at power-suppressed level. So in summary, the two 
formulas we use for the extraction of ^ ulk / alt can De written as 



bulk _ 








In 


2m) 




^ix 


y»it - 








In 


~2i^) 


ff ( 


v, ix 



)bulk (u+ , ) + d bul k) 

logQf(n + §)-^yf (|) Y h (u+l) + dt , 



(2.55) 



where the separation of the dressing phase contribution into the bulk and alternating terms merely 
follows from the one for the coefficients 7 n in Eq. (|2.42p . 

The final step then boils down to evaluation of Eq. (|2.55[) with the help of (|2.49|) and (|2.5Up . This 
together with the expression for the coefficients d n defining the dressing phase and the identity 

1 I dx ( g\ n iut 1 I dx ( x \ n iut 



gives the system of equations ()2.44p with prescribed K n 
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2.3.1 Bulk energy 

Now we turn to finding the bulk energy making use of the representation (|2.47p . The source term 
defining it is rather simple and reads 

^ lk = (log fj - ^) < + <J< + C b + o{lh 2 ) , (2.57) 

where were previously given in Eqs. (|2.45p and (|2.46p . respectively, while the last contribution 

is 

C = z-8 nt i + ^-^S n 3 . (2.58) 



One could obviously absorb the term ~ <5 ni i present in K^ ub into the coefficient in front of in (|2.57p 
since the latter is localized at n = 1 as well, and then deal separately with the component ~ 5 n ^ in K^ h , 
which to an extent defines the only new problem. The decomposition (|2.57[ I2.58[) is nevertheless con- 
venient since it systematically separates contributions in the 1/rj 2 expansion that appear with different 
transcendentalities at a given order in weak coupling. 

If we are interested in the (bulk) energy only, _E bulk = 2 + 4g7 bulk , then only the coefficient 7 bulk 
matters, and, by means of Eq. (j2.47|) . we know that the former admits the representation 

^ - ^ (no - w) j MhM! _ V 7 . + p 3 (2 , 9) 

up to o(l/ rj 2 ) corrections. Therefore, from the knowledge of the solution to the BES equation, that is, 
from the coefficients 7^ or equivalently their generating function 7 (i) = S n >i(2n)7 J n (t), one obtains 
the expression for the bulk energy up to order 0(1/t] 2 ). As we can see it takes the form 

£ bulk = 2,4 (log fj - JL) + B + C/rj 1 + . . . , (2.60) 

where A = 2<?7 is the cusp anomalous dimension and B the twist-two virtual scaling function [27j . 

B = 2 + tgSjj =2 + 2 T MfkMMrt - ^ _ (2 . 61) 

Jo t e — 1 

where we cast it in the representation for the function B equivalent to the one given in [27]. The weak 
coupling expansion of the new coefficient C is given up to five-loops by 

C = -4 9 Sf + 4, 3 7l = "% 4 + ~ ^ + + 64 C l) 9 W + • • • • (2.62) 

Re-expanding the bulk energy in terms of the spin S we induce odd powers in the inverse spin in 
addition to the general structure in terms of the jy-variable, 

^buik = 2 AlogS + B+ C n/ sn + ■■■ ■ ( 2 - 63 ) 

n=l,2 

Here C n is a polynomial of degree n in log S = log S + 7e with coupling constant-depending coefficients. 
For the reason that the degree of the polynomials C n remains bounded at any value of the coupling, the 
bulk energy is quite gentle from the point of view of the interpolation between the regimes of weak and 
strong coupling. The explicit expressions for the polynomials C n are not quite illuminating and will not 
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be quoted here. They immediately follow from (|2.60p and the rj = rj(S) relation for twist two case, see 
Eq. (|2.28|) . Up to five loops, all what is needed are the expressions for up to n = 4, 

„ „ 2vr 2 o 22vr 4 . /146vr 6 >2 \ 7 /7096vr 8 32vr 2 Cl , ,\ o 



45 

„ 2tt 4 5 74vr 6 7 /326tt 8 on . . \ 9 

7 4 = 4C55 6 - + 84 C7 ) / + . . . , 

and similarly for 5j®, 



(2.64) 



/327r 6 C 3 2087r 4 C 5 1407r 2 C 7 , \ „ 

+ + + —3-^ + 48 C 3 3 + 3528C9 J g 9 + ■ ■ ■ 



hi 



4ir\ 4 , ^43vr6 _ \ 6 ^892vr 8 8tt 2 C 3 2 



2(3 3 / ' ) -' 



J 7 | = —g* - ^g* + — - 12Cl )g 6 - ~ - 208C3C5 U 8 + • • 

6 y 45 y V 630 / V 14175 3 



V-lO^+^-^ + lMft),' 



180 315 V 141 75 
Here we have also included coefficients that will be useful later. 



(2.65) 



2.3.2 Alternating contribution 

Finally, let us address the alternating contribution to the ABA energy. The source terms for this one is 
more involved but can be derived in a straightforward manner from (|2.50p and (|2.55p to take the form 



K a\t _ 



dt J n (2gt)s+(2gt) 



T Jo t 



-1 



1 
2Tvr 



dx 
x 



^— $ — I — 1 Y h (u + 



where s + (t) is found from the infinite series 

*+(*) = ^2(2n)s 2n J 2n (t) 



with the expansion coefficients that read 



1 



S2n 



dx 
x 



n>l 



2n 



2)i 



Y h (u + 



(2.66) 



(2.67) 



(2.68) 



Amr] 2 

up to corrections smaller than 1/rj 2 at any order in weak coupling. 

The system of equations (12.44f) for 7^* appears closed once one realizes that to compute k^, to 
arbitrary order in the coupling, one merely needs to know Y^(u + i/2) at leading order at large spin. 
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The latter function can be considered as being known once the solutions 7^ to the BES equation together 
with its subleading components #7 have been constructed!^! Then, after plugging the coefficients K alt 
into the system of equations ()2.44p and solving the latter, one gets the expression for 7 alt . 

Instead of solving the system of equations (|2.44|) with the inhomogeneous terms (|2.66|) . we apply 
the formulae (|2.47p in order to extract the coefficient 7 alt which controls the energy. The expression so 
obtained is parameterised by the solution to the BES equation, i.e., 7 solving (|2.44p with (I2.45j) . and 
the coefficients (|2,66p . It simplifes after using the fact that s+(2gt) in (|2.66[) can be decomposed over 
the even Bessel functions (|2.67p and then the equations obeyed by the solution 7 . Proceeding along 
these lines we derive the following representation for S alt 



E 



alt 



1 

2ivr 



dx 
x 



£ 

n>l 



2(2n - 1)7 



2n-l 



2n-l 



+ E2(2n) 72 n ^ 

n>i ^ y 



2n 



(2.69) 



which holds again up to o(l/r) 2 ) corrections. In this formula, the contour of integration goes around the 
cut along u 2 < (2g) 2 present in x{u) and Y^(u + |). 

To evaluate the integral (|2.69p at weak coupling, we can use that 



Y h (u + 



exp 



^ + 2^ 72n _i 

n>l 



2n-l 



n>l 



In 



(2.70) 



up to o(l/ rj 2 ) corrections (at any order in weak coupling). This representation follows from the consider- 
ation analogous to the one which lead to the integrand in Eq. (|2.24p . We see that due to the presence of 
the factor exp = (— l) s we immediately find that Y^(u + §) ~ (— l) s /rj 2 and thus E &lt ~ (— l) s /rj 2 
confirming and explaining the naming scheme we adopted for this contribution. 

Now, recalling that 7„ ~ g n for g ~ 0, we observe that we can compute -E' alt at weak coupling by 
expanding those terms in (|2.7(jp that involve inverse powers of x, e.g., (g/x) n , and then perform the 
contour integral in (|2.69p by the method of residues. For instance, the leading order contribution at 
weak coupling arises from 
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+ O(g 4 



(2.71) 



which in light of the dependence jf ~ g implies the £' alt vanishes at one loop. Proceeding to higher loops 
is straightforward, though slightly more tedious, and after some algebra one could derive the following 
expression 

alt_ ("1) S 



E 



8<? 2 7i7i72 + 12 5 2 7 3 7i + 8 9 2 7l ( 7 4 + 7 2 2 )(73 + §7?) - 8<? 2 7 2 (74 + 7 2 2 ) + • • • 



(2.72) 



which suffices for evaluation up to five-loop order F^l 

Finding the weak coupling expansion is now straightforward by substituting the large-// expansion 
for the coefficients j n , i.e., j n = 7 logj? + ^7 + 0(l/r] 2 ), together with the explicit weak coupling 
results (I2.64P and (]2.65p . Eventually we find i^ alt to admit the following generic form 
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alt 



Y,9 2n Pn + o{lh 2 



n>2 



16 This remark becomes obvious after looking at the Eq. (|2.70|l where Y^(u + ^) is given in terms of j n 
construct Y^(u + ^) at order 0(1/?? 2 ) we only need to use that 7„ = 7* log 77 + + O^/rj 2 ^ 



(2.73) 



Manifestly, to 



17 Note that the validity of p. 721) up to five loops also relies on the fact that the coefficients 7JJ are actually more 

n + 2 

ill JJICU^C VJ1 fn ' v y , CjV^C^fL 

2. 



where the naive scaling applies, 7J 1 ~ g. The latter properties allowed us to remove some terms in (12.721) . as 
~ g 2 737273 for instance, that eventually enter at 0(g 12 ) instead of 0(g 10 ). 
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where P n is a polynomial of degree (re — 2) in the log fj variable. As we shall see in a later section, the 
structure (|2.73|) is actually correct to all loops, though this is not transparent from (|2.72|) . This pattern 
contrasts with the simplicity of the bulk energy, which, even after re-expanding in the spin S, displays 
a logarithmic enhancement in log S that is strictly bounded in power at any coupling!^ We want to 
emphasize again that the corrections we are neglecting in (|2,73p are smaller than 1/rj 2 to any order 
at weak coupling only. These corrections are indeed expected to exhibit the same kind of logarithmic 
enhancement as the one accompanying the above 0(l/rj 2 ) contributions, and their resummation at 
finite coupling may considerably complicate the picture. At strong coupling, the contributions controlled 
by (|2.69p could not capture the leading large-spin expression of £? alt anymore 1^1 This result is anticipated 
to be the case at strong coupling in the regime where perturbative string theory is applicable and 
estimates for subleading large-spin corrections suggest that they scale as ~ 1/logS, see [30 | l3T j [32] . 

After this general remarks we come back to our calculation and present the weak coupling expansion 
for E 3 ^ up to five loop order. The polynomials entering Eq. (I2.73D are given by 

P2 = 0, 

-P3 = ^— log 77, 

P 4 = -128Cs log 2 fj + — — log fj + 16vr 2 C 3 , 

15 

128vr 4 , , , 9>\, 2 f~ > 2 8992vr 6 \ /320vr 2 C 5 896tt 4 C 3 \ 

P 5 = log 3 fj + (1280C 5 + 128^ 2 C 3 ) log 2 fj + (J68C3 2 - J log fj - ( — ^ + — ^ J • 

(2.74) 

According to our preliminary findings, we encoded the fact that i? alt does not receive corrections till 
two loops by starting the summation in Eq. (12.730 with re = 2 even though Pi = 0. The reason for this 
is that in distinction to the one-loop contribution P\ which vanishes exactly at any spin, the two-loop 
alternating contribution is non- vanishing in general but turns out to be suppressed as 1/S 4 at large 
spin. This can be easily verified starting from the exact (i.e., finite spin), two-loop, twist-two anomalous 
dimension [18J@ We also checked along the same lines that our result (]2.74p is in agreement with 
the three- loop and four- loop- ABA anomalous dimension obtained in Refs. [19\ 135] . The five- loop ABA 
anomalous dimension, though available |42j . is too cumbersome and difficult to handle at large spin to 
perform the comparison. 

Before we close this subsection, let us give one more representation for the alternating energy ()2.69|) 
which will be useful for our subsequent analysis at strong coupling. It reads 

^ = l ~j dE ^ u + D y h(n + §) + o(l/r? 2 ) , (2.75) 

where E^(u) is the energy of the scalar excitation (hole) computed in [481144] . The contour of integration 
in (I2.75D encircles the cut running from u = —2g — i/2 to u = 2g — i/2 (as shown in Fig. [7]) present 

18 The structure (|2.69p on the other hand reflects the fact that _B alt is non-linear in the coefficients y n ~ log fj, see 
Eq. (US]). 

19 Strictly speaking this phenomenon is not restricted to the alternating component of the ABA energy, because as we 
increase the accuracy in 1/rj we will generate corrections quite similar to the ones controlled by the right-hand side of (|2.69[) 
but with Yh(u + raised to some power. These corrections have no reason to alternate with the spin anymore and thus 
should contribute to the regular (i.e., non-alternating) part of the energy. One may expect however that there exists a 
more general decomposition of the ABA energy into a 'simple' bulk contribution and a 'complicated' remaining piece, 
with the latter part being identical to the alternating component of the energy to the leading ~ l/?? 2 order but otherwise 
receiving regular contributions beyond it. A possible example for such a remaining piece can be obtained from the right- 
hand side of (I2.69[) with Y^(u + i/2) — > log (1 + Y\ 1 (u + i/2)). This would be suggestive of a TBA-like interpretation of the 
decomposition between the presumed all-order- in- 1/rj bulk and remaining parts, which is however far beyond the scope of 
our analysis. 

20 Methods for separating regular and alternating contributions from exact representations usually involving nested 
harmonic sums can be found in or adapted from |62| . 
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both in Yh(u) and Eh(u). The proof of the representation (|2.75p goes as follows. First, we notice that 
the term in square brackets in the right-hand side of Eq. (|2.69p can be written as 



2.r 



dx 



£ 

n>l 



IX 



2 1 1 



n>l 
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(2.76) 



and then we use the formula 

S h (n + i) = l + ^7 2 „_ 1 



■ \ 2n— 1 / • \ 2n 

" -£*(t) 

7 n>l v y 7 



(2.77) 



The latter can be easily derived from the representation given for E^iu) given in [48] . Plugging it back 
into Eq. §ZM\\ immediately yields (|3.122|) . 



3 Finite-size corrections 

In this section, we will construct and analyze the expression for the finite-size corrections to the energy 
of the GKP string. 

To understand and motivate our proposal it helps recalling some distinguished features of the spec- 
trum and dynamics of excitations on the latter background. To start with, let us address the simplest 
set of excitations propagating on the GKP string, the holes. These naturally arise within the framework 
of the asymptotic solution to the Baxter equation. It will suffice for our purposes to restrict the consid- 
eration to one loop order. What we are looking for is the solution to Eq. (|1.13p valid for the spectral 
parameter scaling as u ~ O(S ). The latter condition merely reflects the value of the argument of the 
Baxter function set to calculate the anomalous dimension, i.e., u = ±|. The method is heavily based 
on the observation that the two terms entering the left-hand side of the Baxter equation (11.13P possess 
different asymptotics as S — > oo. The reason behind it becomes very transparent when one notices that 
the right-hand side behaves as ~ S 2 due to the appropriate scaling of the large-spin twist-L transfer 
matrix 

L-2 



%) = -7? 2 (3-1) 



where rj = S + Introducing the ratio of the Baxter polynomials 4>(u) = Q(u + i)/Q(u), one realizes 

that in order to match the scaling behavior of the left-hand side it should be either large or small. 
Therefore, it is legitimate to neglect one of the terms. This procedure yields two first-order finite- 
different equations instead of one second order. Then, the one-loop large-spin solution immediately 
reads 

. Til- 1 2 r(— iu + i6j) 
Q+ u = c +V - 2lu lh - l r(i y . , L J) , (3.2) 

where c+ is an overall normalization constant. From the regularity of the general solution Q(u) = 
Q+(u) + Q-(u), where Q-{u*) = Q + (u)* and the real-hole assumption, i.e., 9m(4) = 0, we get the 
condition 

1 + Y h (5 k ) = 0, (3.3) 

for the ratio 

Y h (u) = • (3-4) 

Q+(u) 
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Equation f|3.3[) is the requirement for the absence of poles in Q(u) and yields the following quantization 
conditions when one substitutes the explicit solution (|3.2|) in it, 

+v c+ni+is^l^-is,)- - ' (3 - 5) 

where the condition on the ratio of the unknown constants emerges from the total momentum condition 

(H2HD 

The equations (|3.3l 13,51 13. 6[) can be interpreted as the Bethe-Yang equations for the factorized 
scattering on top of the GKP string. Indeed the Y function apparently factorizes as 

L-2 

Y h (u) = e w e^l[S{u,S j ), (3.7) 

i=i 

with the two-to-two S matrix 

q( „ ^_ r(m-^-)r(i-»u)r(i + ^-) 
b{u > dj) ~ r(^ - in) v{\ + m)r(i - ' 



and the propagating phase 



tu 



e iP,(u) = (3 _ 9) 

The interpretation of the quantity 2 log 77 as defining the physical length then suggests to intro- 
duce the reciprocal momentum p(u) = 2u as it appears conjugated to the former in the progagating 
phase (|3.9p . The form of the latter remains problematic however since it bears an extra phase in (|3.9p 
that does not seem to admit a clear interpretation. The S-matrix also is a bit surprising in form since 
it is not a function of the rapidity difference. Both phenomena can be traced back to the presence 
of the dressing factors in the Baxter equation, which are typical to spin-chain problems. In a relative 
set-up for the periodic Toda chain, these factors are absent and the quantization conditions have the 
standard Bethe form where the S-matrix is a function of the rapidity difference only and the propagating 
phase is the exponential of the momentum times a parameter that plays the role of an effective length, 
see [631 EH EH] and references therein. These two features of the propagation and scattering on top of 
the GKP string are possibly related to the fact that this background is not quite homogeneous [H], as 
observed at strong coupling in explicit perturbative computations on the string theory side |31| . Though 
puzzling, these do not prevent us from trying to apply standard recipes to our problem. Note finally 
the presence of the total spin-chain momentum in (|3.7p . It endows the current consideration with an 
extra lever arm allowing us to vary possible boundary conditions for the boson, and deviate from the 
conventional case of zero quasi-momentum for periodic boundary conditions. 

The structure (13. 7p is quite generic and holds to given accuracy in the large-spin expansion for any 
value of the 't Hooft coupling, the factorization being guaranteed by the approximation and/or the 
fact that the problem we solve is equivalent to an open spin-chain problem. This way one can define 
an all-loop momentum and all-loop S-matrix by considering the all-loop twist-L problem at large spin. 
The so-obtained S-matrix was analyzed recently at both weak and strong coupling in Ref. [66] and 
matched successfully against the finite-gap theory string predictions!^] Here we are mainly interested 



21 The result of [66] at strong coupling applies to the giant-hole kinematics, i.e., for two holes with momenta ~ j > 1. 
In the opposite regime where the momenta are small, the holes scatter as the asymptotics excitations in the 0(6) sigma 
model \M- 
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in the propagating part of the Y function, as this is precisely to this expression that the Y-function 
reduces when evaluated in the vacuum, i.e., in the absence of small holes, and which corresponds to the 
twist- two configuration. 

The all-loop expression of the Y function for a scalar (hole) in the vacuum can thus be parameterized 

as 

Y h (u) = (-l) s e iP ^ u \ (3.10) 
where we used the fact that exp (ii?) = (— l) s holds for the twist two sector (and integer spin S). Here 

P h (u) = 2p h (u)logfj + 2Sp h (u), (3.11) 

where Ph(u) = 2u + 0{g 2 ) is the all- loop momentum of the hole excitation constructed in [38] and dp^u) 
is a new quantity which should account for, and deform at higher-loop, the O(log r/) correction present 
in (|3.9p . Note that Ph(u) is entirely determined by the vacuum (i.e., twist-two) solution. It accounts for 
the effective propagation of the scalar in the background of covariant derivatives. A concise expression 
for it will be given later. 

To define the Liischer corrections we would like to continue Yh to its mirror kinematics and then 
write down the finite-size correction to the energy in the form 

El s =~J d P r™(u)Y™™(u) , (3.12) 

where p™ irror (u) is the momentum of the hole with rapidity u in the mirror kinematics and is the 
total number of scalar flavors, i.e., = 6 for the case at hand. As we shall find from our subsequent 
analysis, in analogy to a relativistic theory, the dispersion relation in the mirror kinematics of the 
scalar is identical to the real one. This is the double Wick rotation symmetry advocated in [33], with 
p™ irror (u) = Ph(u) and E™ lvvoI (u) = E^(u). This shall immediately lead to the following form of the 
mirror Y-function 

Y h mirror (u) = (-1) 5 e-^W , (3.13) 
where £h{u) is the mirror of — iPh(u) entering (|3.10p . The former quantity will be given as 

£ h (u) = 2E h (u) log fj + 25E h (u) , (3.14) 

where Eh(u) is just the energy of the hole with rapidity u in the real kinematics and SE^iu) is the 
mirror of — i6ph(u), see Eq. (|3.1ip . that we need to construct. The latter is an extra contribution that 
reflects this curious feature of the large spin background that we already mentioned above. It plays an 
important role in the explicit computation of ()3. 12[) as it indeed provides the integration measure in Eq. 
(]3.12p that makes the integral convergent to any order in the weak coupling expansion. 

Of course the above 6 scalars are not the only excitations on top of the GKP string and the Liischer 
formula would not be complete if these new degrees of freedom were not included into the analysis. Of 
particular importance are the lightest excitations, with their twist equal to one, as they turn out to 
provide the leading contribution, ~ 1/S 2 at weak coupling, to the vacuum energy. This class of extra 
excitations includes 4 + 4 fermions and 1 + 1 gauge fields. Introducing both of them into the game can 
be done by extending the previous considerations to higher-rank subsectors of the ABA system. The 
fermions can be unscrambled, for instance, from the st(2|l) subsector, while the gauge fields require 
a more thorough in-depth analysis [67] • The complete list of excitations was studied in |48} 149] . In 
Appendix El we provide for reader's convenience a brief outline explaining how these emerge from the 
ABA. This consideration serves as a preamble to the consideration that follows below where we explain 
how to construct the propagating phases for these excitations in the real kinematics. Then on a case- 
by-case basis we pass to the mirror kinematics and eventually be able to evaluate the complete Liischer 
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formula for leading finite-size corrections at large spin 22 ! 

£ FS = " E Y [ dpf" m (u)Yr^(u) + . . . , (3.15) 

where the summation runs over all species of elementary (i.e., twist-one) excitations, * =h, f, gf, with 
n* the number of particles *, and with the dots standing for higher Liischer corrrections as well as for 
corrections from composite (higher- twist) excitations. For all excitations, we shall find that the double 
Wick rotation symmetry applies, with p™ 1TTOT (u) = p+(u) and E" 11IIOV (u) = E*(u). 

3.1 Excitations in the real kinematics 

In this subsection, we shall generalize our previous discussion about the effective Bethe Ansatz equations 
to the case of arbitrary excitations. These equations can all be cast into the form 

l + n(«*) = 0, (3.16) 

where u± is the rapidity of the ^-excitation. Our goal is to provide explicit expressions for all Y-functions 

n(n). 

For a collection of K+ excitations (of a given flavor), with * ^ h, we have 

n(«) = a* n S, b r>, U4j) II u^) , (3.17) 

3=1 3=1 

where 

a* = S^ e (u,u). (3.18) 

It turns out that cr* = — 1 for the fermionic excitation and <t* = 1 otherwise. Here S^ e (u,v) is 
the bare S-matrix for scattering of the excitations * and *' carrying rapidities u and v, respectively. 
It is different from the effective S-matrix controlling the quantizations conditions because the latter 
also receives contributions from the induced scattering, originating from the background of covariant 
derivatives parameterized by the distribution of main roots u^j, as reviewed in the Appendix [Bj The 
latter distribution gets deformed in the presence of impurities that generate the induced scattering. The 
explicit form of the S±± is not really relevant for our present discussion however, because at the end all 
we want is to evaluate (|3.17p in the vacuum, that is, with the help of the twist-two distribution of roots 
obtained in the absence of impurities, 

n(«)| waiJm = <r*nsSr(«,uy) J (3.19) 

3=1 



where the distribution of roots u^j solves the twist-two ABA equations. These roots are encoded in the 
Baxter polynomial Q(u) = Ylj=i(. u ~ n 4,j) that solves the twist-two Baxter equation (|1.13|) with L = 2 
and K/± = S. The resulting y-functions (in the vacuum) are always parametrized as a pure phase 

Y+iu) = a+i-lf e iP *W , (3.20) 



where P+(u) = 2p+(u)logr] + 5p*(u) + 0(l/r] 2 ), with p*(u) being the all-loop momentum of the *- 
excitation constructed in [38] . 



22 Our notations below are slightly abusive, since in conventional treatments the degeneracy factor n* associated to each 
contribution would actually be part of the quantity Y' 1 mlrror (M). 
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For the reader's convenience, let us quote the expressions for scattering of excitations on themselves. 
The corresponding bare S-matrix is trivial for fermions 

S£ are = 1 , (3.21) 

meaning that there is no bare scattering between fermionic roots. For the gauge fields and their bound 
states we found that it takes the form 

n=l x ' 

However, none of these will be relevant below since in all expressions provided below we will not take 
into account the presence of extra excitations in the bath. They will only apply to the vacuum state. 

3.1.1 Scalar 

We start with the scalar (i.e., the hole) excitation that we have already extensively explored before. 
To recapitulate, it is associated with a dual solution of the equations for magnons. As such it is not 
surprising that 

Y h (u) = Y magnon (u) = Q n f~At { rt ■ ( 3 - 23 ) 

Within the strip —1 < 3?m(u) < 1, which contains the real kinematics, we can use the fact that 
Q(u ± i) = Q±(u ± i) to leading order at large spin and then simplifies with the help of (|2.3ip . It leads 
to 

(3-24) 

which we also encountered before. We want to evaluate (I3,24p in the vacuum, i.e., without any excita- 
tions. This simply means that we should plug in the twist- two solution to the Baxter equation into this 
expression. Using Eqs. (|2.32p and (I2.33P we find it to admit the form 

7 h («) = (-l) s e I W, (3.25) 

where 

f°° dt ^{ut)e^U2gt)-ut [°° dt sm (ut) e*/ 2 j(2gt) 

Ph(n) = 4nl ° gr? + 4 /o T e^I 2 J J —1 ' (3 ' 26) 

and we recall that log f\ = log r] — ip(l) = log r/ + te- Note that along the way we also used the condition 
c_/c + = (— l) s , derived previously in (|2.34p . 

The one-loop estimate of (|3.26p can be obtained making use of the observation that the function 7 is 
small at weak coupling, thus leaving only the two first terms in the right-hand side of (|3.26p relevant as 
g —7- 0. The latter two simplify after noticing that Jo(2gt) = l + 0(g 2 ) at leading order which eventually 
leads to 

after applying the identity 

n f°° dt sin (ut) e*/ 2 -ut T(± + iu) n ,. lX , n on , 

2 T e l = ^log . '- + mi)u. (3.28) 

Jo * e _1 r(| - iu) 

This confirms the result p.9p we reported in preamble to this section. 
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3.1.2 Fermion 



We shall now turn to the fermionic case. As we adverted before, these excitations can be found by 
considering the sl(2|l) subsector of the ABA system of equations [5]. Here we are mostly interested in 
the phase accumulated while a fermion propagates through the twist-two background, characterized by 
the distribution of roots u^j. The latter quantity is controlled, as demonstrated in Appendix [HI by 



Y t (x) = -(-l) 5 e iPf ^ =- jj 



x — xj 



, X — X ■ 
3=1 J 



(3.29) 



This expression illustrates one virtue specific to the fermion: its natural rapidity is the Zhukowsky 
variable x, and not u. To cover the whole x plane we need two copies of the u plane, to which we 
refer to as the small |x| < g and large \x\ > g fermion domains. In all cases u is related to x via the 
Zhukowski map u = x + g 2 /x, but its solution, x(u), is chosen to be the large one for the fermion with 
\x\ > g, such that x ~ u at large u, and respectively the small one for the fermion with \x\ < g, such 
that x ~ j 2 /m at large u. The fermion possesses a real value of the momentum p when x is real which 
implies that u 2 > (2g) 2 , independent on whether |x| is greater than g or not. In particular, a fermion 
at rest with p = or equivalently x = is small with a rapidity u ~ oo. 

In practice, it is more convenient to present expressions for the Y-function parameterized by the 
rapidity u. It follows that we need to provide two representations, one for the large fermion with rapidity 
u and another for the small fermion with the same rapidity. To distinguish between the two, we shall 
use the script 'If and 'sf, respectively. Of course these two expressions represent the same function of 
x, but in the two disjoint domains \x\ ^ g. By continuity they should match on the circle \x\ = g which 
is realized in it-space by the interval u 2 < (2g) 2 . The two expressions below are such that they both 
have square root branch points at u = ±2g and are glued together by an analytic continuation along 
their common cut along u 2 < (2g) 2 . 

Let us consider first the large fermion, i.e., carrying the rapidity \x\ > g. In this case it is convenient 
to rewrite the ratios involved in ()3.29j) in the form 



JPlf(u) 



M> s n 



u — u\ JL 1 — g 2 /xx j 



-.11 — u. 
J= l 3 ]=l 



Hi 



g 2 /xx] 



(-1)' 



: Q(u-h) 
Q(u + i) 



A_(-u)l 1//2 



(3.30) 



where the last equality is a straightforward rephrasing done with the help of the definition of the Baxter 
polynomial and dressing factors. We want to evaluate the above expression at large spin with the rapidity 
variable u belonging to the strip —1/2 < 9m(«) < 1/2, since the latter contains the real kinematics 
u 2 > (2g) 2 . In this strip, we can approximate Q by Q± depending on the sign of the imaginary part of 
its argument and get 



3 iP lf (u) 



("I)' 



■ Q-(u 



Q+(u + 



A_(«) 



1/2 



(3.31) 



Now it is an easy exercise to plug in the large-spin solutions (12.32H and (I2.33j) into Eq. (|3.31j) and after 
a short derivation conclude that 



„ , , A . _ . f°° dt sin (ut)J (2gt) - ut n [°° 

Pi f (u) =4«log7? + 4 / — v ; p f ; 2 / 

Jo t e — i Jo 

Here we also used the following integral representation 



dt sin (ut)j(2gt) 
1 e*-l 



dt 



sin (ut)'f-(2gt) . 

(3.32) 



A+(n) 



n 1/2 



exp 



n>l 



72n-l 
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2n-l 



exp 



r °° dt 

i I —sm(ut)~/_(2gt) 
lo t 



(3.33) 



33 



where x is the large solution u = x + g 2 /x. Here the equality after the second equality sign holds under 
the assumption that u is real and such u 2 > (2g) 2 , and thus the same restriction applies to Eq. (|3,32l) . 
The expression (|3.32[) is all we need for computing P\f(u), at any coupling, up to corrections explicitly 
inversely suppressed in spin. 

The small fermion domain can be deduced by analytically continuing (|3.32p through the central cut 
along u 2 < (2g) 2 . Indeed this corresponds to the continuation to the domain \x\ < g. It can be done 
making use of the equations obeyed by the function j(t). We do proceed differently, however, starting 
over from Eq. (I3.29P and then writing 



e iP a{ (u) 



III 



1 — x/x^ 



x/Xa 



A+(n) 



■nl/2 



A_ 



it 



exp 



n>l 



2n-ln 



(3.34) 



where x is now the small solution of the map u = x + g 2 /x, with x ~ g 2 ju at u ~ oo. Here we simply 
factorized T[jX~j/xJ = exp (it?) in (13.29p . which is equal to (— l) s for the twist two sector, and then 
applied several identities established earlier. For the real rapidity u, away from the cut u 2 > (2g) 2 , we 
can also rewrite ()3.34p as 

(3.35) 



f°° dt 

P s{ (u)= / -sm(ut)^(2gt). 
Jo 1 



For illustration purposes and preparing ourselves to the analysis that follows, let us construct the 
weak coupling expansion for Y\f(u). To this end, we neglect 7 in (|3.32p . approximate Jo(2gt) by 1, and 
find that the one-loop expression is simply given by 



Y]i (u) = -(-l) b rj 



S„4iu 



r(i 



r(i + iu) 2 



(3.36) 



where we used the result 



dt sin (ut) — ut 



-1 



T(l + iu) 
ilog-^ r{ + 2#u. 

r(i — iu) 



(3.37) 



This is the expected form for the fermion which has conformal spin 1. It comes with an extra minus 
sign as compared to the findings of |47j which merely is due to the statistics. For a small fermion, we 
simply find Y s f(u) = (— l) s , which is essentially equivalent to the limit u — > of (|3.36p . 

3.1.3 Gauge field and bound states 

Last but not least, the gauge field excitations and bound states thereof, are embedded into the ABA 
equations in the form of stacks of strings [67], The details of this embedding are recalled in Appendix [Bl 
By applying the fusion method to these complex solutions, we obtain the following representation for 
the function Y^{u) associated to a bound state of I gauge fields, 



Y t {u) = n 



-e] 



xjl 



g 2 /x [ 



>X 3 



Xj 1 



Q(u-T> 



A_(u+f)A_(« - f) 
A + (n+f)A + (n-f) 



1/2 



(3.38) 



where 



x[u 



u ± ^£ and m 



1 + 



At large spin and for u residing within the strip 



-f < 9m(u) 



< 



2 • 



we can evaluate (|3.38p with the 



help of Q±(u) determined in Eqs. (|2.32p and (|2.33p neglecting corrections 
here. Then Yi(u) if found to admit the expected generic form 



I/77 2 which are irrelevant 



Yt(u) = (- 



(3.39) 
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where 



dt sin (ut) e~ a / 2 J (2gt) - ut o f°° dt r ^ n _ a/2 



P t (u) = Aulogfj + 4 / — — — v 2/ — sin(ut)e 



o 



1 _ e ~* e* -1 

(3.40) 

It is easily verified that to leading order in large spin, Pe(u) = 2pg(u) log fj + O(log fj) with pi(u) being 
the momentum found in Ref. [48] . 

The one-loop expression is obtained from (|3.40p along the same lines as before. We get after a brief 
back-of-the-envelop calculation 

Y £ (u) = (-1) V* j^gg , (3.41) 

i (si + my 

where sg = 1 + 1/2. This is the anticipated result for excitations with conformal spin sil^l 
3.2 Mirror kinematics and double Wick rotation symmetry 

For real physical kinematics the energy and momentum are both real. For bosons this holds provided 
the rapidity u is real while for fermions it requires the reality of the corresponding Zhukowski variable 
x. The essence of the double Wick rotation is to look for a transformation u — > u' such that 

E^u') = ipf rrm (u) , p*(u') = iEf™(u) , (3.42) 

become simultaneously purely imaginary. We will explicitly construct this transformation in the case of 
scalar, gauge field and bound stated and find that it leads to 

Ef TT °\u) = E+(u) , pf Iim (u) = p*(u) . (3.43) 

These equalities are equivalent to the statement of the double Wick rotation symmetry of the dispersion 
relation proposed in |44| . 

Then, from the knowledge of the transformation u — > v! and the expression for the y-function 
derived earlier, we shall construct the amplitude for the propagation of an excitation in the vacuum in 
the mirror kinematics. We shall see that it takes the form 



lfe^W = (-i) 5 exp 



2E* (u) log fj - 25 E* (u) + O ( 1 /rf 



(3.44) 



where £*(tt) = —iP*(u'), with P*(u) introduced in (|3.20p . and where r) = S to leading order at large 
spin. Both the energy E*{u) and anomalous part 5E*(u) shall be found solely from the asymptotic 
solution to the Baxter equation. While the remaining 0(1/ rj 2 ) contributions in fj3.44f) are beyond of 
our current accuracy and are expected to receive finite-size corrections. They will not be considered 
here. The integral representation for the energy and its weak coupling expansion were already derived 
in [48t I44j . such that the new result here is for 5E+(u). This anomalous contribution parallels the one 
found in the real kinematics, as clear from the fact that 

SE+iy) = -i5p*(u') . (3.45) 



23 Note however that these twist-^ excitations are not quite identical to the excitations (i.e., holes) of thel spin-s^ Heisen- 
berg spin chain at the one- loop level. This is correct for the dispersion relation or the vacuum contribution to Yi(u) but 
not for the two-to-two S-matrix of such excitations. There is no contradiction here with the findings of [67] where only the 
case of length 3 was treated, implying a single excitation on top of the vacuum for us. The only exception to what have 
just been said is the gauge field, I = 1, where the mapping to the conformal spin chain with s — 3/2 is valid for any length 
(i.e., for any number of excitations for us). 

24 Note that the mirror transformation depends on the flavor of the excitation. Moreover it is different from the one for 
magnons, which should not be too suprising since dispersion relations have different analytical properties in both cases. 
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One immediate consequence of (|3.44j) at weak coupling is that 



y™ r (n) ~ 1/r/ 2 ~ l/S 2 , (3.46) 

at large spin, for all twist one excitation. In other words, the finite-size corrections to the energy shall 
be suppressed with the spin, as expected. This is entirely due to the asympotics 

E,(u) = I + 0(g 2 ) . (3.47) 

Corrections to it lead to enhancement by powers of log fj, with the maximal power increasing with the 
loop order. 

Our treatment of the fermions is far less satisfactory than for the bosons. The reason being that 
we unable to unravel the mirror transformation u — > v! in this case. Put differently, we did not realize 
how to achieve simultaneously purely imaginary values for Ef(u') and Pf(u') by transforming u to u' by 
a simple analytical continuation. The situation becomes however more transparent at strong coupling 
where the transformation can be accomplished, but only because the analytical properties 

in the u-plane of energy and momentum simplify substentially in this limit. It is found then that 
the fermion satisfies the double Wick rotation symmetry (|3.43p to all order at strong coupling. This 
strategy of considering the strong coupling regime as a governing principle also underlies our analysis 
of the the mirror transformation for bosons, as we shall illustrate for scalar below. Nevertheless, unlike 
the latter cases, the situation for fermions remains obscure beyond perturbative analysis, whence the 
exact analytical structure is restored. Fortunately the strong coupling analysis is very suggestive and 
allows one to conjecture an expression for the mirror Y function valid for any value of the coupling, 
perturbative or not. This will lead us to an educated guess for it which we shall employ afterwards 
for our weak coupling analysis. In this case, as well, we shall find that the representation (|3.44[) holds 
true. Our guess has also a bestowal property that £^ mrror (-u) = Ef(u) for any coupling. Then to put the 
fermion and boson on an equal footing we shall complement it with p]? lirror ('u) = Pf(u), implying that 
the double Wick rotation symmetry holds for fermions at any value of the coupling. Altogether this will 
allow us to evaluate the fermionic contribution to the Liischer formula. 

In the following, we shall make extensive analytic continuations of typical integrals defining P*(u) 
as found in previous sections. These entail crossing certain cuts in the complex rapidity plane. This can 
be done efficiently by using integral equations for the functions J±(t). With this application in mind, 
these equations are cast into a more convenient form 

[°°dt f 00 dt sin (ut) 7 (2gt) f°° dt sin (ut) J (2gt) - ut 
J J^(uth-(2gt) = -J o - ef _ x +2nlogr ? + 2y o — _ , 

"°° dt f°° dt (cos (ut) - J (2gt)) 1 (-2gt) f°° dt (cos (ut) - J (2gt))J (2gt) 
T cos(uth + (2gt) = -j o _ + 2^ — _ . 

(3.48) 

Both of these results are valid up to 0(l/r] 2 ) corrections and, more importantly, for u 2 < (2g) 2 . These 
equations (I3.48P can be deduced from the system (I2.44p . with the known homogeneous terms K n = 
k;^ log 77 + 5k^ + 0(l/ry 2 ), see Eqs. (|2.45| I2.46[) . after observing that the coefficients 7„ entering the 
Neumann expansion (|2.8|> admit the integral representation 







f°° dt r°° dt 

l2n = 2(2n)\ -J 2n (2gt) 1+ (2gt), l2n _ x = 2(2n - 1) / -J 2n . 1 (2gt) 1 _(2gt) , (3.49) 

Jo 1 Jo t 

for n = 1. 2. ... 1^1 These integrals then follow straightforwardly from the Neumann expansion (|2.8p and 
25 One has also 

r°° dt 

jJ (2gt)ry+(2gt) = Q. (3.50) 
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the orthonormality of the Bessel's function; 



After plugging (|3.49p into (|2.44|) . one easily arrives at the equations (|3.48|) by performing suitable 
linear combinations implementing the change of basis J n (2gt) — > (sin (ut), cos (ut)) with u 2 < (2g) 2 , see 
Ref. [HI [25] for more details. 

Notice finally that Eqs. (|3.48p imply that the two integrals 

r°° dt f°° dt 

/ -sm(uf) 7 _(2<7i), / -cos(ut) 1+ (2gt), (3.52) 
Jo 1 Jo 1 

can be analytically continued into holomorphic functions of u from the interval u 2 < (2g) 2 to the strip 
— 1 < 9m(u) < 1. This continuation is different from the one we would obtain starting from the interval 
u 2 > (2g) 2 and transforming the integrals to the form 

r^ S m(ut)^(2gt) = Ti f°° ^ e ±M 7 _(2 9 i) , [°° ^ cos (ut) 1+ {2gt) = f°° ^ e ±iut 7+ (2gt) , 
Jo 1 Jo 1 Jo 1 Jo 1 

(3.53) 

where the upper-/lower-script refers to u in the upper-/lower-half plane. The latter continuation indeed 
comes with a cut along u 2 < (2g) 2 while the former does not. There is no contradiction however, since 
the two integrals in (|3.52j) . which are both well defined for real u, are not analytical over the real li-axis. 
Though they define continuous functions of u, some of their derivatives jump at u 2 = (2g) 2 and thus 
the analytic continuation from inside and outside the interval u 2 < (2g) 2 can be inequivalent. 

We are now in position to attack the transition to the mirror kinematics on an excitation-by- 
excitation basis. 

3.2.1 Scalar 

We start our discussion with the scalar excitations. The transformation that supports the continuation 
to the mirror kinematics in this case is merely given by the shift u — > u + i. However this has to be 
done with due care since one encounters cuts in the rapidity plane along the way both in the energy and 
momentum of the scalar (see [JB]). To motivate the transformation alluded to above and understand 
how to deal with the issue of cuts, let us have a look at strong coupling where the dispersion relation 
becomes relativistic in a certain regime. The region in question is the low-momentum domain where a 
scalar excitation possesses the following energy and momentum 

Eh(u) = m cosh 8 + . . . , Ph(u) = m sinh 6 + . . . , (3.54) 

where the rapidity 9 = ttu/2 assumed to be fixed, while 5 > 1 and m ~ exp (— irg) is the mass gap [43] 
(see Eq. 13.1151 below for explicit expression). The ellipsis in (13.541) stands for corrections suppressed 
as m 3 at fixed 6, which incidentally break Lorentz invariancer 7 ! From the representation (|3.54j) it is 
clear that the double Wick rotation is obtained through the shift 9^9 + iir/2, as is well known. It 
corresponds to u — > u + i resulting in 

p™ irroi (u) = -iE h (u + i) = m sinh 9 + . . . = p h {u) , (3.55) 

26 We assume that we can commute summation and integration. This is not guaranteed by the convergence of the 
Neumann series (|2.8[l alone, even if it converges over the whole complex t-plane, which is the case here. The important 
point is that the expansion coefficients y n solving the system of equations (|2.44p are well suppressed at large n, for any 
given value of the coupling. 

27 The Lorentz invariance of the problem arises at strong coupling and at low energy only [43] . 
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Figure 3: Rapidity plane for energy and momentum of a scalar excitation. The double Wick rotation 
is implemented by the shift u — > u + i and depicted by the arrow. Composing two such rotations yields 
the particle to antiparticle transformation (E,p) — > (—E,—p) related to the crossing symmetry. The 
'mirror' line in the lower-half plane is associated to the transformation (E,p) — > (—ip,—iE). Up to 
the presence of the cuts, this is the same pattern as the one observed for energy and momentum of a 
relativistic particle in the #-plane, with 9 = ttu/2. 

and similarly for 22™ irror (u). A more general representation given in [38] that includes neglected contribu- 
tions in (|3.54|) shows that the same conclusions apply beyond the validity of relativistic approximation. 
Namely, the energy and momentum simultaneously achieve purely imaginary values under the shift 
u — > u + i when starting from the mirror kinematics and, as it turns out, exhibit the double Wick 
rotation symmetry (|3,43p . The latter representation, though still being limited to u 2 < (2g) 2 , suggests 
a generalization for any u: the rapidity u is carried over to the mirror kinematics by u — > u + i assum- 
ing it is performed by going in between the first two branch points in the upper-half plane located at 
u = i/2 ± 2g, respectively. The situation is illustrated in Fig. [H Notice that the path to the mirror 
kinematics closes off at weak coupling when the two branch points collide. This explains why it is im- 
possible to perform the double Wick rotation directly from the weak coupling expression for the energy 
and momentum. Not surprisingly, it opens up at strong coupling where this continuation is expected to 
make sense from the string theory perspective. This parallels, of course, the analysis for magnons. 

It is not difficult, starting from the general representation proposed in Ref. [IE], to verify that the 
identities in (|3.43p are fulfilled for the hole excitation with the energy and momentum defined in the 
mirror kinematics obtained via the procedure advocated above. The essential point here is to use Eqs. 
(|3,48p at intermediate steps for analytic continuation to higher imaginary values of the rapidity inside 
the strip — 2g < 3f?e(u) < 2g. Since we are interested in the transition to the mirror kinematics for the 
Y function itself rather than only for the momentum Ph(u), we shall explicitly illustrate the procedure 
for the quantity Ph(u), which also includes the anomalous component Sp^u) in addition to Ph_(u). Thus 
we define 

£ h (u) = -iP h (u + i) , (3.56) 
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with the shift being performed through the strip — 2g < 9?e(it) < 2g, see Fig. [3l 

To derive a convenient representation for £h(u) we can start from Ph(u) introduced in Eq. (|3.26p 
and perform the shift on the right-hand side of Eq. (|3.56p in two steps. First, we compute Ph(u + i/2) 
for u 2 < (2g) 2 . This is done by substituting u — > u + i/2 directly into (|3.26|) . since the latter integral 
representation initially covers the strip —1/2 < 9m(ii) < 1/2. This procedure yields 

f°° dt f°° dt 

P h (u + i/2-i0) =2i\og(j)/x) + -n-i / — e" m * j(2gt) + 2 / — sin (ut)j-(2gt) , (3.57) 

Jo t Jo t 

where we used the first equation in (j3.48j) . which produces the last term in the right-hand side of (13.571) . 
and also employed the explicit value of the integral 



00 dt 
T 



e~ lut J (2gt) 



-1 



ITT 

logx + ^(l)-y, (3.58) 



with u = x + g 2 /x and |x| > g. Note that — iO prescription on the left-hand side of Eq. (|3,57p remind us 
that this quantity is the boundary value (from above) of the original one given in (|3.26|) . It is important 
to stress however that the representation (|3.57p cannot be used readily if u is given even a small positive 
imaginary part, as the first integral in the right-hand side of (13.57P (would possess an exponentially 
increasing integrand and thus) requires proper analytical continuation. To remedy to it, we can use the 
fact that 

r°° dt f°° dt r°° dt r°° dt 

-i -e-»* 7 (2 9 t)+2 / -sin(nt) 7 -(2 9 t) = i / %»* 7 (-2 9 f)-2i - cos (ut) 1+ (2gt) , 
Jo 1 Jo 1 Jo 1 Jo 1 

(3.59) 

since j(t) = 7+(t) + J-(t), and then cast (|3.57p into the form 

f°° dt dt 
i£ h (u-i/2 + i0) = 2ilog(r)/x)+ir + i — e mt j(-2gt) - 2i / — cos {ut)j + (2gt) , (3.60) 

Jo t Jo t 

where by construction 

i£ h {u-i/2) = P h (u + i/2), (3.61) 

if u 2 < (2g) 2 . Now we can go to second and last step and complete the transformation by shifting 
u — > u + i/2 in (|3.60p . This is easily done at the level of the last term in the right-hand side of (|3.60p 
by using the second equation in (]3.48p . while for the first term we can use an identity similar to the one 
given before in Eq. (|3.58p . which reads 



L 



50 dt 

T 



e lut J (2gt) 1 



e*-l 



lo A (^)+V>(l) + y, (3.62) 



where u = x + g 2 /x and |x| < g. It is important to realize that this result holds for |x| < g since in the 
process we have indeed crossed the boarder \x\ = g and thus can no longer consider |x| > g as initially. 
Combining all results, we finally find that 



f' c 

£ h {u) = -2\og{g 2 /r,) +2^{1) - A 

Jo 



dt (cos (ut) e*/ 2 -J (2gt))J (2gt) - t/2 



t e* —1 

+ 2 r d2cos(ut)e^ 2 - M 2gt)^_ 2gt) ^ 
Jo t e — 1 

This is the sought expression that determines £h(u) to all loops, up to corrections vanishing at large 
spin. It is easily found in particular that 

£ h {u) = 2£ h (u)log7? + ... , (3.64) 
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after expanding (|3.63p with ^(t) = 7 (t) logr? + ••• and comparing the obtained expression with the 
general representation given in [IS] for Em (u) r 8 l Since, on the other hand, £h{u) = 2£ , ™ irror (u) log fj + . . . 
we verify our claim that E^ 1TTOI (u) = Eh(u). 

As we advocated before, the representation (|3.63|) is all we need for evaluation of the Y function at 
large values of spin and any coupling. The weak coupling limit is by far the simplest one to derive. It 
follows from neglecting the second line in (I3.63p . since 7 is small, while replacing the Bessel function 
Jo(2gt) in the first line with its Taylor expansion 1 + 0(g 2 ). This immediately yields 

r mirror (u) = ( _ 1} S gxp _^ {u) = (-1) 5 *ff + 0(g 6 ) , (3.65) 

?] z cosh [ITU) 

where we applied the integral 

f°° dt cos (ut) e*/ 2 -1 -t/2 , / 7T \ , « 

We note that the leading large-spin component of Y due to the energy E^(u) = 1 + 0(g 2 ) has produced, 
according to ()3.65|) . only the overall l/rf dependence in it at a given order at weak coupling. Most of the 
expression comes therefore from the anomalous component 5Eh(u). Notably, the latter is responsible for 
taming the large rapidity asymptotics in Y 1 [ nirror (n) and guarantees, as we shall see later, the convergence 
of corresponding integrals in the Liischer formulaic Finally, it also involves several factors of 't Hooft 
coupling, YJJ mrror (u) ~ g 4 , immediately implying that the finite-size energy starts at two loop order. 



3.2.2 Gauge field and bound states 

The continuation to the mirror kinematics is slightly more involved for gauge fields and their bound 
states. Here, in complete similarity to the hole excitations addressed in the previous section, it would 
also help to start the consideration with the strong coupling regime in order to figure out the correct 
continuation procedure. For the sake of generality, however, we shall present another suggestive path, 
which does no entail limitation to the strong coupling expansion. The transformation consists essentially 
of two steps. In the rapidity space, they combine to form a loop such that u — )■ u at the end of the 
day. Nevertheless, this requires changing sheets of the Riemann surface where physical observables are 
defined, thus making the transformation quite non-trivial. 

The first step is rather similar to the one adopted for the scalar. The energy Ei and momentum 
Pl for a bound state of I gauge fields both have their first branch points in the upper-half plane at 
u = ±2g + i£/2. In order to reach the mirror kinematics we need first to go through them, see Fig. HI 
Once this is done, we find ourselves on a new sheet where the energy and momentum both have a set 
of I cuts centered around the origin and spaced from one another along the imaginary axis by i. The 
uppermost cut, i.e., the one from where we entered to this sheet, lies at 9m(«) = i£/2, while the one at 
the bottom of the stack is at 3?m(ii) = —il/2. On this sheet, energy and momentum do not have more 
singular points than the ones associated with the cuts previously discussed, and, in particular, they are 
both analytic at infinity where they take the values 

Ee(u = 00) = £ , pi{u = 00) = i£ . (3.67) 

These expressions are exact at any value of the coupling)^. For the particular case of the gauge field, 
i.e., £ = 1, it is equivalent to the relation E g f(p = i) = 1 that was first deduced in [43j as a consequence 

28 The comparison also requires the use of the second equation in (|3.48|l evaluated at u = and restricted to the leading 
component ~ log 77. 

29 This contrasts it with a relativistic theory where the damping is produced by the large rapidity behavior of the energy. 
30 Note that have we decided at the very beginning to cross into the lower-half plane this would lead to the identity 
Ee(p — —i£) — I which is complex conjugate of Eq. (|3.67[1 . 
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Figure 4: Double Wick rotation for the gauge fields and bound states. 

of the Goldstone theorem for broken symmetry. We stress that the points u = oo involved in (|3.67p is 
completely different from infinity that could be reached by arriving to the real kinematics. The latter 
point indeed corresponds to the energy and momentum being real and taking arbitrarily large values^]. 

However at the moment we are only half way through since we still did not strike our final destination 
of mirror kinematics. As can be seen from (|3,67p the energy is essentially real here, while to perform 
the double Wick rotation we want it to be purely imaginary. Nevertheless, we made significant progress 
compared to our initial point. This can be illustrated by the following argument. We started with a 
particle at rest, in the real kinematics, which corresponds to Ei(u = 0) = me and pe(u = 0) = 0, with 
m,£ being the mass of the excitation. Then we transported these along the imaginary line u = i£, with 
£ real, up to the point £ = oo where ()3.67[) is fulfilled. Along the way, the energy was lowered from mi 
to £, since mi ^ £@ So the next step naturally presents itself and suggests that what we would like to 
accomplish upon completing our path is to reduce the energy further, i.e., all the way to Ei = 0. We 
can use the fact that the point u = oo is regular on the Riemann sheet where we currently reside and 
continue our continuation along u = i£ further with £ now ranging from — oo to 0. This way we will get 
to the desired location Ei = 0. 

In other words, in order to complete the transformation, the second, which is also our last, step 
consists in crossing through the cut lying at the bottom of the stack of i cuts present on the sheet 
alluded to above. This way we pass yet to another Riemann sheet, which is pretty much the same as 
the original one that we started from. On this sheet we find the mirror kinematics precisely at the same 
position as the real kinematics. Hence, as announced before, the transformation is a loop from u to u, 
up to the change of sheets. It is displayed in Fig. [4] 

What is remarkable is that once all the above steps are performed making use of analytical expressions 
for the energy and momentum, one comes to the same conclusion as for the scalar excitation, namely 

Ei{u) -> ipf VVOI (u) = ip t {u) , pi{u) -> iEf" ov {u) = iEi{u) . (3.68) 

Hence we again recover the double Wick rotation symmetry (j3.43j) . in agreement with [44]. The proof 
of (|3.68p is essentially the same as the one given earlier for scalars. 

31 At this infinity, one finds Ee(p — oo) ~ logp, a result first derived at strong coupling in [50] and generalized to any 
coupling in [48] . 

32 This inequality is saturated at zero coupling as then mj = I. At strong coupling we know that me = y2i, which is 
of course bigger than I. It is a bit of assumption that me increases monotonically from £ to y/2£ as the coupling grows. 
Though confirmed by the first correction at both weak and strong coupling, see [48] . the latter claim has however not been 
proved for generic value of g. 
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Since we established the precise form of the path in the complex plane that we have to follow to get 
from the real to mirror kinematics, we can start from the representation (|3.39p for Y^(it) with (|3.4(jp and 
continue it to Y£ mirror (it). Going along these lines we find that the latter quantity is of the form (|3.44p 
with the exponent given by 



£ e (u) = -2 log tf/rf) - 2(1 - 2)^(1) - 4 



+ 



2 / -(cos(^)e- fl / 2 -J (2 5 t)) 



° dt (cos (ut) e-^- 2 )*/ 2 -J (2gt))J (2gt) + (I - 2)t/2 

T e*-l 
1+ (2gt) 7 _(2#) 



l_ e -t e *-l 

(3.69) 

At large spin, Eq. (|3.64p derived for scalars is applicable here with obvious substitutions h — > £. 

The expression (|3.69|) can be easily estimated at weak coupling. As was done earlier for scalars, we 
can indeed neglect the second line in (|3.69|) and evaluate the first line using J (2gt) = 1 + 0(g 2 ). We 
get immediately that 

ymirror (n) = ( _l)^r(| + iufT^ - iuf + 0(g 6 ) , (3.70) 

where we used the following value for the integral 

2 r* . gWrH^Mg = log (r(| + m)r(| _ m)) _ (e _ mi) , (3J1) 

As expected, the Y function gets stronger suppression when the mass of the excitation increases, and 
for the twist-£ bound state, with mass = £ + 0(g 2 ), it scales as l/r] 2e ~ 1/S 2e . The dependence on 
the coupling, Y^ mirror (ii) = 0(g 4 ), is however the same for all excitations. Finally, we note that in the 
particular case of the gauge- field excitation, i.e., for t = 1, we find 

Yf™(u) = (~1) S 2 *X , + 0(g e ) , (3.72) 
?7 Z cosh (TTU) 

which is thus identical to the one for scalar in Eq. (|3.65|) . 
3.2.3 Fermion 

Finally, we come to the discussion of fermions. In this last case to be studied, the mirror transformation 
is the most problematic and our understanding of the procedure does not really go beyond the realm 
of the strong coupling regime. Fortunately it is all we need in order to make an educated guess for the 
mirror of If. So let us start with the strong coupling analysis. 

We shall consider the case of a small fermion, i.e., the one carrying the rapidity |x| < g. At strong 
coupling, it implies that the momentum of the fermion is of order O(l), i.e., p s { ~ g° for 5 > 1. This is 
the perturbative string regime where both the energy and momentum are found to admit an expansion 
in inverse powers of the 't Hooft coupling 1/g at fixed u = u/(2g). The expansion is well-defined for 
any u in the complex plane except for the branch points u = ±2g where the energy and momentum 
diverge. For illustration, we have [48] 

- i 1 \ 1/4 
U + 1 N ' 



E s{ + p sf = +..., (3.73) 

with the ellipsis standing for terms being suppressed by powers of 1/g but enhanced by the inverse 
powers of (u 2 — 1). What happens as we get too close to the points u = ±2g is that the strong coupling 
expansion needs to be resummed to all orders. If done properly, the singularities are resolved and the 
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Figure 5: Double Wick rotation for the fermion in the strong coupling regime. 

behavior around this points gets replaced by E s i(u ~ ±2<?) ~ p s i(u ~ ±2g) ~ g 1 ^ 4 . This is the transition 
to the near-flat-space regime that we shall not discuss here any further. In other words, we will avoid 
approaching too close to the branch points, which, as we recall, are the only singularities in the u-plane 
of the small fermion. 

The mirror transformation is easy to elucidate for the perturbative string regime of the fermion. 
This corresponds to a clockwise rotation around u = 2g, see Fig. [5j So, that similarly to the gauge field 
and its bound states, the path that we sought for is a loop u — )■ u and it turns out to yield 

E a f ->■ ip B f , pa -)■ iE si . (3.74) 

This is especially clear from ()3.73p where such a rotation leaves the result invariant up to the phase 
exp(i7r/2) = i. More generally it would not be difficult to show that the symmetry (I3.74p is satisfied 
to any loop order in the \jg expansion. We note finally that by passing four times around the point 
u = 2g, we will go through the sequence of transformations 

E af (u) -)• ip s i{u) ->■ -E s{ (u) ->■ -ipsf(u) -)■ E af (u) , (3.75) 

in agreement with the fourth order nature of the branch point u = 2g apparent in Eq. (13. 731) . 

Something goes wrong at this moment however because we know from the exact expression for a 
small fermion that the point u = 2g is a genuine square-root branch point, see Ref. [48] . Moreover, 
rotating once around it simply amounts to changing from the small to large fermion domain, 

E s{ (u) E u (u) , (3.76) 

as illustrated in Fig. [6j This cannot be a contender for the purported mirror transformation, since the 
outcome is not purely imaginary, and evidently contradicts (|3.74p . There is no real contradiction here, 
rather it is a predicament which is easy to unravel. Indeed, everything we said above concerned E B i and 
p s f at strong coupling, for which the natural variable is the rescaled rapidity, u = u/(2g). The function 
E\i(u) that we should reach at the end point of the loop has an infinite number of cuts in the upper-half 
plane, as opposed to E s [(u). At strong coupling and in the variable u, all these cuts condense into a 
single cut, the one with fourth order branch points apparent in (|3.73p . This was first elucidated in |24| 
in the analysis of the analytical properties of the solution to the BES equation. In our case, it means 
that we never really get to the sheet where the energy becomes described by En(u). Note that if we 
were aiming at attempting this, we would have to resolve the singularities at u = 2g first, thus restoring 
the proper cut structure and then navigate between them. 

The next question is where exactly does the mirror transformation land us at the end? This is the 
question we did not answer. From our discussion above it seems clear that the mirror transformation 
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requires a loop u — >• u crossing some of the aforementioned cuts. Unfortunately, all our attempts 
to deform the loop by crossing more cuts along the way have failed to produce a consistent mirror 
transformation. The fact that at strong coupling the infinite tower of cuts condense into one could 
imply that the exact mirror transformation cannot be obtained by crossing only a few of them. If it 
is the case one would have to figure out how to implement the continuation through an infinite tower 
of cuts, which goes far beyond our current abilities, see Fig. [6l The situation is however not entirely 
hopeless, because, after all, we have control on the mirror transformation at strong coupling. This means 
that we should be able to construct the mirror of Y s f (u) ~ exp iP s f(u) at least in the latter regime. It 
is worth a try because, as we shall see, it produces an expression for £$f(u) which is acceptable and 
meaningful at any coupling. 

Our consideration to turn P s f(u) into the yet mysterious £ B f(u) goes as follows. The starting point 
is the identity 

P h {u - I + iO) = P lf (u + iO) - P sf (u + iO) , (3.77) 

that relates Ph and Pif. In this equation everything is known except for the new quantity P s f(u + iO) 
that is defined as 2 

P sf ( u + i0) = -2ilog(^-) -TT-i [°° ^e iut 7+ (2gt), (3.78) 

VW Jo t 

where x is the small solution of the Zhukowski map u = x + g 2 /x, i.e., such that \x\ < g. The 
equation (|3.77p . with P s f as above, is a direct consequence of the expressions (|3.26p and (|3.32p for Ph 
and Pif, respectively. Now we can use the fact that the small and large fermion are connected for 
u 2 < (2g) 2 , i.e., that P l{ (u + iO) = P sf (u - iO), to find that 

P h (u -i + iO)= P s{ ( u - iO) - P si ( u + iO) , (3.79) 

when u 2 < (2g) 2 . This introduces P s {(u) into the game. What is interesting in the equality (I3.79P is 
that its left-hand side is exponentially small at strong coupling, as long as u stays away from u = ±2g. 
This is a simple consequence of the fact that, as we discussed before, the regime u 2 < 2g 2 happens to be 
non-perturbative for the scalar excitation at strong coupling, thus resulting in the asymptotic behavior 
Ph ~ exp (—ng). Therefore, at strong coupling, the relation (|3.79p simplifies to 

P sf (u - iO) = P s{ (u + iO) , u 2 < {2gf , (3.80) 

up to exponentially small correctionsjfl In other words, the two functions P s {(u) and P s {(u) are, at the 
perturbative level, analytically related to one another through the cut along u 2 < (2g) 2 . It immediately 
follows that after going around the loop u — >■ u, describing the strong coupling mirror transformation, 
we get 

P ri (u) -> P B{ (u) . (3.81) 

Indeed, according to (I3.80p . after passing in between the branch points u = ±2g, the function P s f(u) 
turns into P s f(u), to any order in the \ jg expansion. Thus the strong coupling mirror transformation 
for the small fermion leads to 

firf(w) = -»P rf («) , (3.82) 

with P B t(u) given in (|3.78p . If we are only interested in the image in the mirror kinematics of a fermion 
with real momentum, then we can restrict ourselves to u 2 > (2g) 2 and write 

/ a 2 \ dt 
S s{ (u) =-21og ( — J + »7T- / - cos {ut) 1+ (2gt) . (3.83) 
\i] x / Jo 1 

33 We emphasize again that these corrections grow however if u gets close to ±2g, such that we should avoid these points, 
in agreement with our previous discussion of the mirror transformation in the perturbative regime. 
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Figure 6: Continuation from the small to the large fermion domain is shown by the path going from 
the lower to upper bank of the cut. On the other hand, the path to the mirror kinematics is not clear 
away from the strong coupling regime. 

In a nutshell, what we have done by neglecting Ph(u — I + iO) in (|3.79p is precisely to remove an 
infinite tower of cuts depicted in Fig. [6] that was obstructing the mirror transformation. Thus we left 
the lower-half plane of the small fermion, which is free of singularities and controlled by P s f(u), and 
enter, after the above trick, into a 'new' upper-half plane with similar analytical properties governed by 
i£ s {(u) = P sf (u). 

Now let us emphasize an important point, which also happens to become our working hypothesis, 
that the expression (I3.83P is perfectly suitable for direct non-perturbative generalization. As it stands, 
it is perfectly well defined at any value of coupling and even obeys a very suggestive identity. Namely, 
when expanded at large spin, it leads to 

£ si (u) = 2E st (u) log fj + 0(log° r,) , (3.84) 

at any coupling, in perfect agreement with the putative double Wick rotation symmetry for the fermion. 
If we assume furthermore that p™ irror (tt) = p s {(u), which we know is satisfied at strong coupling, the 
situation becomes completely symmetric with the one observed for the bosonic excitations. Finally, it is 
interesting to note that a natural attempt to put fermions and bosons on an equal footing, as inspired 
by our findings for bosonic excitations, would allow us to conjecture that 

f°° dt 

£rf(«) = 2Iog»7 -/(«)- / -cos(ut) 7+ (2gt), (3.85) 

Jo t 

with f(u) being a spin-independent function of rapidity u. This ansatz would immediately warrant the 
validity of the large spin asymptotics (|3.84p . Our preceding analysis implies us that the strong-coupling 
mirror transformation yields 

/(«) = 21o g r^-*7r > (3.86) 

with x being the small solution to u = x + g 2 /x. Our proposal is to assume that ()3.86p is actually valid 
at any coupling, not only g — > oo where it was rederived. We shall see that it leads to very sensible 
quantitative results at weak coupling in complete agreement with perturbative expectations. 
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Now being equipped with the expression for £ s f(u), we find 11101 (u) via the formula 

Y s f TTOT (u) = -(-lf e - £ ^ . (3.87) 

One peculiar feature of l^™ irror (u), as compared to the bosonic case, is that it displays a singular behavior 
at u = oo, which corresponds to x = 0, i.e., the fermion at rest. Indeed, one easily verifies from (|3.87p 
and (|3.83|) that Y i miTTOr (x) ~ 1 /x 2 at small x. It originates from the term 2 log x in ()3.83p and is therefore 
justified at any coupling. The coupling-dependent piece is indeed controlled by the integral involving 
the function 7+(i) in (|3.83|) which is always regular at x = 0. This is apparent from the representation 







jCos(ut) 7+ (2gt) = 2Y^l2n(—) ~x 2 - (3i 
' n>l ^ a 



The above mentioned singular behavior of the y-function requires a prescription for integration at small 
momenta in the Liischer formula, as we shall see and elucidate in the following section. Finally, we notice 
the presence of in in (|3.83|) which after exponentiation compensates for the extra minus sign for the 
fermion present in (|3.87[) . This way the expression for Y™ irmT (u) supports the form (|3.44|) with E s f(u) 
and 5E s f(u) being real functions r 4 l 

The last step is to derive the expression £\{(u) for the fermion in the large rapidity domain which 
is the pertinent domain at weak coupling. This is quite straightforward to do and follows from the 
analytical continuation of Eq. (|3.83p to |x| > g region. Using the equations satisfied by 7(i), we easily 
find that 

r t„\ ^,„(S'\^' , rdt (c s(ut)-.M2gt))M2gt) 

Jo t e — i Jo 1 

where u 2 > (2g) 2 and with now x being the large solution of u = x + g 2 jx. It has a form very similar 
to the ones derived for the bosons, which is very encouraging. 

Finally, the evaluation of the mirror y-function at weak coupling can be done along the same lines 
as for the bosons. Namely, neglecting the second line in ()3.89p and replacing Jo(2gt) — > 1, we obtain 

y™ ( n) = (~1) S 2 .*X , + 0(g 6 ) , (3.90) 



after using the integral 



dt cos (ut) — 1 / ttu , 

~ log —tj-t , (3-91) 



o 



t e — 1 V sinh (ttu 



together with x = u+0(g 2 ). The previously observed small-x singularity manifests itself in the divergent 
behavior ~ 1/u 2 of (I3.90P at vanishing u, at this order of the coupling. Except for that, we note a striking 
similarity with expressions deduced by a rigorous analysis for the bosons, both in the dependence on 
the coupling ~ g 4 and in suppression at large rapidity. 

To conclude this section, let us add a comment on the mirror y-function for the small fermion. The 
latter does not play a significant role in our analysis at weak coupling. Rather than performing the 
entire analysis anew in this domain, one could instead find the result directly from ()3.90p . Namely, one 
first performs the small momentum expansion by using the relation p = 2u + 0(g 2 ) for the large fermion 
and then keeping the leading term at small p from (|3.90p and finally reexpressing it in terms of the 
rapidity variable of the small fermion p = 2g 2 /u + 0{g 4 ). This would be observed to match the weak 
coupling expansion starting directly from (|3.87p and (|3.83p . yielding Y™™ OT {u) = (—l) s u 2 /r] 2 + 0{g 2 ). 



34 These functions are real for u 2 > (2g) 2 
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3.3 Liischer formula 

We are now in a position to evaluate the Liischer formula at weak coupling. We recall that it controls 
the finite-size corrections to the energy of the twist-two operators. Its form was previously motivated 
in (|3.15p and argued to be given by 

£ FS = -e y ^ y h mirror (p) - 2 J £ y g f ror Cp) - s J ^ y f mirror (p) + — , (3-92) 

where the ellipsis stand for contributions from bound states of gauge fields and/or from subleading 
order Liischer corrections, both suppressed by higher powers of the spin (at weak coupling) @ The main 
ingredients entering this equation are the Y functions in the mirror kinematics 

n mim >), (3-93) 

for the elementary twist-one excitations, * = h, gf, f. These were constructed in the previous subsection 
to leading order at large spin and all orders of the 't Hooft coupling. The obtained asymptotic expressions 
are sufficient for our purposes and we do not need to worry about possible finite-size corrections to the 
Y's themselves. These effects should enter at the next order of the large-spin expansion and would 
require us to develop systematic TBA equations for their determination. 

Before we turn to the weak-coupling analysis of Eq. (|3.92p , some comments are in order concerning 
the choice of the integration contours there. When parametrized in terms of the momentum, it is clear 
that each integral in (|3.92p should be taken for real p to run from — oo to +00. For the scalar and gauge 
fields contributions the same integration contour is maintained as one passes to the more convenient 
rapidity u as the integration variable. However, the fermion case is a bit more tricky. The physical 
rapidity here is x and the contribution from the fermions reads 

E p = -8 f — ^ Y f milIor (x) , (3.94) 
J 2tt ax 

with x ranging from —00 to +00. But this integral should come with a prescription for integrating 
around the point x = on the real axis, corresponding to p = 0. Indeed, as we already observed, the 
function Y f mirror (x) is singular at small x and behaves as 

y f mirror (x) ~ 1/x 2 + 0(x 2 ) . (3.95) 

We choose to avoid this singularity by integrating around the point x = either by going slightly 
above or below this point in the complex x-plane. The two choices are equivalent; their difference being 
proportional to 

I dx^Yr mivvor (x), (3.96) 
J dx 

evaluated around x = 0, which vanishes since there is no residue at this point. The latter property also 
guarantees that the integral (13.94p with the =tiO prescription around x = is real. 

When written in terms of the rapidity u, the integral (|3,94p receives contribution both from the large 
and small fermion domains, corresponding to \x\ ^ g, respectively. At weak coupling, it is however more 
efficient to evaluate integrals involving quantities computed for large fermions. The reason behind it is 
that the large fermion domain covers the regime with p = 0{\) for g ~ which is more important than 
the one implied by the small fermion domain, corresponding to p ~ g 2 for g ~ 0. To re-express (|3,94p 

35 Higher Liischer corrections coming from the scalar sector are no longer small at strong coupling because the mass of 
the scalar becomes vanishingly small as g — s> 00. But as long as the scalar mass is of order one, which is the case to any 
given order in the weak coupling expansion, the first Liischer correction is suppressed by 1 /rf ~ 1/S 2 and subleading ones 
by extra powers of the spin. 
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in terms of Y 1 ]? lirror (w) and p\f(u) solely, we can deform the contour of integration in Eq. (|3.94[) such that 
it resides entirely in the large fermion domain > g. There are no obstructions encountered along the 
way since both pt(x) and l^ mirror (x) are analytical inside the disk \x\ ^ g, except for the double pole at 



x = in Y™ 1 ™" 1 



(x) which is bypassed. Then one finds the representation 



FS 



du dp\i 
2tt du 



V 1 



(3.97) 



where the contour of integration runs from u = — oo to u = +00 avoiding the cut along u 2 < (2g) 2 
by going either slightly above or slightly below it. This is the representation we shall use for our 
computation at weak coupling. 

Finally, let us note that the weak coupling expansion performs a kind of Laurent expansion in the 
integrand of Eq. (|3.97p which is valid away from the various cuts of p\f(u) and Y^? 1TIOT (u). In particular, 
the cut centered around u = 0, which we evade in (|3,97p . decomposes into contributions of the type 
~ g 2n /u 2n . As a consequence, at weak coupling, all contributions obtained by expanding the integrand 
in Eq. (|3.97[) in powers of g 2 are regularized by integrating with a ±z0 prescription around u = 0. This 
is how all integrals involving fermions entering the weak coupling expansion of (I3.97h will be understood 
henceforth. A similar prescription was adopted in [49 1 [68] when evaluating integrals relevant to the 
computation of gluon scattering amplitudes using the OPE for light-like Wilson loops [44] and which 
display an intriguing resemblance with the ones considered presently. 



3.3.1 Weak coupling expansion 

The result (|3.92j) for the finite-size energy correction admits an infinite series expansion at weak coupling 

FS _ ("I) 5 



E 



(3.98) 



where V n a polynomial in log rj of degree n — 2 with log r] = log r\ + 7e as before. This structure is an 
immediate consequence of the general form for the Y functions that we established earlier, 



(«) 



-l) s exp (-2E*(u) log fj - 25E*(u)) 



(3.99) 



with the scaling E*(u) = l + 0(g 2 ), valid for all twist-one excitations, explaining the overall suppression 
factor in E FS ~ l/?/ 2 ~ 1/S 2 . The weak coupling expansion of the exponent ()3.99p obviously generates 
powers of log fj increasing with the loop order and ultimately leads to the representation (|3.98p . 

We will report in this section on the explicit five- loop calculation of the finite-size energy (|3.98l) . 
with details deferred to Appendix |Dj In order to illustrate the problem at hand let us consider a 
particular set of contributions from the leading logarithmic effects. These are the coefficients in Eq. 
p.98p accompanying the maximal power of logry at a given order in 't Hooft coupling, i.e., 



c n _ 2 log n 2 fj + 



(3.100) 



where n = 2, 3, . . . . For these, the expression derived from the Liischer formula (|3.92p is rather compact 
and the emerging integrals are generic for its weak coupling expansion. The c n coefficients are generated 
by (|3.92p upon the substitution 



\r mirror 



V 1 
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(«) 
(u) 

(u) 



(-1)' 



(-1)* 



(-1)' 



rj 2 cosh (ttu 

^¥ 

rj 2 cosh 2 (ttu) 
rj 2 sinh 2 (iru) 



■ exp 



^g 2 {^{\+iu)+^{\-iu) 



exp 



exp 



4 5 2 (v>(f + iu) + v(! 



2V(1)) log fj 
2^(1)) log fj 



(3.101) 



4:g 2 (tp(l + iu) + V(l - iu) - 2ip(l)) log fj 
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together with the fact that = 2u + 0(g 2 ) for all of these twist-one excitations. Of course the 
form (|3.101|) does not produce correct subleading-log contributions beyond 0(g 4 ), however, this is 
irrelevant for the exploratory calculation we are after. Plugging (|3.10ip into Eq. (|3,92p and comparing 
with Q5J3SD , we find 

Cn = 6li n) +2/J ) +8/ f (n) , (3.102) 



where 
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n\ 
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-4) n 7T 
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+ iu) + tp(l -iu) 
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cosh 2 (ttu) 




(-0(1 + iu) + ^(| - iu) 


-2^(1))" 


cosh 2 (ttu) 




(?jj(l + iu) + V(l - iu) 


-2^(1))" 



4 n) = -— — I du 

/(«) = _ ±-2UL I du W 1 "*> 1 ^ ^ , (3.103) 



f n\ J sinh 2 (71-u) 

The contour of integration is chosen to be the real u-axis for the scalar and gauge fields, and similarly 
for the fermion but with a ±z0 prescription to go around the point u = 0, as discussed before. For any 
n = 0, 1, . . . the integrals (|3. 1Q3[) are well-defined and real, yielding a definite answer for the expansion 
coefficients c n . For low enough n, their evaluation is feasible but lengthy, see Appendix [D] for details, 
and eventually gives 

16vr 2 / 64vr 2 \ /448vr 4 256vr 2 

c = , Cl = — , pa = - 256C 3 + ^— , c 3 = - — — + 1024C 3 



3 ' V 3 J V 4 5 3 

(3.104) 

This finding suggests that the degree of transcendentality of the coefficients c n increases with the loop 
order, in agreement with general expectation for twist-two anomalous dimension. It is also interesting 
to point out the vanishing of the constant cq. This property is a consequence of the identity 

i r i 

du = tt = — du =- , (3.105) 

sinh 2 (ir(u ± iO)) 2 J cosh 2 (ttu) 

which is obtained by shifting the contour of integration to the upper-/lower-half plane, u — > u±i. Then, 
due to the exact balance between the total numbers of fermions and bosons the coefficient cq vanishes. It 
is not entirely clear to us whether this cancellation should have been expected or if it merely reflects the 
enlarged symmetry algebra of the spectrum of excitations at weak coupling. A similar phenomenon was 
observed for the (BMN) Liischer formula in Refs. |37tl38| where the cancellation occurs up to four loops. 
There, the effect relies on the supersymmetry of the BMN vacuum. In our case the cancellation takes 
place at the first order (i.e., two loops) only and is presumably restricted to the leading power corrections 
~ 1/r/ 2 that we are considering in our study. For contributions suppressed stronger by the spin, it would 
indeed be rather surprising to find that the cancellation persists, since bound states entering into the 
game do not possess super-partners thus distorting the equality of fermionic and bosonic degrees of 
freedom. Nevertheless the phenomenon of cancellation turns out to be rather generic among the various 
integrals that contribute to the polynomials V n in ([3.98P 

Summarizing, by taking into account omitted radiative corrections to Eq. (|3.10ip . we can immedi- 
ately complete the polynomials in ()3.98p . The analysis outlined in Appendix[D]provides the polynomials 



36 Indeed, the weak coupling expansion of y t mlrror (M) generates many periodic factors that simplify among bosons and 
fermions after using the contour deformation trick adopted above. However, it also generates combinations of derivatives 
of i/i-functions that are not periodic and do not lead to a complete cancellation, see Appendix [Ci and [Dl 
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up to five loops, 
V 2 = 0, 

^ 167T 2 , 

Vz = — log T] , 



P: - ( 256C 3 + log 2 fj + log fj + 16vr 2 C 3 



/448vr 4 _ 256vr 2 \ , / 896tt 2 C 3 1664vr 4 128tt 2 \ , „ 

n = - ( -^g- + 1024C 3 - — j log 3 r? + f 2560C 5 + ^ + -^g— + 256C 3 + — j log 2 f? 

,/ 1W(!/ j 8992tt6 2 \ _ ^ 320tt 2 C 5 896tt 4 C 3 ^ 

(3.106) 

An attentive reader would notice the absence of ~ log ?] contribution in V3. Again, this stems from 
the cancellation between the bosons and fermions in complete similarity to the identity V% = cq = 
examplified before. 



3.3.2 Finite-size vs wrapping corrections 

We are now in position to compare the finite-size corrections to the GKP vacuum, we are advocating 
here, and the wrapping corrections to states built on the BMN ground state. We first recall what we 
have learnt and what we aim at proving. 

On the one hand, we have shown that the twist-two ABA energy can be decomposed as 

^ABA = ^bulk + E elt ^ (3.107) 

where i? alt is the alternating component of the ABA energy that scales as i? alt ~ (— \) s /rj 2 at large 
spin, up to powers of log 77 in the numerator, and so to any order at weak coupling, see Eqs. (|2.73|) and 
(I2.74p . The remaining component, E hnlk , does not alternate with the parity of the spin S and has a 
much softer large-spin expansion, given in Eq. (|2.60p up to 0(l/r? 2 ). On the other hand, motivated by 
the picture of the twist-two background as defining the ground state of the (long) GKP string, we have 
constructed the finite-size energy E FS valid at leading order in large spin in the form of the Liischer 
formula (13.92j) . We have seen that this quantity alternates with S mimicking S alt and its evaluation 
at weak coupling have shown that it scales as E FS ~ (—l) s /rj 2 , and, more generally, has the structure 
similar to the one observed for E alt , see Eqs (|3.98p and ()3.106p . 

As alluded to in the Introduction, our prescription for correcting the ABA predictions _E ABA and 
recovering the complete twist-two energy E is to perform the substitution 

E alt E FS , (3.108) 

such that one replaces (13.107P by 

E = E hulk + E FS . (3.109) 

Though we do not have a compelling physical reason at weak coupling why this is the proper procedure, 
the similarity of the alternating ABA energy and finite-size corrections, together with the evidence 
that the ABA result has to be corrected, makes the substitution (|3. 108|) as a working hypothesis. We 
will later provide more physical evidence in favor of (|3,108p coming from the strong coupling analysis. 
However, first we will check that our substitution rule passes quantitative tests at weak coupling. At 
a qualitative level, one could regard the sum (I3.109P as a decomposition of two different classes of 
contributions, shown in Fig. [TJ The bulk energy should resum those corrections associated to diagram 
with trivial topology as in Fig. Q] (d). The ones controlled by the finite-size energy are associated 
with graphs possessing non-trivial topology, with the first Liischer correction being associated to those 
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with excitations wrapping once around the GKP string, see Fig. [T] (c). Since our understanding of the 
decomposition (|3.1U9|) is so far restricted to the leading contribution at large spin, it should hold true 
up to corrections suppressed stronger than 1/rj 2 ~ 1/S 2 to any order in weak coupling. 
We notice that as a corollary of (|3.1(J9|) we should recover the identity 

E FS = E &lt + E WT£Lp , (3.110) 

since by construction the wrapping corrections complete the ABA result, E = P ABA + ^wrap_ Q ne Q £ 
the main results of this paper is to establish that the relation (|3.110|) is indeed correct up to five loops 
at weak coupling (allowing us to extrapolate this claim to any order in g). One might be concerned 
that though both E and P alt alternate with spin this turns out not to be the case for E wrap , as we 
shall see. (This is the case even for the leading order contribution to i? wrap in the large-spin limit.) 
However, this fact is not surprising at all if one recalls that only even values of the spin are physical, 
i.e, corresponds to state with zero quasimomentum consistently described by the long-range ABA. The 
conclusion therefore is that (|3.11(J|) should be obeyed for even spins only. The fact that the extension to 
odd spins proposed here does not match the one realized in the construction of the wrapping corrections 
should not be so puzzling after all. The extension to odd spins of the (BMN) Liischer formula of Ref. 
|38j was done under special requirements that do not seem to have a simple physical interpretation in 
our picture. Reciprocally, it appears that odd spin configurations are somehow associated, in our case, 
to a change in the boundary conditions imposed to excitations propagating around the GKP string; a 
property which bears no similarity with the steps adopted in |38j. 

Keeping this caveat in mind, let us now come to the quantitative test. In to order to simplify the 
comparison, let us introduce a representation for the wrapping corrections similar to the ones found for 
the alternating and finite-size energies (12.731 [3798]) . Hence, we set 

^wrap = ^ P™rap g 2n + ^/^ , ( 3 .1H) 

^ n>2 

where P™ rap is a coefficient (weakly dependent on spin) emerging from the large-S expansion of the 
n-loop wrapping energy that governs its l/?7 2 -scaling behavior. Testing Eq. (|3.11Qj) is then equivalent 
to verifying 

V n = P n + P^ p , (3.112) 

with V n and P n known explicitly from Eqs. (|3.106p and (|2.T4[) up to five loops. We first recall that 
since both V n and P n are polynomials of degree n — 2 in log f], the same should be the case for i-*™ rap , at 
least for even values of the spin. We will see that this is indeed the case up to five loops, independent 
of the parity of the spin, as was already elucidated at four loops in [39] by the large-spin expansion of 
available wrapping corrections obtained in Ref. [38j . 

The two particular cases of (I3.112j) yield V2 = Pi and V3 = P3, since the wrapping effects start 
at four loops for twist-two operators, i.e., P^^ = 0. Comparing (I2.74p and (I3.106P we see that these 
relations are indeed satisfied, and what is more they are correct for both even and odd spins. To 
proceed further we need p™ rap an d p™ rap . Starting from the exact (finite-spin) expressions for the 
four-loop and five-loop wrapping corrections to twist-two anomalous dimension found in Refs. [38j and 
[42] . respectively, we derived the following expressions, 

pr p = - (128C3 + ^fV) iog 2 ^, 

^5 ( + ( 1024 Cs — )v) !og V 

( A 512vr 2 C3 /1664tt 4 128vr 2 \ \ 2 , 2 2 , 

+ ( 1280C5 + 3-^ + + 256C3 + —3— J a\ log 2 fj + (768C 3 2 + 128^3^) log 7? , 

(3.113) 
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with the signature factor a = (— 1) . As expected they are polynomials with correct degrees of log 77 at 
a given order in 't Hooft coupling. For even spins, i.e., a = 1, we checked the agreement at four loops 
with the large spin result of [39] , since log fj = log S = log S + 7e at the given accuracy. The five-loop 
result seems to be new, however. What is important here is that, after plugging (|2.74|) . (|3.1U6|) and 
(|3.113p into Eq. (|3.112j) for n = 4, 5, we immediately verify the validity of (|3.110p up to five loops 1^1 
We stress that the equality holds for a = 1 only, indicating a discrepancy between our approach and 
the standard one [38] with regards to the extension to the (unphysical) odd values of the spin. It is 
interesting to note that this discrepancy does not seem to affect those contributions accompanying the 
maximal transcendentality at a given loop order and power of log ry, since these are precisely the terms 
independent of a as can be verified by a simple inspection of Eq. (13,113!) . 

As demonstrated in this section, we succeeded in reproducing the exact twist-two energy at large 
spin, including the wrapping corrections not captured by the ABA equations. The five-loop agreement 
demonstrated above gives strong support to the correctness of our proposal for the Liischer formula 
(|3.92p and to the accompanying substitution rule (|3.108j) . Though our result only captures the leading 
large-spin asymptotics, it provides us with a transparent physical understanding of the structure of 
the large spin expansion beyond wrapping order. This is of course a consequence of the fact that all 
one needs in the description is the solution to the cusp equation and its siblings [27] controlling the 
subleading ~ log fj contribution to the large-spin distribution of the Bethe roots. Both of these are 
easily constructed at weak coupling. One can also analyze these at strong coupling, and, in the next 
subsection, we will use them to suggest another check of the validity of the Liischer formula (|3.92l) . 



3.3.3 Remarks on the strong coupling limit 



In this subsection, we will consider the Liischer formula (|3.92j) at strong coupling. Its evaluation is 
slightly more demanding than at weak coupling, because each excitation develops several scaling regimes 
at strong coupling. Namely, the full kinematical range (assuming the reality of the momentum p) 
typically splits in various domains that chime in with those found for magnons at strong coupling, e.g., 
one can find a perturbative regime with E,p ~ g° which is typically probed by the perturbative string 
methods [14J and a semiclassical one with E,p ~ g that can be analyzed using finite-gap theory [50] or 
strong-coupling methods of Ref. [48J. 

For the present analysis, the domain of interest is the one that minimizes the energy. This is easily 
understood by noticing that at large spin each of the individual integrals in (|3.92p can be evaluated by 
means of the saddle point approximation. The dominant contribution therefore comes from the regime 
where E+{p) is minimal, i.e., with p around p = such that E+(p) = m+ + 0(p 2 ) with m* being the 
mass of the ^-excitation. Each integral in (I3.92|) thus scales as 

Ef 3 ~ l/r/ 2m * ~ l/S 2m * , (3.114) 

up to some (fractional) powers of the length ~ 2 log rj. 

It follows from the estimate (I3.114p that the leading contribution to the Liischer formula (|3.92p 
comes from the 6 scalars. This reflects the fact that the scalar mass becomes vanishingly small at strong 
coupling being exponentially suppressed with g [4"3 j 1691 16T] . 



2 l/4 A l/8 



1 + OQ./VX) 



(3.115) 



with VA = 47r g 3> 1. This phenomenon was explained in [33] where the 6 scalars were identified at 
strong coupling with the multiplet of massive excitations of the 0(6) sigma model, which encodes the 



As a side remark, we notice the absence of 0(log fj) terms in both p™ rap and p™ rap . it the latter observation extends 
to all loops, and if the equation (|3.110[) is also correct to all loops, it would turn into a non-trivial equality between E FS 
and -E alt for this kind of corrections. 
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low-energy effective theory of excitations on the long GKP string. In distinction, the masses of fermions 
and gauge fields remain of order ~ g° at strong coupling |14|, H5] , more precisely m,f = 1 and m g f = \/2, 
yielding corrections to the energy suppressed at least as 1/S 2 and thus negligible compared to the 
contribution of holes. 

So let us focus on the scalar sector from now on and let us evaluate more accurately the Liischer 
formula, this time at strong coupling. For scalars, the domain connected to the point p = corresponds 
to the non-perturbative regime where these excitations become relativistic and their dispersion relation 
reads [EH SE] 

Eh(u) = fn cosh 6 , Ph(u) = msinh# , (3.116) 

where = ttu/2 and the mass was introduced in Eq. (13. 115ft . The result (|3.116p is an approximation 
to the exact dispersion relation that holds at strong coupling for u ~ 1 up to corrections suppressed 
by higher powers of m. To compute E FS we also need the expression for the scalar ^-function in the 
the mirror kinematics at strong coupling. Actually we need the expression for SE^u) since all other 
ingredients involved in Y^ nirror (ii) are known. By using the method developed in Ref. [61J, one can 
construct the latter and then to demonstrate that the complete expression reads 

ymirror^ = g-mHooahfl = (_ e -E h (u)R ? (3 m ) 

where R is a parameter interpreted as the effective length of the 0(6) sigma model. Up to corrections 
explicitly suppressed by the spin S, the effective length admits the strong coupling expansion with the 
dominant term behaving as 

fl = 21og^£) +0(1/VA). (3.118) 

What we immediately conclude from Eqs. (|3. 117[) and (|3. 1 18|) is that both 5Eh(u) and Eh(u) are 
proportional to cosh#. We should stress however that this behavior is restricted to this particular 
kinematical window where the 0(6) description applies. It is lost once we include subleading non- 
perturbative corrections ~ m 3 which originate from irrelevant deformations of the 0(6) sigma modell 3 ^! 
Now, plugging the expressions (|3.117p and ([3.116P into the Liischer formula (|3.92p . we find 

E FS = - 6m( 2 ~ 1)S J d6 co S h(6)e- mcosh ^ R = - Qm{ ~ 1)S K x {mR) , (3.119) 

where K\{z) is the modified Bessel function, with the well-known asymptotic behavior K\{z) ~ exp (—z) 
asz-> +00. The equality (|3.119p holds up to non-perturbative corrections ~ m 3 coming from the scalars 
(and associated to irrelevant deformations of the 0(6) model) and up to ~ 1/5" 2 , and higher suppressed, 
contributions due to the fermions, and gauge fields. Finally and more importantly, we also recall that 
it only takes into account the first Liischer correction, associated to the virtual propagation of a single 
excitation around the cylinder, see Fig. [1] (c). 

Summarizing, we find that the energy is given (for even spins) by 

E = £ bulk - —K^mR) + ... . (3.120) 

This is the expected result for the energy of the vacuum in the 0(6) sigma model defined on a cylinder 
of circumference R with periodic boundary conditions |52| . We note that from the low-energy point of 
view the bulk energy is a non-universal quantity that depends on the embedding of the 0(6) model into 
the superstring sigma model. Were the latter model not UV finite, the bulk energy would be associated 



38 So only after including these subleading effects do we recover that we are effectively dealing with a non-homogeneous 
background - an intrinsic virtue of the GKP background otherwise observed in the perturbative string theory descrip- 
tion [MET]. 
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to the UV divergent contributions to the vacuum energy that typically scale with the length to leading 
order at large length: E bulk ~ A 2 R with A some UV cut off. In most treatments, this quantity is 
discarded and the vacuum energy is substracted such that E vacunm = E FS . This is the case if E vacunm 
is computed by means of the TBA equations for the 0(6) sigma model, for instance, where the bulk 
energy is simply zero as it is the one used in the Bethe-Yang approximation. In our case, £^ bulk is given 
as an expansion in the inverse powers of 1/rj with expansion coefficients being function of the coupling. 
Since we discarded in (|3,120p corrections suppressed stronger than 1/rj 2 for r\ ~ S S> 1 only the leading 
and subleading contributions to E huFk should be kept in (|3.120p . One finds for it |22^ 127] 



•i.uik v 7 !, f8nS\ \ A 



^buik = j_ bg _ 2_ + (log S/S) , (3.121) 

to leading order at strong coupling, in agreement with the string theory semiclassical result Q3] • 

Since the finite-size energy is always suppressed at large spin, independently of the magnitude of 
the coupling, weak or large, the first two terms in the large spin expansion of the bulk energy, E hnlk = 
2A log S + B + (log S/S), always provide the leading contributions to the total energy E at large spin. 
Both functions, A and B, are well known for any value of the coupling [I7|l27|l22MI2itl2Sll2gtlI3|II5]. 
The formula (13.120P predicts that the next term in the large spin expansion of E is actually controlled 
by the finite-size energy at strong coupling and scales as 1/S 2m with m ~ exp (— irg). This conflicts 
at first sight with the perturbative string analyses [301 ED E2] that detect a correction suppressed as 
1 / log S at one- loop level — a prediction that was confirmed using integrability in [33] . We will come 
back later to this apparent discrepancy between our result and the one stemming from string theory, 
which as we shall see turns out to have a simple explanation, but first will compare the result (I3.120j) 
to the ABA prediction at strong coupling. 

We have seen that our formula for the finite-size energy correctly reproduces the expectation from 
the 0(6) sigma model. This agreement relies on the substitution f|3. 108 j) to reconstruct the exact energy. 
To shed some light on the physics behind this substitution let us compare the finite-size energy (I3.119j) 
with the prediction stemming from the ABA equations. To do so, we replace E FS by i? alt in (|3. 109|) in 
order to get back to the ABA energy (|3.107p and then evaluate the alternating ABA energy using the 
representation (|2.75p . To compute the integral in (I2.75P at strong coupling we first deform the contour 
of integration in (|2.75p and derive an alternative representation which reads 

E &lt = ~ j> d Ph {u - §)Y h mirror (u - i) . (3.122) 

Here the contour of integration goes anticlockwise along the cut u — i/2 G [— 2g, 2g] present in both 
Y h mirror (u) and Ph(u). The fact that (|2.75p and (|3.122p are equivalent is quite straighforward once we 
recall that the mirror kinematics is reached (for scalars) by analytically continuing through the cut 
involved in ()2.75p . Putting it differently, one can pass from the representation (|2.75p to ()3.122p by going 
up and down through the cut. Crossing the latter, one finds 

Y h (u + §) ^ Y h mirror (n - §) , Eh (u + i)^ip h (u-i). (3.123) 

The extra minus sign generated along the way compensates for the change of the circulation of the 
integral such that both contours (|2.75p and (|3.122p are traversed in the anticlockwise direction. 

Before computing (|3.122|) at strong coupling, it is interesting to note that the representation (|3.122|) 
gives a more transparent understanding of the large spin structure of i? alt at weak coupling. Indeed, 
up to the choice of the integration contour, the integral (|3.122|) is remarkably similar to the Liischer 
integral for scalar excitations. For this reason their structure are the same and E alt and E FS admit 
similar expansion at weak coupling (12.73p . (I3.98p . though the expansion coefficients obviously differ. It 
is in particular quite straightforward to derive a close formula for the leading logs in the weak coupling 
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expansion of i? alt . It follows from the same kind of approximation as for the Liischer formula, namely, 
by performing the substituting 



»^(~1) 5 ffi , exp 
cosh ( ITU ) 



4g 2 (^(i + iu) + V(i - »«) - 2-0(1)) log r? 



(3.124) 



into (|3.122p together with Ph(u) — > 2u. 

Getting back to the computation of (|3,122p in the strong coupling regime, we now deform the contour 
of integration in (|3.122j) in such a manner that the integral splits into two, 

E ^ = -- f dp h (u)Y^" m (u) + - f d Ph (u - i)y h mirror (n - i) , (3.125) 



7T / 7T 



where in both cases u runs from — oo to +oo. Since the right-hand side of (|3.125j) originates from the 
contour integral in (|3.122p it should be clear that the shift u — > u — i appearing in the second term 
in (]3.125p is performed outside of the cut [— i/2 — 2g, —i/2+2g] where both functions Ph(u) and y^ nirror (?x) 
possess discontinuities. This shift is therefore different from the one considered previously when the 
transition from the real to mirror kinematics (for scalars) was discussed. In particular, ph_(u — i) ^ 
—iE^{u) and Y^ mrror (u — i) ^ lh("u)13 These remarks are important because they clarify why the two 
integrals in (13.125P have a radically different scaling at strong coupling, the second one being much more 
suppressed as a function of the spinF^I 

To estimate the two integrals in (|3.125p we can use the saddle point approximation as before. The 
first integral is actually identical up to a numerical factor to the scalar contribution to the finite-size 
energy. We already know therefore that it scales as 1/ r] 2m ~ 1/ S 2m at large spin with m given in (I3.115P . 
For the second integral we use the fact that the scaling at large spin is controlled by E^(u — i) asu-s-O 
and apply the recurrence relational 

E h {u -i) = E g{ {u) - i( Ps{ {u + i) + Psi {u - £)) , (3.126) 

which is easily derived from the representation obtained in [35]. Note that the second term in (|3.126p 
is purely imaginary since the momentum of a small fermion is odd in u, i.e., p s f(—u) = —p s f{u), while 
the first term is real and even. Now since E^—i) = E g f(0) = m g f where m g f is the mass of the gauge 
field and because m g f = \[2 to leading order at strong coupling, we deduce that the second integral 
in (|3.125p is suppressed as 1/r/ 2 ^ ~ l/S 2 ^ 2 at large spin, which is thus negligible compared to the first 
contribution. 

Putting all pieces together, we conclude that the alternating ABA energy is given at strong coupling 

by 

E M = - 2 (~ l > K x (mR) + ... . (3.127) 

Comparing this result with the finite-size energy (|3.119p . we observe that the ABA prediction does not 
properly account for the right number of scalar excitations. The ABA energy seems to detect only 2 
out of 6 scalar excitations of the GKP string. This might not be totally surprising when we remember 
that the ABA equations were designed for the BMN vacuum which singles out 2 scalars out of the 6 and 
breaks the 0(6) symmetry down to 0(4) @ To some extent, the result (I3.127P and its interpretation 



39 This is easily visualized in Fig. [3]where the cut separating real and mirror kinematics is drawn outward and where the 
expressions involved in the second integrals in (|3.125[1 are obtained on the second sheet of the Riemann surface, which is 
reached from the upper-half plane by passing below the cut. This obviously does not correspond to getting to the real line 
that lies on the same sheet as the mirror one on Fig [3] 

40 This contrasts with the weak coupling regime where the two integrals yield comparable contributions. 

41 With the shift u — > u — i performed outside of the strip — 2g < 5fte(it) < 2g. 

42 These two scalars are the fields Z and Z, which are the complex conjugate of one another, with Z being the building 
block of the BMN vacuum Tr Z L and Z being a mass 2 excitations typically unstable in the latter background. 
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gives more credit to our substitution rule f|3. 108|) . Stepping away from the g — > oo limit, i.e, at lower 
value of the 't Hooft coupling, one has to include fermions and gauge fields into the game, so tha the 
replacement (|3. 108j) can be interpreted as restoring the proper account of fluctuations on top of the 
GKP string. 

Finally, let us come back to the difficulties previously encountered when comparing (|3.120p with the 
perturbative string theory computations [3Q |, I31 |, [32] . It is clear that the formula we derived in (13.120P 
for the twist-two energy at strong coupling is essentially non-perturbative and thus cannot be compared 
directly against any perturbative string predictions. The point is that our result (|3.120p holds at 
strong coupling iff mR 3> 1 while the string perturbative analysis always assumes that mR <C 1. The 
interpolation between these two regimes cannot be achieved by accounting only for the first Liischer 
correction alone: we recall indeed that the (complete) finite-size energy should include also higher 
Liischer corrections, which incorporate the multi-particle exchanges of virtual scalars, in particular. 
These should scale as 1/S 2nm (with n > 1) at large spin and they are all therefore of the same order of 
magnitude as the leading Liischer correction (|3.119p provided mR ~ 2m log S is not big enough. This 
can only happen at strong coupling when m becomes arbitrarily small since otherwise mR ~ R ~ 2 log S 
which is automatically large. As our consideration originated at weak coupling where m ~ 1, it was 
quite natural to neglect those higher Liischer corrections because they enter at order 1/S 4 and higher, 
and are hence much more difficult to analyze in our framework. What we observe here is that this 
hierarchy gets lost when the coupling gets strong enough due to the quasi-masslessness of the scalar 
excitations. On the other hand, on the string theory side, the perturbative analysis always proceed by 
expanding in strong coupling prior to taking the large spin limit. It means that the strong inequality 
mR 3> 1 cannot be met since m is exactly zero to any order in the strong coupling expansion!^! All of 
what we have just said is in accord with standard observations once we adopt the low-energy perspective 
and think in terms of the 0(6) sigma model. Then, the first Liischer correction is good approximation 
in the massive regime corresponding to mR S> 1, while the perturbative analysis applies in the UV 
regime associated with mR <C 1, that is, when the Compton wavelength of the massive excitations can 
be regarded as infinite compared to the system size. In the latter 'phase' the correlation length is then 
infinite and the finite-size corrections exhibits the scaling ~ 1/R ~ 1/ log S typical of massless theory, in 
agreement with the stringy perturbative predictions |30 |. I31| . [32 | I33j . The details of the transition require 
the full-fledged TBA equations for the 0(6) sigma model, analyzed in [52], but some of its consequences 
on the structure of the large-spin expansion on the stringy side can be derived by renormalization group 
arguments only. These were already given at the end of the introduction, and we refer the reader there 
instead of repeating them here. 

These last remarks conclude our analysis of the subleading large-spin corrections to the twist-two 
anomalous dimension. As we have established, the weak-to-strong transition is to a large extent due to 
the non-trivial dynamics in the scalar sector, which is responsible for the change of scale 1/S 2 — > 1/ log S 
as one flows between the two extremes, i.e., the weak- and strong-coupling regimes, characterized by 
g 2 <C 1/log S and g » ^ log log S, respectively. 
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43 For the string theory calculations mR = which we interpret as saying that the perturbative analysis is valid for 
mR < 1. 
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Figure 7: The sketch of the complex u-plane for the function f(u + ). The crosses denote the square-root 
branch points with a cut running between them. 

A Comment on singular part of the function t(u) 



In this appendix, we give more details on the derivation of the expression (|2.22p for the non-polynomial 
contribution, s(u), to the function t(u) appearing in the right-hand side of the Baxter equation (11.131) , 
We discuss in particular the choice of the integration contour used in the integrals entering Eq. (|2.22p . 

To start with, let us recall that the function s(u) is given in ([2.18]) as a series in inverse powers of 
the shifted Zhukowski rapidities = x(u ± i/2), which we reproduce here for convenience, 



n>l 



w_ 

x+ 



n>l 



9 

ix~ 



(A.l) 



The expansion coefficients s n ,s n , are nonzero due to singularities of the dressing factors A±(u) in the 
u plane. Indeed, when regarded as a perturbative expansion at weak coupling, these factors develop 
poles at u = of increasing order with each power of the 't Hooft parameter g 2 . This is the origin for 
the introduction of the s(u) component of the transfer matrix to compensate non-polynomial effects of 
both sides of the long-range Baxter equation. Therefore, the coefficients s n ,s n are chosen in a manner 
for the sum (|A.ip defining s(u) to be a collection of the polar parts in the Laurent expansion around 



u 



= ±| of the left-hand side of the Baxter equation (|1.13|) . Let us see it in more details. 
Consider a function / that admits a Laurent expansion around u + = u + i/2 = 



/(« + ) 



fn 



n>l y ' 



+ regular . 



(A.2) 



In the particular case under study, we are interested in resumming the polar parts of the Baxter equa- 



tion (|1.13p around u 
can be chosen as 



±4. We shall treat the analysis around u 



f{u + ) = A + (n+) 



Q(u + + |) 



i/2 as a specific example. Then / 

(A.3) 



as this is the only source of singularities around u = — | on the left-hand side of the Baxter equa- 
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tiono For self-consistency of the Baxter equation, the singular part of f(u + ) ought to be exactly 
counterbalanced by s(u) on its right-hand side, implying that s(u) should contain the sum 



s it 



/(«+)!. =y-^-. (A.4) 

J ^ rising / / (n + ) n 
n>l V ; 



Analogous contribution comes from singularities at u = |, with obvious substitutions in all formulas, 
and s(u) is the sum of the former and latter. To get to the expression 12.221 we need to re-expand s(u) 
in the basis of inverse powers of x + rather than u + . To perform this step let us first observe that 

fn = ^-i-v n f{v), (A.5) 

where the integration contour encloses v = in the counterclockwise direction. Then, after plug- 
ging ()A.5p into ()A.4p and interchanging summation and integration, we find 

K^) = ^i-^—m. (a.6) 

Znr J vr — v 

Next, we proceed to the rapidites x + and y via the Zhukowski maps 

X+ +g 2 /x + =U + , y + g 2 /y = v, (A.7) 



yieding 



«(«) = ^-<f — J v 2/ + r (l " 9 2 lv 2 )f(v) . (Ai 

2m J (x+ - y)(l - g z /x+y) 



Its expansion in inverse powers of x + is straightforward and reads 

^2 \ n i 



y 



f{y) . (A.9) 



This is the sought expression, in agreement with the result quoted in (|2.22 j) for the coefficient s n \ | 

We note that the contour of integration should be defined such that it encloses all possible singu- 
larities of f(u + ) emanating from u = —\ as the coupling increases. The poles condensing at this point 
transform into a square-root cut [— | — 2g,— | + 2g] associated with the Zhukowski map for finite 't 
Hooft coupling. Thus the integration path should encircle this cut in the rapidity plane, see Fig. [71 It 
means that the contour C in the y integral ()A.9p should encompass the disk \y\ ^ g. Finally, there are 
potentially extra singularities stemming from the series ~ X^n>i ln{g / x+ ) n defining the dressing factors 
Aj r (u + ). In the weak coupling expansion, these singularities are confined to the second sheet |x + | < g. 
But after resummation, some singularities, like extra poles, can migrate to the first sheet and should 
be encircled by C in the integral (|A.9p . This does not happen however in the case we are discussing] 
such that the contour of integration C in ()A.9P can be chosen to be the circle \y\ = g, for instance, at 
all time. 



44 For convenience we divide both sides of the Baxter equation (|1.13[1 by Q(u). We also assume that all the roots of the 
latter polynomial are real, which is the case in the sl(2) sector we are considering here. 

45 A similar analysis around u — i/2 would provide the representation for the coefficient s„ given in (|2.22[) . 

46 In our case, indeed, these singularities appear in the form of an infinite sequence of cuts in the domain \x + \ < g. They 
occupy the same position in the associated u plane at any coupling and thus never get to the sheet with \x + \ > g. These 
conclusion can be easily drawn from the analysis of |24] , 
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B Twist-one excitations in ABA 



Let us briefly recall the identification of elementary twist-one excitations emerging in the consideration 
of degrees of freedom propagating on the long-range spin chain. The latter can be encoded in the 
Dynkin diagram displayed in Fig. [5J Each node is associated with Kj Bethe roots. These are excited 
by the acting Kj times, respectively, with step-up operators of the Serre-Chevalley basis of the 
superconformal algebra on the vacuum state |0). The ground state \Q) of the magnet corresponds to 
the single trace operator of the M = 4 YM theory built from the scalar fields Z = ^43, 

\Q)=trZ L . (B.l) 

So all excitations propagate on this BMN background. The algebra of superconformal generators and 
their action on the vacuum state can be very concisely encoded by an oscillator formalism of Ref. [70] 
(see also Refs. [5 J and, in particular, Ref. [71] for conventions adopted herd^l and more details) such 
that the step-up operators e + acting on the nodes of the Dynkin diagram read 

4 = a\b\ , (B.2) 
4 = -C36 1 , 4 = c 3t c 4 , ej = Sb\ , 

in terms of bosonic (a a , a/ a , b@, ) and fermionic (ca,c^ a ) creation and annihilation operators corre- 
sponding to sl(2) ® Sl(2) and su(4) subalbegras of psu(2,2|4). Then a generic state of the spin chain 
reads 

\o Ki ,...,k 7 ) = f fl H +j ) Ki+i H-i) KA A (4) Ki P) (b.3) 

y=i,2,3 / 

where the vacuum state is 

\U) = (c 3t c 4t ) L |0) . (B.4) 

The quantum numbers of the resulting state can be easily determined by finding the eigenvalues of 
Cartan generators h p on them making use of the conventional commutation relations [h p , e^] = ±^4 p9 e^ 
with the Cartan matrix A pq . 

The hole excitations were already introduced in the main text. From the point of view of the Dynkin 
diagram the holes arise as "duals" to Bethe roots of the central 4-th node with all other excitation 
numbers set to zero Kj = 0, j / 4. This ABA corresponds to a long-range s[(2) spin chain with K4 = S 
magnons, represented by the light-cone covariant derivatives D + ~ D^, propagating on it, 

(e+) K4 \Sl) =\xB%Z L . (B.5) 

The resulting ABA equations were given in the main text and we quote them again below, 

1 _ A -Kfc)Q4Kfc-0 (B6) 

A+(ui )k )Q4(u4,k+i) ' 

Here we used the notation for the Bethe roots of the central node encoded in the Baxter polynomial 
Q±{u) = nj^i(' u — u 4,j)- This polynomial and Bethe roots were stripped from the subscript 4 in the 
body of the paper since they were redundant. 



'Compared to Ref. [71] we set Bethe-root numbers rij to Kj as in [5]. 
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Figure 8: Dynkin diagram for psu(2, 2|4) spin chain in the BFFBBFFB grading with excitation numbers 
Kj generated by the step-up operators et (|B.2|) acting on its nodes. 

Below we demonstrate how the remaining twist-one excitations emerge from the all-loop ABA equa- 
tions. To this end we have to know the coupling of nested Bethe roots u k ^4 with 114. This is easily 
deduced from the ABA equation for the central node [5] 

, A -( u 4,k) Q*( u 4,k -«) -rr c / \tl o t \ rat\ 

- 1 = A+{u + k)Q ^ k + l) \l Sutej,«u.) I! Srdur,^) . (B.7) 

where we set all excitation numbers for the left wing of the Dynkin diagram to zero, i.e, K\ = K2 = 
K3 = 0. The reason for this is that all twist-one excitations can be extracted from either half of Bethe 
roots (left or right). In the above equation equation, we introduced the scattering matrices of nested 
Bethe roots off the main momentum carrying root U4 

Su(u5,ih) = — — S 74 {u 7 ,u 4 ) = 1 — 9 J XjX ^ , (B.8) 
x 5 - x A 1 - g z /x 7 x i 



B.l Fermions 

To find the scattering matrix for propagation of twist-one fermionic excitations on the background of 
covariant derivatives, we have to include one of the two adjacent fermionic nodes, either 3 or 5. Let 
us focus on the latter. An analysis identical to the one that follows will yield the same results for the 
mirror-symmetric node 3 giving other twist-one fermions. The above identification can be easily verified 
by realizing the a fermion in the background of Zs corresponds to the following state 

e + e +lo\ - ihh^fJ ^Uf^ r AU L - 1 \n\ -t,,), 4 /" 



|fi) = [6 i 6l](a[c 4T )(c dT c 4T )^ i |0) = tr^fZ^ 1 , (B.9) 



where we ignored the overall diagonal Cartan generator fei&j acting on the state since it does no affect 
its quantum numbers. Thus if one sets all Kj = but for then one attains the autonomous sl(2|l) 
sector where in addition to the light-cone derivatives, the operators contain insertions of gauginos. The 
corresponding equations found from the nested psu(2,2|4) Bethe Ansatz then read 

1 A -Kfc)Q4K fc -0 fi c f . 

-1 = 7 , + ^ 777 \\S5A{U 5 j,UA,k) , (B.10) 

1 = Y[S 5i (u 5)k ,U4,j) ■ ( B - n ) 
i=i 

From the above identification of quantum numbers it is clear that the scattering of fermions off magnons 
is described 554, 

Su(u 5tk , U4,j) = S 5i (u^ k ,u 4 ,j) . (B.12) 
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Thus we see that the second equation of the above system of ABA equation for s[(2|l) is the Bethe-Yang 
equation for the phase factor 

Yi(u) = -\\S iA {u,u^ k ) (B.13) 

k=l 

accumulated by the fermion with rapidity u = u§ when it scatters off magnons, see Eq. (|3.29p . 
B.2 Gauge fields and bound states 

Let us finally turn to the gauge field. Again we can consider only half of the Dynkin diagram. Namely, 
the twist-one gauge field F + ±_ ~ Fn corresponds to the following set of oscillator excitation numbers 

K 4 = 2, K 5 = 2, K 6 = l, (B.14) 

such that 

e+(e+) 2 (e+) 2 |tt> = [(6^) W f )](4) V^^lO) = tlFnZ L-l (B15) 

So the scattering of the gauge field excitation on the background of K4 derivatives is described by the 
ABA equations 

A + (u+ fc ) Qi{u^ k + 1) x 5d - x~ k ' 

! = j| u h - u ^ -q x ^ - x _m t (B17) 

j = l u e,j u 5,k j =l x 5,k X 4 j 

-1 = n u "- ~_ U5,3 ■ (B.18) 

j=l U 6,k n 5,i 

One can immediately see that the trivial solutions to the last nested equation are 

1*5,1 = ub,i + h ' n 5,2 = «6,i - \ ■ (B.19) 

These are the stacks describing the gluon field found in the analysis of Ref. [67] . The scattering matrix 
of gauge fields off magnons can be found by fusion of elementary scattering matrices 554 on stacks 

S gf , 4 (n gf , u 4j ) = J] X5 ' k ~_^ , (B.20) 

k=l X5 > fc X 4j 

where we introduced the center-of-mass rapidity u g f = A short calculation immediately yields 

u + —ulfl — g 2 /x~x~ f 1 — g 2 /x~xZr 

= V/ + +1 V/ + • (B.21) 

u ~ n gf 1 ~ 9 l x +x g 1 - 9 / x+x g{ 

Finally, though the follow-up discussion goes beyond the scope of leading twist excitation, let us 
address for completeness of the picture to the bound states of gauge fields that are given by the operators 
(again for half of the Dynkin diagram only) 

D^F +1 _ . (B.22) 
They are encoded into the following quantum numbers 

K A = £ + 1 , K 5 = £ + 1 , K 6 =£, K 7 = £-l, (B.23) 



61 



reducing to the previously studied single gauge field for £ = 1. Solving the nested ABA equation one 
finds the stack structure of the nested Bethe roots which describe the above gluon bound states |67j 

u 5 ,j =u e + i(j-l-±£) , u 7>j =u e + i(j- \t) , (B.24) 

where ui stands for the center-of-mass of the bound state excitation (with U£ = i = u g f) related to uq 
via UQj = U£ + i{j — \ — \t). Fusing the S-matrices of elementary Bethe roots, one finds the one for 
scattering of gluon bound states off magnons 

^ x 5 , k - x 4J ^ 1 - g 2 /x 7 , k x 4j 
Computing the latter yields the result quoted in Eq. (|3.38p of the main body of the paper, 

- x 4 • 1 - g 2 /x l ' Pi x 4j 
= xM-xil-f/xMxi ' (R26) 

where x 1 ^ = x{u ± i£). 



C Mirror F-functions at weak coupling 

In this appendix, we compile various results concerning the amplitude of propagation of excitations in 
the vacuum in the mirror kinematics. As we advocated in the main text, the latter admits the following 
representation 



n mirror w = (-i) 5 ex P 



2E*(u) log 7? - 25E*(u) + 0(1/7/' 



2^ 



(C.l) 



where S stands for the spin, E+(v) is the energy of the excitation with the ^-flavor, and fj = Se* 
at leading order in the large spin. Both E+(u) and 5E*{u) are solely functions of the coupling g and 
rapidity u, and, in particular, they are spin-independent. Our goal is to derive explicit weak coupling 
expressions for (fCTTil for all twist -one excitations, with, in particular, the fermion in the large rapidity 
domain (* = If). The energies E+^u), with * = h, If, gf, are known and their weak coupling expansion 
was computed in Refs. |48l I44j . Thus, we have to complement this data with power-law series expansion 
in 't Hooft coupling g 2 for the corresponding 5E*(u) making use the representations determined in Eqs. 
(|3.63l 13.691 I3.89|) . We shall present our results, however, in a form that is different from (1C.1|) but 
more convenient for the evaluation of the Liischer formula at weak coupling. It would not be difficult 
to extract from them the two components E+(u) and 5E+(u) if needed. Note also that to the accuracy 
in the large-spin expansion, the 0(l/rj 2 ) corrections in (jC.ip are irrelevant for our current analysis. 

In order to avoid redundancy in the representation of the Y-functions for the various twist-two 
excitations, it is instructive to cast them in the following factorized form 

ymirror (u) = ( _ l} 5 ^^^mirror^ (C2) 

where f*(u) is a i-periodic function of u, while ?/™ irror (it) is a reduced Y- functions, henceforth called the 
essential factor for the reason that once it is set to one for all excitations, the Liischer formula (|3.92p 
vanishes identically. By definition, the periodic factor satisfies the functional relation 

U(u + i) = h{u) , (C.3) 

for any of the twist-one excitations. Here the shift u — > u + i is performed outside of the strip —2g< 
5ie(u) < 2g, where the function /*(«) happens to have branch-point singularities. One of the reasons to 
factor out the latter quantity in Y-functions (|C.2p is that it possesses the following interesting property 

/h(u) = / g f(u) = -/if(u±§). (C.4) 
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Here, we emphasize again, the fermionic expression is obtained for the large fermion with |x| > g. The 
identity (|C.4p . which can be immediately verified to hold true for leading contribution at weak coupling 
making use of the explicit results 

= /gf(«) = 2 *ff (1 + 0(g 2 )) , h(u) = (1 + 0( 5 2 )) , (C.5) 

r] z cosh (7tmJ rj z sum (7ru) 

is responsible for numerous cancellations between bosonic and fermionic contributions in the Liischer 
formula. When evaluating Liischer corrections up to five loops (= 0(g 10 )) the established conspiracy 
between different flavors makes it sufficient to know one order lower, i.e., up to 0(g s ). This 

property also relies on the fact that y™ irror (^) = l + 0(g 2 ) at weak coupling, as we shall see momentarily. 
Moreover, thanks to the relations (|C4p . it is enough to compute the case of the scalar excitation only. 
On the other hand, the factor y™ rmT (u) has to be determined independently for all twist-one excitations 
anew up to 0(g 6 ). Explicit derivations and final expressions both for y™ irror (u) and fh(u) can be found 
below. Finally, we shall also include, at the end of this appendix, the weak coupling expansion for 
the momenta of various excitations, since they enter the Liischer formula on equal footings with the 
Y- factors. 

C.l Periodic factor 

Let us start with the calculation of the periodic factor shared by all ^-functions. As we stressed above, 
we can restrict our consideration to the case of the hole, since the gauge field and (large) fermion 
excitations follow from Eq. (|C.4p . We split the complete expression (|3.63p as = — log — log y™ irror 
and introduce the first term in the sum as 



, 01 , 2 M o , a f°° dt ( cos ^ et/2 - J ^gt)) Jo&gt) - 1/2 

log AW = 21 og {g h) ~ 2^(1) +4 / — 

Jo 1 e — 1 

_ 2 f°° dt cos (ut) e t / 2 -J (2gt) _ [™ dt ^(2gt)J (2gt) - 2 gil t 

Jo t e*-l 1+[ 9 ' J t e*-l 



(C.6) 



It can be equivalently rewritten as 

ng 2 \ f°° dt cos (ut) e*/ 2 



log/ h (u) = 2 log ( ^ n -^) +4 T 

V??cosh(7r'u)y J 



t e*-l 



J (2gt) - 1 - h+{2gt) 



+ a, (C.7) 



where a is constant w.r.t. u, 



a = -4 r dt J o( 2 ^- 1 + 2 I" dt M2gth(-2gt)+2 gil t 

Jo t e*-l Jo t e*-l ' 1 

and we used the following explicit solution for the integral 



CO 







dt cos (ut) e*/ 2 -1 -t/2 _ f tt 
t e* — 1 ^ 1 cosh (ttu) 



2 I - ot i -=l°g( IZJTTZ^ (C9) 



To get the representation (1C.8P for a we used the relation ^(—2gt) = r y + (—2gt) + 7_(— 2gt) = r y + (2gt) — 
7 _ (2gt). 

Before we move to the analysis of the weak coupling expansion, let us verify first that it has the 
periodicity (|C,3p . To demonstrate this, we notice that making use of the decomposition (|2.8p . the 
expression in the square brackets in Eq. (|C.7p . can be rewritten as a series over the even Bessel functions, 

J (2gt) - 1 - \ 1+ (2gt) = - J^(2 + 2n l2n )J 2n (2gt) . (CIO) 

n>l 
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Then we observe that each integral contributing to the series 

dt cos (ut) e*/ 2 



o 



, , J 2n (2gt), n = l,2,..., (C.ll) 
t e r —1 

is i-periodic. This is evident by shifting the argument by half the period ±| 

r f ° s (( - f *»> e ' /2 j^ gt ) = r f^j^) + 1 r±*~j»M, 

Jo t e —i Jq t e — i z Jo t 

and recognizing that the last term 

f°° dt 

/ -e^J 2n (2 5 t) (C.13) 
Jo t 

is merely two representations of the same function, namely (ig/x) 2n /(2n), for « taken in the lower- 
and upper-half plane, respectively, as was exhibited by Eq. ()2.10p . Thus the right-hand sides of the 
above equations coincide independently of whether the shift is positive or negativj^l. More precisely, 
the two above integrals are related to one another through analytic continuation for u away from the 
cut u 2 < (2g) 2 . So the periodicity of (jC. llj) is guaranteed if the shift u — > u + i is performed outside 
of the strip —2g < 3?e(u) < 2g, where the integrals (jC.lip have an infinite number of cuts running 
between u = ±2g + in/2 with n = ±1, ±2, .... One may worry about the commutativity of operations 
of summation and integration. However, the summation does not move away or change the nature of 
existing singularities, but it induces extra ones though they are confined to the cuts in the strip. 

Let us come back now to the perturbative expansion of f|C.7|) . Note first that the second and last 
terms in the right-hand side of (1C.7|) appear explicitly suppressed at weak coupling compared to the 
leading term. Thus, dominant contribution for small 't Hooft coupling reads 

log AO) = 2 log ( ^r^] + 0(g 2 ) . (C.15) 

\ 77 cosh (ttu) ) 

As we previously explained, we only need to know log/h(n) up to order 0(g A ), i.e., fh(u) up to 0(g 8 ). 
This result can be easily obtained by Taylor expanding the integrand in (jC.lOP and then plugging it 
back into Eq. (j(T7j) . We then find that 



f^ U > = ~2 ^27 V eXP 

r/ z cosh {ttu) 



tt 2 7r 4 (cosh {2iru) - 2) ^. 



(C.16) 



with the coefficients being 



/3 = -2 5 2 (l + 72 ), 7 = -y(3 + 472-274), (C.17) 



and where we used the value for the integral 

f°° , te* 

Jo e 1 coslr (ttu 



00 te*/ 2 cos (ut) tt 2 



dt t ; ' = — — — - , (C.18) 



48 Note that there is also an alternative representation for the integral (IC.11|) 

M2 



f°° dt cos (ut) e i\«V^( _1 ) r 1 a 1 \„-2Trku 

/ T ^TTT ^2n(2fft) = (-1) 2^^r~ l2n ( 4nk 9) e > ( C - 14 ) 



2k 

k>l 



with Iin{z) = (— l) n J2n(«z) being the modified Bessel's function. The representation is valid for 5Re(u) > 2g and obtained 
by closing the integration contour at 00 and performing integration by the method of residues. The expression in the 
symmetric domain 5Re(it) < —2g can be easily obtained by the substitution u — > —u. The periodicity then becomes 
manifest. 
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and its derivatives. Finally, to determine the constant a, one can Taylor expand the numerators of each 
integrand in (|C.8p . and get 

27T 2 27T 4 

a = —g 2 (2 + 7a) + 4C 3 ff 3 (3 7 i - 7s) - -r=-g 4 (9 + 872 - 74) + • • • • (C.19) 
6 45 

To find the numerical values for coefficients accompanying powers of the 't Hooft coupling for these it- 
independent constant, we use the known decomposition 7 n = 7^ log 77 + ^7^ valid up to l/Ty 2 -suppressed 
corrections and substitute the available expressions (12.641 l2~65j) for the coefficients 7^, to obtain 

&TT 2 1 S7T 4 

a = 24C 35 4 log fj + ^-g 2 - ^-g 4 + 0(g 6 ) , (C.20) 

and 

P = -2g 2 - ^-g 4 + 0(g 6 ) , 7 = ~g A + 0(g e ) . (C.21) 

C.2 Essential factor 

Now we move on to discuss the second ingredient of the ^-functions, the essential factor y mirror j n the 
decomposition ()C,2p . The latter reads 

bgtf*~(«) = 2 r ^-^,-(2gt) C Q S{ ut)-2 911 t 
Jo t e _1 

for the excitations of ^-flavor encoded in the conformal spin s+, i.e., s+ = 1/2,1,3/2, for the scalar, 
fermion and gauge field, respectively. 

As a consistency check, let us verify that the essential factor (|C.22p put together with the periodic 
one (|C.6p reproduce the complete expression for ^-function in the mirror kinematics according to (|C.2p 
with (|C.4p . This is quite straightforward for the hole where the sum of (10221) and flUl| immediately 
reproduces the total expression (|3.63|) . The same holds for the (large) fermion, where the shift involved 
in (|C.4p can be performed making use of the representation f|C. 12|) . Finally, for the gauge field, the 
matching is recovered with the help of the identity 

/ -J (2gt) 1+ {2gt) =0, (C.23) 
Jo 1 

which can be understood as emerging from the orthogonality condition for the Bessel functions 

/ -Jo{2gt)J 2n {2gt) = 0, n = l,2,..., (C.24) 
Jo t 

in the Neumann series representation for 7+(i), i.e., 7 + (2gi) = J2n>i ^(^ n )l2nJ2n(2gt). It is assumed 
that the commutation of the summation and integration is legitimate due to tamed large-n asymptotic 
behavior of the expansion coefficients 72 n - 

At leading order in weak coupling, the essential part (|0.22|) is controlled by the one-loop energy of 
excitations. Indeed, this follows from 7_(2gi) ~ Ag 2 tlogf], together with 2g^/i ~ Ag 2 logf], which leads 
to 



yf rro » = exp 



4 5 2 (VKs* + iu) + ^(s* - iu) - 2^(1)) log fj + 0(g 4 ) 



(C.25) 



with ip{z) the digamma function. It is straightforward to determine further terms in the perturbative 
expansion. Namely, with the help of the series representation for 7„(2gt) in odd Bessel functions, 

7 _(2^) = 2^(2n - l) 72n _ 1 J 2n _ 1 (2 5 i) , (C.26) 

n>l 
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and the scaling behavior J2 n -i{2gt) ~ (gt) 2n 1 at weak coupling (with t fixed), one can truncate the 
Neumann series (|C.26p and find the representation 



2/r rt » = exp 



(C.27) 



aip (s*,u) + bip 2 (s+, u) + cip A (s*,u) + 0(g ) . 

Here we introduced the notation 

VVi(s*, u) = ipn(s* + iu) + ip n (s* - iu) - 2^(l)£ n)0 , (C.28) 
with polygamma functions ip n (z) = <9™ +1 log T(z). The coefficients a,6, c, in (|C.27|) are given by 

a = -2gyx, ^ = 5 3 ( 7l - 73) , c = _^_(2 7l - 3 73 + 75) • (C.29) 

This result is obtained by Taylor expanding (|C.26|) up to 0(i 5 ), plugging it back into (|C22|) . and 
evaluating the resulting integrals with the help of the relation 

f°° t^e-i*-^ cos (ut)-5 nfi t 

2/ dt — '— = -t/j 2 n{s, u) , (C.30) 

Jo e — 1 

with the right-hand side being (|C.28p . 

Finally, the weak coupling expansion of the coefficient (IC.29j) is found by substituting the expressions 
obtained before for 7„ = 7 log ry + £7^, neglecting 0{l/rf ) corrections. Using (|2.64| [2765]) . we get 

a = - ( V - + ^-A log fj + 12C35 4 - + 8OC5) 9° + 0(g 8 ) , (C.31) 



b = ( 2g A - ^g e ) logry - ^g & + 0(g 8 ) , c = logry + 0(g 8 ) , (C.32) 



and 

b = (2^ - '- , , 

3 J ' 3 v / 3 

where we also used the fact that 75 ~ when evaluating c. 
C.3 Momenta 

The final ingredient in the Liischer formula (|3.92p is the momentum p+ of the ^-excitation. Analogously 
to the previous consideration, we can decompose these into periodic and essential components as well, 

p,(u)= P r(u)+pT(u), (C.33) 

with the periodicity condition p+ cr (u + i) = p+ cr (u) for u outside the strip — 2g < $ie(u) < 2g. Then, 

similarly to (IC.4h . we have the following relations between the periodic components of the twist-one 
excitations, 

PTO = , *T («) = pT{u ± §) . (C.34) 

As alluded to before, a purely periodic integrand does not contribute to the Liischer formula, thus we 
merely need to know p^ cr (u) up to 0(g A ) only. On the other hand, the essential part pl ss (u) is needed 
up to 0(g & ). 

The expression for the periodic component of the scalar momentum is given by 

per,, f°° dt sin (u^e^^gt) 

Ph («)=-y o T — — • ( C - 35 ) 

Note that it involves the solution to the BES equation only, i.e. 7 , as opposed to full function 7 = 
7 logry + ^7 + 0(l/ry 2 ) entering the Y-functions, see Eq. (|C.22p for instance. So this portion of the 
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momentum depends on the rapidity u and 't Hooft coupling g but not on the spin of the excitation. 
The weak coupling expansion of (|C35p is constructed along the same lines as before. Below, we present 
the result for the derivative of the momentum instead, since this is the quantity that enters the Liischer 
formula. The former reads 



M!M = -2vrV * (l ~ ~ 2 7 rV C ° Sh(2 r ) " 2 U ~ ^ ) + 

du cosh (ttu) V 3 / cosh (ttu) \ 3 / 

The essential part of the momentum admits the following general expression 

eft sin (ut)j1(2gt) 



0(g b ) . (C.36) 



pT (u) = 2u- f^W^e-Vi)*, (C.37) 

and explicitly depends on the conformal spin of the ^-excitation. At weak coupling it is simply given 
by 

= 2 + 4&g«ih(**,») + 0(g 8 ) , (C.38) 
with the notation for the combination of the polygamma functions introduced in Eq. (|C.28p . 



D Liischer formula at weak coupling 

In this appendix, we present details on the computation of the perturbative expansion of the Liischer 
formula by assembling all results derived in the previous appendices. Starting from the expression for 
the Y-functions in the mirror kinematics given in Appendix [Cj it is quite straightforward to perform 
the weak coupling expansion of Eq. (|3.92p . In Sect. 13.31 we found that it admits the following general 
representation 

E ^ = t^j2v n g 2n + o(i/ v 2 ), (D.i) 

^ n>2 

where V n are polynomials of degree n — 2 in the variable log fj = log r\ + 7e . In the main text we provided 
explicit expressions for each of the polynomials V n up to five loops, i.e., for n = 2, . . . , 5, see Eq. (|3. 106[) . 
Below we provide intermediate expression obtained with the help of the mirror y-functions listed in 
Appendix O as well as evaluation of all relevant integrals. 

Making use of the results in Appendix [Cl we find the following expressions for the polynomials V n , 

v 2 = o, 

V 3 = 4Z^ log fj, 

V, = -12C 3 X (1) + (-2£<°> - 12 jW + Att 2 !^) log?? - 81^ log 2 fj , 



Pb = (^£ (0) + 36C 3 J W + (S0C 5 - ^) Z (1) ) (D.2) 



+ 



-C + 6M+ 32Af - 80 - 2^ 2 /C (0) - \2t?J^ + — Z« + 48C 3 X {2) J log fj 
3 5 / 



+ (sKW + 24i7 (2) + 96C 3 X (D _ ^jpA lQg 2 - + 3| J(3 ) bg 3 - 

in terms of the integrals X^ n \ J^ n \ . . . presented below. The latters encode contributions both from 
bosonic and fermionic excitations. The representation (|D.2p already accommodates simplifications ob- 
tained after discarding all terms involving purely periodic integrands, since they cancel out when adding 
bosons and fermions together. The result (|3.106p quoted in Sec. 13.31 follows immediately from Eq. (1D.2P 
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upon substitutions of explicit expressions (QTT31 HXTBl ITXT91 HX2T1 ITX231 11X251 HX271) for the integrals 
X( n ), J'^ n \ .... The general structure of the latters reads 



where summation runs over all twist-one excitations * = h, gf, f rn Here the measure factor M*(u) is 
expressed in terms of the hyperbolic functions, e.g., 

M h , gf (u) — , M f (u) ~ . 2 - , (D.4) 

cosh (7ra) sinn (7rrtJ 

which are suppressed at least as exp (— 27r|tt|) at large rapidity. The combination of polygamma functions 
t/)n(s±,u) was introduced in Eq. (|C.28|) where s* stands for the conformal spin of the ^-excitation, 
s = 1/2,1,3/2, for the hole, fermion and gauge field, respectively. Finally, all integrals in (|D.3|) are 
taken over the entire real axis, with the =LiO prescription for the integration around u = in the 
fermionic integrand^ With the latter specification, all integrals of the type (|D.3p that we are going to 
consider, are (separately) well defined. 

Before listing all relevant integrals and their explicit values, let us briefly comment on the procedure 
of their computation. Most of the integrals can be evaluated using the Fourier representation for 
ip n (s±,u) given in (|C.30p , This immediately applies for k = 1 in (|D.3|h However when the integrand 
displays higher powers of k, we did not find a better way of computing the integrals but by closing the 
integration contour at infinity and then performing the resulting integrals by the method of residues. 
This procedure typically yields an infinite series representation for them that can be calculated for low 
values of k. Notice however, that for the type-X integrals, it is necessary to first use the integration 
by parts in order to make the sum over the residues convergent. This reflects the fact that these 
integrals involve the digamma function ip(s + iu) which possesses growing asymptotics at large rapidity, 
ip(s + iu) ~ log u. Finally, it is interesting to notice the similarity of these integrals to the ones 
appearing in the analysis of the OPE for light-like Wilson loops |444 \$9\ [68] and thus one could adopt 
the same technique exploring their motivic structure and finding corresponding symbols as a systematic 
framework for their calculation (see [72] for explicit results). 

D.l Integrals for leading logs 

The only integrals that we need for evaluating the leading log behavior are of the type 



(^(5 + iu) 


' + ^(5 - iu) 


-2^(1))" 




cosh 2 (iru) 




0K§ + »«) 


1 +-0(1 ~iu) 


-2VU)) fc 




cosh 2 (iru) 






+ ip(l — iu) 


-2^(l)) fc 


sinh 2 (iru) 



T£> = ir I tin 
Xf ] =ir J du 

where k is an positive integer, which runs over k = 0, 1, 2, 3, if one aims at computing all the way up to 
five loops. They are obviously related to the integrals f|3. 103|) . introduced in Sect. 13. 3[ by 

li k) = ^Zi fc) • (D.6) 
k\ 



49 The fermionic integrands arising here are always implied to take the form of the ones in the large rapidity domain. 
Nevertheless, in order to simplify the notation, we label all the former with 'f only, instead of 'If. 

50 We recall that the integration contour for the fermion passes above, or equivalently below, the cut u 2 < (2g) 2 . At 
weak coupling, it implies that one should avoid all the singularities of the type ~ g 2n /u 2n which accumulate at u ~ 0. 
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For the lowest three ks, they read 

T (o) _ 9 T (o) _ 9 T (o) _ _ 9 

4 1, = -4, ^ = - 4 + ¥. &-*-t- < D ' 7 > 

Z?>=16-^, 2g> = 16 + ^-16&, 2?> = -16+^+80,. 

The k = 3 case, which is required for the computation of the five-loop expression, is much more 
involved. As we are only interested in a particular weighted sum of different flavors, we circumvented 
the complication by computing suitable linear combinations of the integrals (|D.5|h For instance, by 
applying the identity 

V>(§ + iu) + VK| - iu) = ^(5 + iu) + if){\ - iu) + 2 1 1 , (D.8) 

one can relate gauge-field integrals to scalar ones up to integrals involving lower powers of digamma 
functions. The latter integral can then be taken by the method of residues and yields 

X (3)_ X (3) =47r2 + 7^_ 48C3 _ (DQ) 

Similar manipulations can be done for the fermionic integral. Namely, by deforming the contour of 
integration, u — > u — i/2, such that 



1 1 



(D.10) 



sinh (ir(u — |)) cosh (nu) 
and using furthermore the identity 

^(1 + i(u - I)) + ^(1 - »(« - I)) = i>{\ + *«) + ^(3 - i«) . ( D - n ) 

u — g 

one finds along these lines 

X f (3) +4 3) = -^- CD.12) 

It turns out that the two expressions (|D.9p and (|D.12|) are all we need for our purposes. Summarizing 
our findings, we can compute the weighted sums immediately from (jTX7l) . (iDlJl) and (11X12]) . They read 

0, 

_4vr 2 
" 3 ' 

o 9 (D.13) 
(2) -^+32 C 3, 

Svr 2 - 96C 3 . 

Sd 15 

Plugging these values into (|D.2h . one immediately reproduces the result (I3.106P obtained in Sect, 
for the leading logs. 



J(°) 


^ 6 Jj 0) 




+ ! 


iX f (0) 


X (l) 






+ 1 


sx f (1) 


X (2) 


= 62f 




+ 1 




X (3) 


= e4 3) 


-1- 21^ 
+ zx gf 


+ 1 


ST f (3) 
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D.2 Integrals for subleading logs 

Turning to the subleading logs, we find two kinds of integrals. The first class has the following repre- 
sentatives 

du 



J h 


-■7 




-■7 




-■7 



(^(3 + iu) + - iu) 


-mi)) k 


cosh 4 (ttu) 




(-0(1 + iu) + -0(f -iu) 


-2^(i)) fc 


COsh 4 (7TU) 




(0(1 + iu) + 0(1 - iu) - 


-20(1))* 


sinh 4 (7ru) 



^ 1^2 1 ia > Z. Lll ~ ; ~ f^Vj^ (D-14) 



where k is a positive integer, running over k = 0, 1, 2 for three, four and five loops, respectively. Their 
values are 

4?r 2 (o) _ 47^ 

3 ' Jf " 3 ' 

8vr 2 22vr 4 _ m 8vr 2 llvr 4 



T (0) _ 


4vr 2 
3 ' 






T-(l) _ 








T (2) _ 




227T 4 

45 


T (2) 


J f (2) = 


32tt 2 
3 +16C3 " 


227T 4 

45 


16vr 2 C 3 
3 



(1 ) = _ 4Cs + 



3 45 ' 1 3 45 ' 

32vr 2 , 22vr 4 32vr 2 C 3 

= + I6C3 + — - 24C 5 , 

3 s 45 3 

- 12C5 • 

(D.15) 

In the present case the weighted sums have an explicit minus sign for the fermion contribution, yielding 



(D.16) 



J(0) 


= 6J® 


+2 4 0) -, 


s.7 f {0) 


= 0, 




Jil) 




+ 2J g ( f 1) -; 


3J f (1) 


44vr 4 
45 ' 




J(2) 




+ 2J" g { f 2) - ; 




887T 4 

45 





The second class of integrals is of the form 
/Cj^ = it J du 



_,L. 

f 7 sinh 2 (7ru) 

Here we need only A; = 0, 1 at four and five loop order, respectively. We found by an explicit computation 

4 0) = -24Cs , = -24 C 3 + ^ , Kf = 24 C 3 - ^ 



(02(5 + iu) 


> + Ml 


-iu)) (0(± 


+ in) + 0(7; — iu) 


-2y,(i)) fe 






cosh 2 ( 


ttu) 




(v> 2 (§ + iu) 




- *«)) ty(§ 


+ iu) + 0(| — iu) 


-2y>(i)) fe 






cosh 2 ( 


ttu) 




(02 (1 + iu) 


+ 02(1- 


- iu)) - 


\- iu) + 0(1 — iu) - 


- 20(1))* 



^( fc ) =n J du -t- "'J i- '.'212 ^ ^ ^ '"1 + ' - /;; ' - ( [) 1 7 ) 

/C**^ = it I du 



-gf - -"S3 "15" ' f ~ S3 ~ 15" ' 
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such that we weighted sums read 



X;(0)^ 6X; (0) +2X; (0) +8X; (0) = _8vr 

15 (D 19) 

= 6Xf +2^) + 8*f = ^ + 32(3 + ^ + ^ + 176 C5 . 



D.3 Other types of integrals 

Last but not least, for the five-loop computation of the next-to-next-to- leading logs one requires com- 
puting of various miscellaneous integrals. We list them below with their final values. First, 



£ h = ir / du * -g- — = -480C 5 , 

J cosh (iru) 

^4(5 + *«) + ^4(5 - »«) , on . 32vr 6 



£gf = vr / ^ -2- — = -480C 5 + -Tz- , ( D -20) 

J COSh (7TU) Drf 

/" , ^4(1 + ire) + ibAl - i«) , , 16vr 6 
Ct= ir I du — , \ 2 \ = 480C 5 



sinh 2 (iru) ' 63 

with the sum equal to 



647T° 

£ = 6£ h + 2£ gf + 8£f = — — . (D.21) 



Next, 

M h =ir 3 J du - - ; ; , — - - IG7r\ :j - 

.3 



COSh 4 (7TU) 



M^^iiu ^(1 + '">+ fell''") _ _ 16A> _ 80(5 + ^ , (D . 22) 



cosh 4 (7rn) 945 

M{= tt 3 J du ^ 2(14 y; ; : — — - 1 6-^ :i - 8o<, - 



^ 2 (l + m)+V> 2 (l-m) 2 , 124tt 6 

■ ,4/ \ = _167r C3 - 8OC5 + — 

smh (7ra) 945 



with the weighted superposition where fermion are accompanied by an explicit minus sign given by 

4967T 6 

M=6M h + 2M Ki -8M{ = — — — . (D.23) 
6 945 

Third, 

5 f (^(\+iu)+^(\-iu)-2^(l)) 327T 4 

A/ h = 7T du -g- = -47T C3 - -TT" - 4 ^5 . 

7 COSh (7Tll) 15 

7 cosh D (7ru) 15 945 

5 y (y,(l + in) + ^(l-in)-2^(l)) 2 32vr 4 191vr 6 

7V f = vr / — g- — = 4tt C3 + -T3- + 4C 5 , 

7 smh D (iru) 15 945 

resulting in 

7647T 6 

A/" = 67Vh + 27V g f + 87Vf = — — . (D.25) 
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Finally, we have the set of integrals 



5 f. COSh (2TTU) (V>(i + iu) + ^ ~ iu) ~ 2^(1)) 167T 4 2 

c h = 7T / aw nrr — : = — = — ^ 0$ + 4Cs , 

j cosh (iru) 5 

5 /• cosh (2iru) (V>(f + ra) + - ra) - 2^(1)) 16vr 4 2 542vr 6 

O gf = 7T / du -g- = Air (3 + 4C5 + ~FT7r~ ' ( D - 26 ) 

y cosh (iru) 5 945 

_ 5 /• , cosh (27TU) + in) + - m) - 2^(1)) 16vr 4 2 271vr 6 

C f = 7T I d« — g- — = 4vr C3 + 4C5 + , 

J smh D (iru) 5 945 

yielding (again with a minus sign for the fermion in this composition) 

10847T 6 

= 60 h + 20 gi -80i= — . (D.27) 
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